KAZAKH NATIONAL UNIVERSITY named after AL-FARABI 



Kenes B.Jakupov 



CORRECTION OF CONTINUUM 

MECHANICS THEORETICAL 

PARADOXES 

(I) 



Almaty 
2009 



UDC 532.516 

BBK 22.5 
D 40 



Recommended for publishing 
RSE Institute of mathematics, informatics and mechanics 



Reviewers: 
Doctor of Physical and Mathematical Sciences, Academician of NAC RK 
T.Sh.Kalmenov; 
Doctor of Physical and Mathematical Sciences, professor M.T.Zhienaliev 



Kenes B. Jakupov 

D 40 Correction of continuum mechanics theoretical paradoxes - 
Almaty: publishing house Tiinorpaqbiifl «K2», 2009 - 270p. 

ISBN 978-601-0288-5 

The monograph describes paradoxes of Stokes and Navier-Cauchy-Lame 
hypotheses. New equiations of viscous fluids dynamics, elasticity theory are 
theoretically justified. Current technologies of building effective numerical 
techniques for combined differential equation systems in true variables "velocity- 
pressure" are given in details. 

The monograph is recommended for students, candidates for master's degree, 
post-graduate students, persons working for doctor's degree and researchers 
specialized in the sphere of fluid mechanics and elasticity theory. 



J\ 



1603040000 
00(05)-09 

BBK 22.5 



ISBN 978-601-0288-5 © Kenes B. Jakupov 

© Publishing house Trniorpaqbiia «K2» 

2 



LIST OF CONTENTS 

FOREWORD 7 

Chapter 1. FALLACY OF NAVIER-STOKES EQUIATIONS. 

NEW EQUIATIONS OF VISCOUS GAS DYNAMICS 9 

§1. Paradoxes of Euler, Leibniz formulae 10 

§2. Paradoxes of defirmative motion of medium elementary 

volume 14 

§3. Paradoxes of cumulative derivation of continuum dynamics 

equiations 20 

§4. Inductive approach 24 

§5. Paradoxes of the first Helmholtz theorem 29 

§6. Paradoxes of Stokes hypothesis 35 

§7. Newton's nonsymmetric strain tensor. Paradoxes of determination 

of viscous normal stress are related to hypotheses on pressure and 

to Pascal law for ideal liquids 43 

§8. Background for fallacious conclusion on strain tensor symmetry 50 

§9. Momentum theorem does not mean symmetry of strain tensor. 

Strain tensor is nonsymmetric 55 

§ 1 0. Paradoxal application of theorems on angular momentum 

(momentum) variation. Fallacy of Navier-Stokes equiations 64 

§11. Analogue of Stokes hypothesis. Antisymmetric tensor strains 68 

§12. New equations of vicous fluids dynamics with Newton 
nonsymmetric strain tensor 

dv. 
ff MH) =-[p + (ju/3-Ju')diw]S ji +ju— '- 69 

§13. New equations of vicous fluids dynamics with Newton 's 
nonsymmetric strain tensor in cylindrical coordinates 

n H = -[p + (1 / 3ju- ju')divv]E + /jS 73 

§ 14. New equations of vicous fluids dynamics with Newton 's 
nonsymmetric strain tensor in spherical coordinates 

n H = -[p + (\/3jU- ju')divv]E + juS 74 

Chapter 2. FALLACY OF PRANDTL EQUATIONS 

IN BOUNDARY-LAYER THEORY 77 

§ 1 . Velocity vector components are determined with the help of 

Newton 's second law or momentum theorem 77 

§2. Action of law of mass conservation upon flow is performed via 

pressure only 78 

§3. Fallacy of Prandtl equations is in incomplete observance of 
second Newton law and due to absence of action of law of 



conservation of mass upon flow via pressure 81 

§4. Laminar boundary layers in compressible flow 88 

§5. Three-dimensional boundary layers 89 

§6. Substitute of Blasius equation 90 

Chapter 3. ABOUT TURBULENT FLOWS SIMULATION 91 

§ 1 . Paradox of Reynolds averaging formula 92 

§2. Approximate calculation of averaging integrals 93 

§3. Application of pulsations equiation for closing of Reynolds 

equations 95 

§4. Turbulent flows simulation with account to friction caused by 

normal strains 99 

Chapter 4. FALLACY OF APPLYING "DARCY law" 

IN FILTERING THEORY 107 

§1. How "Darcy law" was created for multidimensional filtration 108 

§2. On setting boundary conditions for unknown functions u,v,w,p 110 

§3. Boundary conditions of Dirichlet and von Neumann 110 

§4. Boundary conditions of von Neumann type Ill 

§5. Paradoxes of "Darcy law" (7), (8) 113 

§6. Application of "Darcy law" contradicts to energy conservation 

law 115 

§7. In "Darcy law", only in one-dimensional steady-state filtering all 

three laws of physics are observed, (/equation of continuity 

accounts for mass source 121 

§8. Forzheimer model of filtration in saturated porous media 123 

§9. Forzheimer model correction 123 

§10. Numerov model of filtration in saturated porous media 124 

§11. Interest to "Darcy law" is related to applications of complex 

variable theory 1 24 

Chapter 5. FALLACY OF THE "FICTITIOUS DOMAIN 

METHOD" 125 

Chapter 6. CONSTRUCTIONS OF MONOTONOUS SIMILAR 

SCHEMES 136 

§ 1 . Monotonous schemes constructing technology 136 

§2. Monotonous similar schemes 2 order of approximation of. 140 

convection terms in equations with constant coefficient of molecular 

transport 140 

§3. Monotonous similar schemes 2 order of approximation 

of convective terms in equations with variable coefficient of 

molecular transport 143 

Chapter 7. PARADOXES OF DEFORMATIONS THEORY AND 
EQUATIONS OF NAV1ER-CAUCHY-LAME 



ELASTICITY THEORY 146 

§ 1 . Paradoxes of deformations theory 146 

§2. Alternative presentation of relative motion d\x 149 

§3 . Paradoxes of Navier- Cauchy-Lame hypothesis 151 

§4. Paradoxes ofHooke's law for Navier-Cauchy-Lame 

symmetric strain tensor 157 

§5. Equations of elasticity theory for nonsymmetric strain tensor 159 

§6. Explicit scheme 1 60 

§7. Solution of three-dimensional problems of elasticity theory on 

uniform grid 162 

Chapter 8. FALLACY OF APPLYING STAGGERED GRIDS IN 

MUNERICAL CALCULATIONS 166 

§ 1 . Irreducible errors in schemes built upon staggered grids 166 

§2. On additional boundary conditions and on the problem of one-valued 

solvability 171 

Chapter 9. TECHNOLOGIES OF EFFECTIVE SCHEMES 

BUILDING ON A SINGLE GRID 172 

§1. Formulation of the problem of longitudinal flow past a plate 172 

§2. Technologies of building difference schemes in viscous fluid 

dynamic equations on a single grid. Semi-implicit difference 

scheme 1 76 

§3. Semi-implicit scheme realization algorithm 180 

§4. Technology of obtaining difference equations for pressure 182 

§5. Iteration algorithm. Global iterations method 189 

§6. Programs on TURBO-PASCAL 199 

§7. Approximation of convection members without 

"schemati viscosity" 206 

§8. Semi-implicit schemes for solving of three-dimensional equiations of 

viscous incompressible fluid 210 

Chapter 10. PARADOX OF "NEGATIVE" GAS DENSITY. 
EFFECTIVE METHODS OF SOLVING 

COMPRESSIBLE GAS EQUATIONS 219 

§ 1 . Paradox of negative gas density is caused by incorrect realization of 

law of mass conservation 219 

§ 2. Semi-implicit differential scheme of complete equations of 

compressible viscid gas 22 1 

1 . Scheme of continuity equation 224 

2°. Realization of state equation p = pRT 225 

3 . Difference equation for pressure 225 

§3. Global iterative method of calculation of pressure 228 

§4. Indexed family semi-implicit schemes 231 

5 



§5. Technology building of scheme for march equations of compressible 

perfect gas 234 

§6. Method of determinants in explicit iterationless scheme of 

solvingstationary march equations 236 

Chapter 11. UNIVERSAL METHOD OF ACCELERATION OF 

CONVERGENCE OF ITERATION ALGORITMS 240 

Chapter 12. ABOUT APPLYING OF SCHEME OF PREDIKTOR- 

CORRECTOR IN H YDRODINAMICS 249 




Residual concentrations of impurities (marked in red), demolished in the 
flow of air over the four buildings of varying heights. Wind flow is 
incident from the left side of buildings. 



To the memory of professor Kh.I.Ibrashev 
and Academician N.N. Yanenko 
dedicated is this book 



FOREWORD 
The book was written on the basis of given by the author for many 
years lectures on "Continuum mechanics" and "Computational mechanics" 
at the faculty of mechanics and mathematics of the Kazakh National 
University named after Al-Farabi. Content of this book reflects up-to-date 
critical perception of a number of established fundamental principles and 
basic concepts of continuum mechanics, including hydrodynamics, and also 
individual aspects of practical applications. Estraordinary topicality of 
problems considered in this book is the only motive for combining 
discussion about these problems into one book in order to attract due 
attention to these problems, especially taking into account the fact that 
giving of many university courses, as a rule, starts from statement of these 
basic concepts and principles. Let us briefly indicate some of them. They 
are: inadequate application of Euler-Leibniz formula while deriving 
equations of dynamics in strains and other continuum equations; 
presentation of Taylor series (inexact differential) with dividing into 
symmetric and antisymmetric parts (which in some textbooks is formulated 
as the first Helmholtz theorem); incorrect interpretation of theorem on 
angular momentum change of particles system in relation to arbitrary 
volume of continuous medium, which resulted in fallacious theoretical 
statement about "symmetry" of strain tensor in continuous medium. False 
simmetry of strain tensor in particular, and included as a component into 
Taylor's series, strain rate symmetric tensor were taken by Stokes as the 
basis for hypothesis about generalized Newton law, fallacy of which is being 
proved in this book. Therefore, formulated in 1845 Stokes equations are not 
correct, due to this reason they may not be used as mathematical models of 
flows of compressible viscous gas, incompressible fluid with transient 
viscosity. Further on, derived in 1904 from these equations in relation to 
flows in stream boundary layer, Prandtl equations contradict to fundamental 
laws of physics: mass conservation law, second Newton law. Proposed, with 
hydrotechnic purposes simple formula of filtration by engineer Darcy, 
obtained incorrect interpretation distorting its primary physical essence, and 
in the result artificial "Darcy law" as a spacial filtering model was created. 
"Darcy law" also conflicts with Newton laws, energy conservation law, 
which excludes possibility of applying this pseudolaw as the adequate 
filtration model. 



Existing from 1882 problem of closing averaged Reynolds equations 
for turbulent flows by engaging additional equations for second, third and 
higher moments, results, as is known, to divergent nonterminating chain of 
Fridman- Keller equations, other semiempirical models groundlessly use the 
whole number of constants and functional associations, in no way 
possessing universality, moreover, physical intensionality etc. 

No doubt, from the point of ' nonsymmetry of strain tensor, necessity of 
revision of Navier-Cauchy-Lame equations in elasticity theory is logical. 

Thanks to development of computational means, equations of 
continuum mechanics found wide range of applications in numerical 
experiments aimed at solving applied problems, and this was the basis for 
origination of its contemporary branch - computational hydromechanics. In 
connection with this, while stating fundamentals of the course of 
computational hydromechanics, it was first of all considered necessary to 
reasonably indicate danger and inadmissibility of applying staggered grid 
in numerical algorithms (for each equation its own grid is used, shifted for a 
half-step from the main grid). Insuperable defect of algorithms with 
staggered grid is the fact that number of unknowns exceeds the number of 
finite-difference equations, and defining of extra unknowns inevitably 
results in appering of unavoidable computational errors that distort desired 
solution of input equations. Along with this, for numerical solution of 
compressible gas equations, difference schemes are applied, and 
calculations under these schemes give negative values of gas density, which 
certainly indicates nonapplicability of such schemes. 

Setting the above-mentioned problems and providing basis for 
possibility of their solving constitute content of this book. 

It is a pleasant duty to express my gratitude to the participants of city 
seminars of mechanics, International scientific conferences and congresses 
on mathematics and mechanics, held at the Kazakh National University 
named by Al-Farabi in 2006-2009 years, to participants of All-Russian 
Conference on Mathematics and Mechanics, dedicated to 60th anniversary 
of Mechanics and Mathematics at Faculty of Tomsk State University on 22- 
24 of September in 2008 year in Tomsk City, for a discussion of the 
presented results. 



Chapter 1. FALLACY OF THE NAVIER-STOKES 

EQUIATIONS. NEW EQUIATIONS OF VISCOUS GAS 

DYNAMICS 

Equations of viscous fluids motion were obtained by French 
scientists: in 1821 - Navier, in 1831 - Poisson and in 1843 Saint- 
Venant. In 1845 the great English physicist Stokes assumed existence 
of linear dependence of strains from components of strain rates 
symmetric tensor and on their basis derived viscous gas dynamics 
equations, later called Stokes equations as opposed to Navier 
equation for incompressible viscous fluid with constant viscosity. 
Due to absence of proper theoretical justification of the given by 
Stokes relationship, L.D. Landau called it Stokes hypothesis. Stokes 
hypothesis leads to paradoxical phenomena that contradicts physics 
of strains, which caused doubts in correctness of deriving of Stokes 
equations on its basis. Thorough surveys revealed a vital error in the 
hypothesis. It is in the fact that Stokes from the Taylor series 

v(r +Sir,t) = v(r,t) + S tSr , transformed to look like 

v(r +Sr,t) = v(r,t) +S<Sr + &£r , given to it by Helmholtz, 

used for determination of viscous strains tensor only one part of 

displacement tensor S = S + & , in particular symmetric strain rates 

tensor S , at this neglecting its antisymmetric part &, which 
characterizes motion component, in particular no less important 
rotary motion medium. It goes without saying, that such neglect is 
paradoxical and alogical. Error of Stokes hypothesis is to a certain 
degree connected with rooted in theoretical physics erroneous 
postulate, according to which continuum strains tensor is always 
symmetrical, (derivation of strain tensor symmetry from moment of 
momentum change theorem is contained in almost all textbooks and 
monographs /l/,/2/,/3/,/4/ etc.) 

Insistent need in theoretical justification in the general case of 
asymmetry of continuum strain tensor, including asymmetry of 
viscous strain tensor, appeared, and, as a result of all this, 
justification of Stokes equations fallacy. Results of researches carried 
out in this Chapter are new equations of flow dynamics with 
variable viscosity with Newton nonsymmetric strain tensor. 



It was easy to see that problems of Stokes hypothesis and 
symmetry of strain tensor are closely connected with paradoxical 
applications of fundamental formulae, in which connection their 
revision here is extremely needed and advisable. 

§1. Paradoxes of Euler, Leibniz formulae 

III describes hydrodynamic derivation of individual time 
derivative for moving volume T of continuum, limited by 
surface <T : 

±lll®ST = -?-ltfct>ST + fl®(v,n)Sa , (1) 

at T ot T a 

referred to in III as Euler formula. 

Hereinafter, according to III, the following denotations are used: 

symbol Sr = Sxi + dyj + dzk - arbitrary infinitesimal segments, 

performed in space at the given moment of time in the given point of 

continuum, symbol dr = dxi + dyj + dzk - infinitesimal 

displacements of fluid elements, derivating within infinitesimal time 

interval dt, v = ui + vj + wk - velocity vector, also 

V = v x i + V 2 j ' + v 3 k . In 111 the following formula is derived 



d 



IIP^ =tff^-ST + fl®(v,n)Sa, (2) 



dt j ; j j ; j dt 

which is used in 13/ as Leibniz formula. The question arises: which of 
these formulae is correct? Answer to this question may be given only 
on the basis of defining a derivative as function increment ratio limit 
to argument increment. For differentiable integrand Ov by x, y, z, 
according to Ostrogradsky - Gauss (divergence) theorem, is brought 
to volume integral: 

^- jjfoSr =jj^ + div(Ov))Sr (3) 

dt J r JJ J r JJ dt 

Further on, formula of velocity of volume dilatation of 
infinitesimal continuum volume St from III will be used very often: 
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dSr „ ,. ^ 

= <?r • gmvv ( 4) 

Rigorous deduction of the given formula is given in III, new proof 
(4) is given here in §2. Paradoxes concerned therewith, that subsist 
two approaches in choice of field of integration T . 

First approach. It was applied by Sedov L.I. in 111, where in view 
of (4) consider, that size T of continuum is locomotive field, that is 
time function r = r{t) , as sum x = \\\Sr of individual sizes 

T 

St = Sr(t) , which depend on time. Let in instant t + At 
locomotive size T occupies status r'=r(t + At). Hence, Sedov L. I. 
in 111 makes up increase ratio 

^ \\JOSt =lim^" { JjJfcC*. y,z,t + At)Sr- jj\0(x, y,z, t)Sr(t) } , 

ut Ar->0 At i 

whence Leibniz formulae comes in limit (2). But in this expression, 
which was made by Sedov, only function O and size t' = r(? + A/) 
taken on instant t + At , wasn't taken into account dependence 
St = Sz(t + At) . If take into account this dependence from time, that 
according to Sedov increase ratio must be constitute in form 
d_ 

dt a»->o At 



\\\®dr =lim-J- -JjJfSK*. y,z,t + At)Sr(t + At) - (5) 

St->0 At T , 

jjj®(x,y,z,t)Sr(t) I 



whence expression comes in limit 

— IjjWr =lim— i\ I f $(*> y>z> t+ ^ St ^ +&)- ( 1 f °(*> y. z. 0^(0 } 

flf r Af->0 Af r , r 

= lim — "if/J^K*' y> z > l + Af )<5r(f + Af ) - jJJ^C*. v, z, f + At)Sr(t + At) + 

At-»0 A? r , r 

+ |JJ0(x, y,z,t + At)Sr(t + At) - j]J®<X y, z, t)Sr(t) ^ 

T T 

= lim — ^\\\®(x,y,z,t + At)ST(t + At)+\\\®(x,y,z,t + At)ST(t + At)- 
A( ^o Af J /_ J T JJ / 
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l\\®(x,y,z,t)ST(t) ~i^(v,H)Sa+l\\j(0ST), (6) 



in view of, that in field t'—t elementary size equals 

_ -, Ar 

8r(t + At) = (Ar,n)Scr, therefore arise speed lim — = v h 

Ar^O At 

t'—t — > (T , ft — ort of outer normal to £7 . 

Well, from (6) follow formula 

J t JJJ&* = JJjf (*>*) + #lXv,n)&r, (7) 

that, apparently, differ from formulas Euler (1) and Leibniz (2). And, 
if begin then and consider, that in locomotive size r'= Tit + At) on 

instant £ + At and coordinates of particles x, y, z will trade places on 
x(t + At), y(t + At), z(t + At) , that increase ratio must look like 

— JTJo<5^- =lim — -i\\[®(x(t + At),y(t + At),z(t + At),t + At)Sr(t + At) - 

Ut At-yO Af , 

- jjjO(x,y,z,t)ST(t) I 

T 

whence after transformation arrives in limit 

— JJJo^r =lim— i\jj®(x(t + At),y(t + At),z(t + At),t + At)Sr(t + At)- 

at A(->o A? , 

— JjT®(jt(> + At), y(t + At),z(t + At),t + At)Sr(t + At) + 

T 

+ jTf®(jt(f + At), y(t + At), z(t + At),t + At)Sr(t + At) - 

T 

-jjjO(x,y,z,t)ST(t) ^ 

T 

= lim — "ijJj^W* + A 0, y(t + At),z(t + At),t + At)Sr{t + At) + 



At 



f jj[0(jc(f + At), y(t + At),z(t + At),t + At)Sr(t + At) - 0(x, y, z, t)Sr(t)] 



3 
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= JJO(v, n)Sa + j]J— (<D St) (8) 



Well, from (8) follow other formula of time derivative 

— Ufa St = JjJ— (OSt) + <$®(v, n)Sa, (9) 

t t a 

that differ from formulas of Euler (1), Leibniz (2) and (7). 

Second approach. It was applied by Loitsyansky L.G. in III. 
Consider, that size T is stated field of continuum, that is not time 
function t'=t = const. Paradox consists in next, that field of 
integration T must be sum of individual sizes St = ST(t) , which 

depend on time in view (4) r = \\\(>t = const . 

Increase ratio in second subsist will so 

— ([jWr =lim — HJJfaW' + A 0, y(t + At), z(t + At), t + At)Sr(t + At) - 

at A(->o At 

-jjj®(x,y,z,t)ST(t) I (10) 

T 

where field of integration not change in time in view r'= r = const . 
From (10) after analogous transformation arrive in limit 

— jjjWr =lim — -i[\\®(x(t + At),y(t + At),z(t + At),t + At)5r(t + At) - 

at T Af->0 At 

t)ST(t) ^ ffflirr^ 



z(t+At),t+At)Sr(t+At)- ®(x,y,z,t)Sr(t)] = JJjA(o^r) 



■d_ 
dt 

T 

In result of stated size r'= r = const arrive formula 



d rrr^ „ cccd 



Si OST = Sh (OST) ' (11) 



dt JJJ JJJ dt 

T T 

that is possible let differentiation under sign of integral. 

Thus, in depend on subsist to calculation total derivative by 
time, arrive 5 different expressions: (1), (2), (7), (9), (11). 
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The fact in III, 111, 131, 141, 151 and the like for opening expression 
— Oc>r - like, apply formulas (3) or (1 1), that in this textbooks, 

T 

constitute foundation of deductive method of derivation of 
differential equations of dynamics, energy balance, continuity 
equation and others. 

Exactly by means of deductive method came in well-known 
differential ratio for moment 

F>p £ ] _F, /1 p ] -tJL[r.*,] = ft 

from which derive mistaken status about symmetry of strain tensor of 
continuum. 

In inductive method derivatives — Oc>r - like aren't 

T 

applied and prove disparity zero: 

[r,A-[r,pF]-X^[r,;r,]*0, 

at ~r { dxj 

thereby prove status about dissymmetry of strain tensor of 
continuum. 

§2.Paradoxes of deformation motion of medium elementary 

volume 

As is generally known, strain rates symmetric tensor 5" 
III, 121, 131, 14 appears in continuum mechanics. Goal of the paragraph 
is proving that nonsymmetric displacement tensor S possesses the 
same properties. In the theory of elementary volume deformative 
motion and with the purpose of establishing viscous continuum 
rheological dynamics laws, Taylor series is widely used 

dv dv dv 
v(r + Sr,t)=v(r,t)-\ Sx-\ Sy -\ di , 

dx dy dz 

at the fixed moment of time t. (Strictly speaking, all theoretical 
mechanics uses Taylor series as principle of virtual displacement of a 
point). Taylor series in matrix-vectiral form looks as follows 

14 



v(f + Sr,t) = v(f,t) + SSr, (1) 

which has a matrix hereinafter referred to as the displacement tensor: 



f 



Qu 


Qvt 


QlA 


C^JC 


<3y 


& 


(3v 


cH? 


c?v 


Qjc 


<3y 


5z 


c5w 


c5w 


e?w 



~\ 



V 



Sjc Sy 



Sz 



<5F r 



f<5x:^ 


r r i 


<5y 


j 


L^ J 


\^k j 



y 



Having denoted elements of this matrix 
S xx =duldx,S yx =duldy,S zx = duldz,S yy = dvldy,S xy =dv/dx, 

S a =dw/dz,S yz =dw/dy,S xz =dw/dx,S zy =dvldz, 

it is possible to represent displacement tensor, by analogy with 

known in hydromechanics strain rates tensor S , as follows: 

^ V V <7 ^ 

'-' xx 1 ^ yx^ zx 

*? <7 *? 

> ~ J xy l ~ J yy KJ zy 

*? V <7 

y '-' xz yz zz j 

Following III, let's introduce three velocities e x ,e , e, of specific 
elongation of fluid elementary vectors <5^(<St,0,0), 
8r 2 (0,Sy,0),Sr 3 (0,0,&), located along axes of rectangular 
coordinate system with beginning in some point M: 

e x = —^(&),e y =—^(Sy),e z =--{&) (2) 

&c dt dy dt di dt 

and three velocities of quadrantal angles shift Y xy > Yyz > ) ' u between 

axes, which were equal to n 1 2 before deformation. Cosines of 
these angles are derived from scalar products 

cos r xy ={Sr l ,5r 2 )l{S>cdy),cosy yz =(dr 2 ,cr 3 )l(&dy), (3) 

cosy zx =($r 3 ,$r l )/(Sxcfc) 
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Time derivatives of these angles are denoted as 

£ xy = -dr xy I dt, s yz = -dy yz I dt, s^ = -dy zx I dt 
The following relation is assumed the basis for defining of kinematic 
meaning of components of displacement tensor S 

±(SF) = sA = Sv, (4) 

at dt 

which results from equalities: 

d d ., ., dr, dr 

— (<y ) = — {r l —r) = = v(r x ) - v(r ) = Sv 

dt dt dt dt 

Meaning r l =r + Sr, formula (1) is represented as 

& = SSF (5) 

Substituting (4) into (5), obtained is equality 

— (SF)=SS : (6) 

dt 

Assuming in it Sr logically equal to <5^(&,0,0), Sr 2 (0, <5y,0), 
Sr 3 (0,0, dz) and projecting to three axes, we'll find 

U&c) = S xx & , Udy) = S yy Sy , ^-(&) = S zz &, 
dt dt dt 

which, in accordance with (2), gives a relation 
S xx =e x ,S yy =e y ,S zz =e z , e x =du/dx,e y =dv/dy,e z =dw/dz (7) 
which results in equalities of diagonal displacement tensor 
component S corresponding to velocities of elongations of line 
elements, located on reference axes and having a beginning in the 
given point of flow. 

Let's calculate derivatives by t from both parts of equalities (3): 

dy xy 1 d(dF v dF 2 ) d ( 1 

-sin/ - = — — + (Sr L ,Sr 2 ) — ( ), 

dt SxSy dt dt SxSy 



dy yi= 1 d(Sf 2 ,Sr 3 ) | ( ^ d ( 1 , 

yz dt SzSy dt dt SzSy 



sin/^ = J-^^^ + (^,^)-(— ) (8) 
dt SxSk dt dt SkSz 
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These equalities in III are applied at the moment of time t = t , 
corresponding to initial strainless state of elementary volume, when 



n 



r xy =r yz = r zx = -,(S? 1 ,s? 2 ) = o,(SF 2 ,sf 3 ) = o,(SF 3 ,sf 1 ) = o 

Will have from (8), because sinuses of right angles are equal to 
111: 



£yz =-Or yz /st = 



1 


d 


<5xSy 


dt 


1 


d 


SySz 


dt 


1 


d 



(<^i,^ 2 ), 

(Sr 2 ,<5r 3 ), 
(Sr 3 ,Sr^ 



(9) 



SzSx dt 

Using (6) and rules of calculation of triple scalar - vector products 
and products of matrix by vector, we have 

at at at 

= (SSr l ) y Sy + Sx(SSr 2 ) x = S xy SxSy + S yx SySx = (S xy + S yx )SxSy, 

— {5r 2 ,5f 3 ) = (S + S )dzdy, — {df 3 ,Sf l ) = (S x , + S a )5z8x, 
dt dt 

Whence, due to (9) derived equalities 



'■Sxy +^'yx' f 'yz 



'■S zy +Syz> £ ZX 



■■s„+s. 



(10) 



In III, 121,131 and others, conforming to formula of the first 
Helmholtz theorem, strain rates tensor is introduced 

r 



Du 

Dx 



V 



if av 

2^ajc 

i_ 

2 



1 ( Du 

-2\dy 

Du i 



Dy ) 

Dw Du ^ 

1 

Dz ) 



Dx 



Dv | If Du 

Dx J 2 ^ Dz 

Dv 1 C Dv 

Dy 2 \^Dz 

1 [ Dw Dv 

2 I Dy Dz . 



Dw^\ > 

DX y 

aw ' 

Dy y 
Dw 



Dz 



(11) 



Along with this antisymmetric tensor is introduced 
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In the overall equality 5=5 + >yis derived. 

Using other representation (first Helmholtz theorem) 

— (5f) = caxSr + SSP , co = —rotv 

dt 2 

in III relations similar to (7) and (10) are obtained: 

S xx =e x ,S yy =e y ,S !Z =e z ,e x = duldx,e y =dvldy,e z =dw/dz, 

— (Sf l ,Sr 2 ) = (S xy +S yx )SxSy, — (Sf 2 ,Sf 3 ) = (S yz +S yz )SzSy, 



d_ 
dt 



(Sf, ,Sf,) = (S 7X + S XT )SzSx, 



(12) 



Due to symmetry of strain rates tensor the following equalities are 
derived 



S xy ~ S yx ~ 2 £ X y ' S yz ~ S zy ~ ^ ^ ' S & ~ S K ~ ^ ^ 



(13) 



Paradoxical is the fact that from both (13) as well as (10) similar 
equalities are derived 

dv du . _dv 8w . _dw du 

XV "^ ' V? ' 7X 

dx dy dz dy dx dz 

despite of striking difference of tensors 5 and 5 from each other! 

Let's consider velocity of relative volume expansion of medium at 
its given point 

9 = --(St), 

St dt 

St - elementary "liquid" volume of medium calculated by triple 
scalar - vector product of elementary coordinate vectors 
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S\ (Sx,0,0),Sr 2 (0, Sy,0),SF 3 (0,0, Sz) : St =( Sr 1 , [Sr 2 , Sr 3 ] ). 
Calculating material time derivative, we obtain 

^ = J_-^(^) = J_A ( ^ l5 [^ 25 ^ 3 ] )= 

<5r J? St dt 

St at St at 

St at 
For elementary parallelepiped St = Sx Sy Sz and by known 

property of unit vectors on coordinate axes i , j ,k , and also due to 

representations 

5r x = Sxi + j +0k , Sr 2 = Oi + Syj + Ok , Sr 3 = Oi + j + Szk 

derived 

[S? 2 ,S? 3 ] = 1 ^Sr 2 ^ ] = J_ t j £] = — 
St Sx Sy Sz Sx Sx 

VSr 3 ,Sr x ] ^ J [Sr l ,Sr^ = k 
St Sy ' St Sz 

Due to which the previous equality will transform into the following 

= — -(St) = — (i,-(Sf l )) + — U, — (Sr 2 )) + — (k,-(SF 3 )) 
St dt Sx dt Sy dt Sz dt 

Using again equality (6), according to which 

— (Sf 1 ) = SSr l , —{Sr 2 )= S Sr 2 , — (Sr 3 ) = S Sr 3 
dt dt dt 

we'll define desired relation 9 in the form of 

St dt ax Sy Sz 

where scalar products are equal to 

— (I,SSr 1 ) = S xx =^,— (],SSr 2 ) = S yy = — , 
Sx ox Sy oy 
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= divv (15) 



— {k,SSr 3 ) = S a = — 

01 dz 

Finally derived 
A 1 d - - - du dv dw ,. „ 

£r<ft ** w zz dx dy dz 

i.e. velocity of volume dilatation of elementary volume of medium at 
its given point is equal to divergence of a velocity vector at this point 
or sum of elongation of per unit length velocities 

6 = e x +e y +e z 

Hence, with the application of displacement tensor S established 
was the formula used above (4) §1. It directly results from (14): 
1 d(Sr ) 

St dt 

or in Euler variables 

d{Sr) d{Sr) d{Sr) d{Sr) „,. _ 

vu h v v w = omivv 

dt dx dy dz 

and expresses rapidity of change of elementary medium volume in 
time with motion assigned by it. 

§3. Paradoxes of integral derivation of equation of dynamics 

continuum. 

In III, 111, 131, 141 apply integration by size T , thus equations of 
dynamics continuum derivate by means of deductive method. 
Exactly applying of deductive method reduced to mistaken 
derivation about symmetry of strain tensor. Dispatch item is law of 
exchange impulse (or number of movements) for system points with 

masses TWj and speeds V ( - : 

d N N - 

where P[= (^V i m i ^)- impulse (number of movements) particle 
with mass ^/ , locomotive under operation of resultant force f i , in 
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which entered all forces, operate on particle with number i, as 
outward as inward forces of interaction particles between 
themselves. By summation these inward forces by third law of 
Newton are cancelled by pairs (see 111). 

In deductive method law (1) is applied to arbitrary size T with 

surface (J , outward normal to which indicated with ?l - Law of 
exchange impulse is written down for arbitrary size T in form: 



d 



jjjvpSr = jjJFpSr + jJ7r n Sa ( 2) 



dt 

t t a 

In III transformation of left part (2) realized by formula (1 1) §1: 

r r 

and in 121 applied Leibniz formula. Apparently, in integral derivation 
of equation of dynamics continuum arise problems with calculations 

— vp^T or derivative from moments — [J\v/?]<5r and the 

r r 

like, on which was indicated in like 5 formulas in §1. Successful 
applying of Leibniz formula or (11) §1 concerned with next, that 
equations of dynamics continuum had to correspond to theorem 
about exchange impulse or 2 law of Newton. 

First, that paid attention, in superficial integral in (2) is taken 
into account only forces, operate on particle, situated on surface (J 

of size T ', that is not taken into account inward efforts 'Mft OCT , 

operate on particles of surface (J ? of individual size St , keeping 
mass / m i = dm, = pSz of aggregate particles. Taking into account 

these forces integral expression must have starting type 

J t W^ = JF^ T + JM^,, S(T (2a) 

T z T cr gt 
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Further, size OT A operate on neighbouring size oT B with force 
f A = cn;r S<7 and, vice versa, size OT B operate with force 

f B = <£j 7T Sa on size OT A , than by 3 law of Newton take place 
equality 



that is so have to be basis (2). From rectilineal set expression (2a), if 
apply Leibniz formula to left part or formula (11) §1, derive equation 
of dynamics continuum in efforts, that is not necessary transition to 
formula (2). 

In theorem about exchange of number of movements (impulse), 
which applied in III to size T like manner (2): 



d_ 
dt 



[[[[r , v/xfr ] = jjj[r , FpSr] + <$[r, n n 8a\ (3) 



equality like (*) isn't take place. In accordance with aforecited 
reasonings rectilineal set theorem about exchange of moment of 
number of movements take type 

j f l\\[r,vpSr]= \\\[r,FpSt] + jJJ[F,fr>] (3 a) 

i r r a Sl 

From (3) derives equality, used in III and other, 

[r,p^-]-[r,pF]-Y^[r,7f.] = 

dt ~[ oxj 

with one's following from here oppositions about symmetry of strain 
tensor of effort continuum. From rectilineal set expression (3a) prove 
ratio 

[r,pA-[r,pF]-[r,t^-] = 0, 

dt ~j ox j 

from which follows equation of dynamics in efforts 
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dv_ = p | y dTTj 

dt j~{ dxj 

therefore, don't question about symmetry of strain tensor of efforts, 
that took place in wrong formula (3). 

Other paradox consists in next, that for derivation of fundamental 
mechanics equations of continuum, not any necessary in integral 
formulas like (2), (3). Show this. 

Let mass of individual size St equals dm = pSz , F - density 
of body forces, 71 n - effort, p — density of medium. Theorem about 
exchange of impulse necessary formulate directly for size oT , take 
into account, that main surface force, operating on surface O ' g T of 

size St equals jj^n^ , main body force, operate on size St , 
equals FStn : 



d_ 
dt 



(v 3n) = FSn + <x\7r n S<j (4) 



<*>>> 



Further by theorem by Ostrogradsky - Gauss make transition to 
volumetrical integral and smallness of individual size St prove 

III ( — - + — - + — ^-)Sr = ( — - + — - + — ^) \\\St = 

St 8x 8y 8z 8x 8y 8z 8r 

ox oy oz 
therefore, equality takes place 

%Zte-A 8 A + SJL )S r (5) 

a^ dx fy dz 

Substitute (5) and Sm = p5l in (4), find 
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— (vpSr) = FpSr + (— ^ + — x + — ^)£r (6) 

a? ox qy dz 

Left part (6) is transformed on base of formula (4) §1 and equations 

of continuity dpi dt + pdivv = : 

dv _ _ _ <iyC> 

/? <f>r + yOV <97mw + v e>r — ^— = 

efr * dt 

= Fpdr + ( — - + — - + — L )£>r, 
dx dy dz 

where from after cancellation OT prove classic equation of dynamic 
continuum in efforts 

dv - dn x dn dn 7 
p — = pF + — - + — - + — z - 
dt dx dy dz 

§4. Inductive approach 

Inductive approach (from the particular to the general) is free from 
the given above disadvantages because notion of continuum particle 
is used. "Liquid" individual volume St is formed by these particles 

with masses Ifl^ and velocities V; . (If inductive approach was used 

from the very beginning of development of continuum mechanics, 
there would be no problem of strain tensor symmetry! In this 
connection this physical approach is detailed here.) 

It is assumed that due to medium continuity, St containes the 

sum of particles V"m ; = dm , weight- average velocity v of volume 

St is defined as relations 111 

v =^ i m i v i l^m; =^ i m i v i I dm , in the same manner weight- 

i i i 

average force affecting St is determined as 

F = 'Y^ m i F i FYjtli =^,tn i F i I dm, whence vdm = yw.y. , 
i i i i 

FSn = ^m,F. , p - density p = dm I St , dm = pdt = ^J n i . 

i i 
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Theorem on system of mass points impulse change is applied to 
elementary volume of continuum St, containing aggregation of 

particles /m j , and not to final volume r=f|f<5r, strictly 

' T 

speaking, this is the essence of the inductive approach: 

~rX m i^i = Y,Fi m i + j^T^nk^k , (1) 

at i i ad[ k 

where M n OCT — ^^nk^k - resultant surface force, affecting 

k 

elementary surface &7 = V<7 t , at this OCT E Cg T , CT^ — site of 

k 

surface Scr , which is occupied by particle m^ , that is under 

action of strain K n i , <J ^ = <fc| 8<J . 

Formula (1) in equivalent representation takes the above-used form 
— (v dm) = FSm + ^7r n Scr , (2) 

which results in, as was shown, continuum dynamics equation under 
strains. 

But it is also possible to make it in a more simple way, without 
applying theorem on integral mean value, and using passage to the 

limit, St = SeSycZ — >0, if, elementary volume 

St = SxSydz, ^ s taken as parallelepiped with faces, which are 

parallel to coordinate planes, where differences are equal 

Sx = x 2 — x x , Sy = y 2 - y\ , & = z 2 - z r 

Impulse change law for parallelepiped with sum of particles 
/ t rn t =om = poT shall be represented with account to bulk and 

surface forces affecting St = Sxdydz . In the result, in application 
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to parallelepiped, theorem on impulse change will take the form 
which is similar to (2): 

— (vSm) — FSn + (7r xl + 7r_ xl )Sy&+ (3) 

dt 

+ (7T y2 + 7T_ yl )3cdl + (7T z2 + 7T_ zl )3cSy , 

where, just as previously, 

v&n = Ym^Fan = Yjn i F i ,n n S<j = JX,<7, 

/ i k 

Comparison of (3) with (2) shows that surface integral for 
parallelepiped is calculated by its faces 

<r & = (#&), U (Sy&) 2 U (SxSz\ U (Sx&) 2 U (Sxfy\ U (<&&) 2 

In the result, the following value of integral is obtained 

§Z **°k = fe + ii_ xl )fy& + {n yl + n_ yl )8xdz + 

<?8t k 

+ (.7T z2 +7T_ zl )Sxdy 

Using formula 

ddr ,. - e 

= divv or 

dt 

and equality of strains 

X- X i = -n x \ , %_ yX = -n yl , n_ zl =-x zl , (4) 

having divided by St = §X§ySz > we derive formula 

p v v(— — + pdivv) = pF + (7T x2 — 7Z xl )l <5>c + 

dt dt 

+ (7T y2 -7T yl )l dy + (7T z2 -7T zl )l & (5) 

It is obvious that passage to the limit in (5) Sx — » 0, dy —> 0, 
fa _». o, gives continuum dynamics equation under strains 



26 



dv ^.dp .. _ - dn x dx dfc 7 
p — + v(-^- + pdivv) = pF + — - + — - + — z - (6) 
dt dt ' ' dx 8y dz 

When continuity equation is equal to zero (no sources and drains) 

1- pdiw = 

dt 

this equation (6) exactly matches with equation (10) from §3. 

Theorem of Ostrogradsky-Gauss and theorem on integral 
average value are not used here. 

Energy balance equation is derived in a similar manner. For 
elementary parallelepiped §j = SxSyfiz of continuum, energy 

conservation law is formulated as follows: 

d ,2 

— [(E + |v| /2)Sm] = (FSm,v)+[(7r x2 ,v {Xi} ) + (7r_ xl ,v (Xi) )]SySz + 

+ [0^2 » Vf 2 j + {7T_ yl , V €i j]SzSx + 
+ [(7T z2 ,V^ 2 j + (7T_ zl ,V^]ScSy- (7) 

- ~iq (X2) -q (x0 )Sy& + (q (y2) -q (yi) )&Sx + 
+ (q {Z )2- ( l(z)i)S )l dy ~}Q3n 

where V^. ~", Tfl = L,Z - values of velocity vector on areas dydz 
in sections Xj , x 2 etc., q = q^ x) i +q, y \j +q, zy k - heat flow vector, 

(is+lvl I2)dm - total energy of volume St, (FSn,v) - 
bulk force power, 

[fe>^ 2) ) + (^_ xl ,v Ui) )]<5y& + 
+ [(x y2 ,v (y2) ) + (7r_ yl ,v (yi) )]&S>c + 

H(x z2 ,v (Z2) ) + (7r_ zl ,v (z0 )]$>cSy 

sum of powers of surface forces affecting pairs of faces 
Sydl, Sl&C, SxSy of parallelipiped, 
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- iq (X2) -q (Xi) )dy& + (q (y2) -q iyi) )&Sx + (q (Z2) -q (Zi) )Sxdy 

- heat flows through these faces, Q5m — energy source or power 

consumer in volume §T ■ 
Having divided both parts (7) by §t = SxSySz an d subtending 

parallelepiped to point fa — > 0, Sy — > 0, 5l — > 0, we derive 
energy balance equation in continuous medium: 

p — (E + \v\ 2 1 2) + (E + Ivl 2 / 2)(— + piz'vv) = (8) 

~ _ 5^^ 5^^ 5^^ 
= p(F,v) +— (tt x ,v) +— (n y ,v) +— (7r z ,v)-divq + pQ, 

ox oy oz 

from which, in accordance with basic Fourier conduction law 

q = —AgradT and dE = c v dT , for Newton strain nonsymmetric 

tensor /8/ obtained is heat influx equation 

dT 3 3 <9v 1 

pc v — = div(lgradT ) + /uY V( — L ) 2 - pdivv -{- /u- /u'){divv) 2 
dt -tf-p 1 5x 7 - 3 

(Heat energy crossing area OCT in unit time is equal to {q,n)d<7). 
Given above inductive approach possesses direct link with 
fundamental laws of physics. It is necessary to note that continuity 
equation (6) was derived in III from law of conservation of matter 
with the help of inductive approach. According to law of 

conservation of matter, mass g m = pg T = ^ m , of individua i volume 

i 

St is a constant Sm =const, due to which 

d3n =Q d(pST 1 = 0p diST 1 + ST dp =0 ((?) 

dt dt dt dt 

dSr _ ,. _ 

Due to = or • divv (9) results in 

dt 

pSrdivv + St = O , 

dt 

having reduced St, we come to continuity equation 
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— — + pdivv = 
dt 

The purpose of the above-given two paragraphs was: to show 
disadvantages of using deductive method and accuracy of inductive 
approach. 

The point is that in deriving false symmetry of continuum strain 
tensor, deductive method is traditionally used. 

§5. Paradoxes of the first Helmholtz theorem 

It is known III that Helmholtz, excluding pressure p from the 
Navier-Stokes incompressible fluid equations III, obtained equations 
for velocity rotation 

„.^,-t /3vv dv.~ ,du <9w\-r ,dv du.7- 

rotv = ( )i +( )j+{ )k 

dy dz dz dx dx dy 

Bearing in mind that dimensionality rotv coinsides with 
dimensionality of angular velocity cb , Helmholtz tried to connect 
them with the help of velocity formula for solid body points /I/: 

v = v +[<»,(r -/-,)], (l) 

where \? , r^ velocity and pole radius vector in relation to which 
body instant rotation occurs at the moment. Effort to express 
components of angular velocity CO = CO x i + CO y j +CO z k via 

components of linear velocity v =ui + vj + wk led Helmholtz to 
solving in relation to CO x , CD y , C0 Z in linear algebraic system 

u=u Q +co y (z-z )-a z (y-y ), 

V = V +0) Z (X-X )-(0 X (Z-Z ), (2) 

w = w Q +o) x (y-y )-o) y (x-x Q ) , 

determinant of which is equal to zero, therefore the given system 
with assigned 

du =U-Uq,Sv=V-Vq,Sw = W-Wq, 

Sc = x-x ,Sy = y-y ,Sz = z-z 
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has infinite number of solution sets, which first of all results from 
equality 

I ^ ~~ ^0 | = | 1^' ^~ *b^_\ = |<5)|^ — ?^| sin a; 

system consistency results from orthogonality condition v — Vq to 

r — Fq: (v,r — Fq)=0. Indeed, if co x is selected as reference 

arbitrary variable, solution set for system (2) is represented in the 

following form 

dv Sz by dw 

a> 7 = — + co x — ,co v =co x (3) 

ox ox ox ox 

Helmholtz {Loytzyansky in/1/?) acted absolutely differently. He 

differentiates (2) by x,y,z, assuming that v Q , f\j , (0 , are constant 

and as the result he obtained the following formulae for components 
of solid body angular velocity: 

dv dw du 

dz ox dy 

dw du dv 

° ) x=—^y=—^ Z =—^ (4) 

dy oz ox 

Paradoxical, but use of method (4) Helmholtz (Loytsyansky?) to the 
solution of linear equations system 

N 

i=\ 

gives values of desired quantities in the form 

db 
x. = — '-, j = l,...,N,i = l,...,N , i.e. obtained are N*N of values of 
da tj 

desired unknowns xj despite of prescribed N. It is obvious that they 
may form any combinations of the following type 

m mj I "' 

where OC m - arbitrary numbers. In this case, in systems (4), this fact 
may be used as follows. 
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Multiplying upper line by n - lower line by m, adding and 

dividing them by m+n ^ , will find values 

m dw n dv 
co x = , 

m + n dy m + n dz 

m du n dw 

( °y = a ~^~' (5) 

m + n dz m + n ox 



m dv n du 



co z = 



m + n dx m + n dy 

Helmholtz (Loytsyansky?) from this infinite diversity of angular 

velocity components used only one aggregation resulted from (5) 

with m=n = \\ 

dw dv du dw 
(rotv) x =— — — = 2co x ,{rotv) =- — = 2co y , 

dy dz dz dx 

. _ dv du . 

(wtv) z = — -— = 2a} z , (6) 

ox dy 

where known expression of velocity rotation was used. 

It is obvious that with no values of m, n , expressions (5) 

co x ,G) y , co z do not coinside with solutions (3) for system (2). Even 

if to take as value of unrestricted variable 

1 .dw 5v N 
co x =-( ), 

2 dy dz 

substituting it into found solutions (3) for system (2) 

dv dz dv 1 .dw dv^dz 1 ,dv du^ 

co z = — + o) x — = — + -( )— *-( ), 

dx dx dx 2 dy dz dx 2 dx dy 

dy dw 1 ,dw dv. dy dw I ,du dw. 

co =co x — = -( )— *-( ) 

dx dx 2 dy dz dx dx 2 dz dx 

we'll receive evidence that Helmholtz formulae (6) ) are not 
solutions of system (2), QED. Paradoxical fact, using incorrect 
solution (6) Helmholtz instead of formula (1) refers to expression 

v = v +[— rotv,{r -r )], (7) 
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which has nothing in common with formula of solid body velocity 
(1), since (6) is not the solution for system (2). 

Due to the fact that Helmholtz formulae (6) are not solutions for 
system (2), this expression is far from being equivalent to formula 
(1) of solid body particle velocity! (This formula is erroneous, it is 
not suitable for calculation of solid body particles velocity.) 
However, this expression became later a prototype for Helmholtz 
first theorem, on which we will focus in more details. 

Helmholtz first theorem is derived from approximation formula, 
presenting members with first derivatives of Taylor series in the 
following form 

v(r -\-Sr,t} — v(r,t) -\- S Sr (8) 

Helmholtz, by analogy with incorrect expression (7) transformed (8) 
to equivalent form (refer III) 

v(r + Sr,t) = v(r,t) + — [rotv,Sr] + SSr, (9) 

having introduced strain rates tensor S . Formula (9) is the content of 
the Helmholtz first theorem. Comparing (9) with incorrect expression 
(7), Helmholtz declares S Sr to be a deformation displacement. If 

to bear in mind that expression v = v +[— rotv, ( - r , 

separately does not have any physical sense at all, as opposed to 
solid body point velocity V = v + [d>, (r — r )], then formulated 

notion of deformation displacement is very doubtful from the 
physical point of view. 

More contensive is displacement SSr in (8) and because of 
results (7), (9), (10) §2, this product may be declared a deformation 
displacement. Let's show that Taylor series (8) may be assigned with 
infinite number of forms containing rotv , thus, it would be proved 
that there is the infinite number of "deformation displacements" of 

Helmholtz type S Sr . With this purpose we'll represent Taylor 
series (8) in projections on coordinate axis: 

/ = c c- \ / . du „ du „ du „ 
u(x + ax,y + oy,z + oz,t) = u(x,y,z,t)-\ ox + — oy + — oz, 

dx 8y ' dz 
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dv dv dv 
v(x + &c,y + Sy,z + &,t) = v(x,y,z,t)-\ Sx-\ Sy-\ dz, (10) 

dx dy dz 

dw dw dw 
w(x + dx,y + 5y,z + dz,t) = w{x,y,z,t) + — 8x-\ dy + — dz 

dx dy dz 

Equivalent transformation of Taylor series (10) has the following 
form: 

, - c - N . N du b-\ .dv du. _ 
u(x+ax,y + (fy,z + <x,t) = u(x,y,z,t) + — ck ( )<5y + 

dx b dx dy 

b-l,du dw.„ .b-ldv ldu.„ b-\dw ldu.„ 

+ ^(—-—)& + (^ — + T —)dy + (^ — + T —)az, 



b dz dx 



b dx b dy 



b dx b dz 



i e = = ^ i ^ b-l.dv du.„ dv 
v(x+ac,y+dy,z+cz,t) = v(x,y,z,t)+ ( )ox+ — oy- 

b dx dy dy 

b — \,dw dv. _ ,b—\du 1 5v N _ ,b — \dw 1 5\> _ 

-( )& + ( + )& + ( + )&, 



£> dy dz 



b dy b dx 



b dy b dz 



t a c c- \ /■ x b-l.du 6w. Q 
w(x + &, v + by, z + oz, t) = w(x, y, z, t) ( )ox + 

b dz dx 

b-l.dw dv. . dw . ,b-\du 1 dw, , ,&-15v 1 dw, . 
+ ( )# + — & + ( + )bx + { + )by, 

b dy dz dz b dz b dx b dz b dy 

By analogy with presentation of Helmholtz (9) Taylor series (1 1) 
may be represented in vector form: 

b-\ 

v(r+Sr,t) = v(r,t) + [rotv,Sr] + S h Sr , (12) 



where in Sj^Sr (product of line S^ by column dr ) there is a 



matrix 

du ' 

dx 

r 



b — ldv I du 

+ 

b dx b dy 



b—ldu 1 dv 

+ 

b dy b dx 

b—ldu 1 dw 

+ 

b dz b dx 



dv 

~dy~ 



b — \ dw 1 du 

+ 

b dx b dz 

b — \ dw 1 dv 

+ 

b dy b dz 



f b-\ dv 1 b\v^ dw 
+ - 



b dz b dy 



dz 



fj~\ 



or 



Sx 
by 

ySZj 



I 



Jj 



(13) 
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b-const, b ^ 0,1 b l< oo. When b-2 from (12) the first theorem of 

Helmholtz (9) is derived, since S 2 = S , i.e. tensor 5 2 is equal to 

strain rates tensor, with b=l universal formula (12) transfers into 

Taylor series in initial representation of (8), since S l = S , i.e. 5", is 

equal to motion matrix. Everywhere above, velocity rotation has 
components 

dw dv . du dw . dv du 

(rotv) x = — - — ,(rotv)= — - — ,(rotv) 7 = — -— (14) 

dy dz dz dx ox dy 

Therefore, expansion into Taylor series (10) may be given another 
equivalent form: 

du 

u(x + St, y + Sy, z + &,t) = u(x, y, z, t) H &c — 

dx 

.dv du.„ .du dw.„ dv _ dw _ 

- ( )dy + ( )8z + — Sy + — Sz, 

dx dy dz dx dx dx 

v(x + Sx, y + Sy, z + Sz, t) = v(x, y, z,t) + ( )Sx + 

dx dy 

dv .dw dv^ 3m c dw 
H Sy - ( )Sz H Sx H f5e, 

^y dy dz dy dy 

w(x + Sx,y+Sy,z + Sz,t) = w(x, y, z,t) - ( )Sx + 

dz dx 

.dw dv. _ 5w „ 9m _ dv _ 

+ ( )<^H &H <5ch Sy (15) 

<?y & & dz dz 

By analogy with (12) projections (15) are reduced into expression 
v(r +Sr,t) — v(r,t) + [rotv,Sr] + SrS (16) 

where 5 - is nonsymmetric displacement matrix from §1, 3f S 
contain products of lines by columns. Helmholtz called the following 
formula a vorticity 

a> = ll2rotv , (17) 

basing on false presentation (7), which was discussed above. 
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As is known III, IAI, velocity rotation in hydrodynamics is 
connected, first of all, with velocity circulation: 

j(v,dr) 

(rotv,n)= ]im ( 18 ) 

and components rotV in the given in (14) form are derived in 
particular from this determination (18), therefore formula 

Q = rotV must be used for vorticity. If to proceed from 
determination of Helmholtz, deformation displacements are 

SfoSr , contained in universal formula (12). 

§6. Paradoxes of Stokes hypothesis 
1°. Stokes hypothesis was based on wrong Helmholtz formula 

v — v +[l/2rotv,(r—r )] 

§5 describes this fallacious Helmholtz formula, on the basis of 
which extension into Taylor series is formulated in the form of 
Helmholtz first theorem /I/: 

vif + Sr,t) = v(r,t)+l/2[rotv,$r] + SSr 

Further on, Stokes, having admitted in this Taylor series the 

summand S Sr as deformation displacement occurring under 

action of surface forces, theoretically assumed that tangential stress 

are proportional to doubled values of components of strain rates 

tensor S : 

thanks to which, tangential stress happened to be symmetrical due to 
symmetry of matrix S . Stokes neglected the equivalent member of 
Taylor series (this is the major error of his hypothesis): 
\/2\rotv,Sr ^Je &Sr . 

His hypothesis Stokes (another name - generalized Newton law 
III), formulated for flow of viscous fluids and gases in tensor form 

7r c =-(p + 2/3judivv)E + 2juS, (l) 

In index notations components of strains tensor n c are represented 
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in a shorter form 

dv- 
n ii(c) =~(P + 2/ ?>[Jdivv ) + 2ju— L ,i= 1,2,3, 

dx i 

dv. dv , 

J * 

and symmetry of strain tensor in the form 

^ijic) = n ji(c) ' i * j, i, j = 1,2,3 (3) 

Noticed is the circumstance that diagonal S tensor elements enter 

normal strains with doubled value, which is connected with slanting 

dv 
of tangential stress (2) into Newton friction law 7t m = JU — — . 

dn 

For instance, Newton friction law with v s =u,n = y,s = x gives 

du 
formula n yx i H \ = jU — , and according to Stokes hypothesis (2) it 

dy 

comes out that n , , = u( 1 ) . Apparent mismatch is 

y{) dy dx 

dv 
connected with inclusion of paradoxical excess member jU — . (By 

dx 

the way, inconsistence of Stokes hypothesis was indicated by 

L.D. Landau.) 

„ dv t 

Here arises question about adequacy of member 2/u in normal 

ab- 
stain in the sense that selection of coefficient "2" is justified only by 
Stokes hypothesis (1), according to which the following artificially 
created equality must be realized 

K Mc ,=2S fi J*j (4) 

In the preceeding paragraph § 5 for Taylor series 

v(r+Sr,t) =v(r,t}-\-SSr 

was given an infinite number of equivalent formulations, containing 
velocity rotation, 
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b-\ 

v(r + Sr, t) = v(r , t) H [rotv, Sr] + S b Sr (5) 

b 
If to follow logic of Stokes hypothesis, strains selected on the basis 
of considerations of their proportionality to components of tensor 

Sjy , will be defined in the following form 

7r xx = -(p + b/ 3/Lidivv) +bju — ,n =-(p + b/ 3/jdivv) +bju — , 

dx ' dy 

7t zz =-(p + b/ 3/jdivv) + b/j — , 

dz 

b — \dv 1 S«. b — 1 dv 1 dw. 

7T = jU{ 1 ),7T = jU{ 1 ), 

b dx b dy b dz b dy 



,b — ldu 1 <9v x b — ldu 1 5w, 

v 1 )> 7t xz = M ! J 

b dy b dx b dz b dx 



.b — ldw 1 du. ,b—ldw 1 dv. 

7V ZX = jU{ 1 ), TV = f-iy 1 ) 

b dx b dz b dy b dz 

For b ^ 2 tangential stresses are nonsymmetric, symmetry occurs 

onle with b=2, when S 2 = S. It is also obvious that when b=l 

matrix (6) transforms into displacement matrix S 1 = S . 

There is the infinite number of viscous flows (only some of them 
are given here), where symmetry of Stokes strains (3) contradicts to 

dv s 
fundamental Newton friction law Tt m = jU — - , where n crosswise 

dn 

to ? direciton, v s = (v, s) - projection of velocity on ort direction 

?, ^ ns - tangential stress. 

2°. Paradox of Stokes hypothesis 5*2 = S in Poiseuille and 

Couette flows 

In laminar flow between parallel tight and rigid channel walls 
(Poiseuille flows), velocities are defined as solution of equation 
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d P d2u t+u\ n 

— = ju-—, u(±b) = 0, 

dx ay 
in the form of dependences 

u = -(b -y ),{— = const), v =0, 

2ju dx dx 

{here V\ =u,v 2 =v,v 3 = w,x { =x,x 2 =y,x 3 =Z). 

In this flow, longitudinal tangential stress towards x by formula (2) is 

d« dv dp 

equal to 7t yx i c \ = ju( 1 ) = — y. Due to the fact that there is 

dy dx dx 

no flow along crosswise direction V = , which means that 

dv 

derivative is equal to zero = , according to Newton friction 

dx 

law crosswise tamgential stress towards y will also be equal to 

dv 

dx 

of cross flow towards y, since v = along the whole channel. 
According to Stokes hypothesis about symmetry of tangential stress 

^yx(c) = ^xy(c) . there is a paradox, which is in the fact that 
crosswise tangential stress is not equal to zero, because 

dp 

7T xy(c~) = 7T yx(c) = — J ^ » which conflicts with fundamental 
dx 

Newton friction law, according to which, as was proved above, 

crosswise tamgential stress is equal to zero TZxyiu) = 0- 

Thus, in Poiseuille flow crosswise and longitudinal tangential stress 
are not equal to each other, i.e. their symmetry in relation to left 

principal diagonal of strain tensor 71 c in the form (2) does not take 

place. This contradiction is quite similar to stated below paradox 
with tangential stress by Stokes (2) in Hagen-Poiseuille flow in pipe. 
For Couette flow, nonsymmetry of tangential stress is defined 
similarly. 
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zero: ^xy{H) ~ J^~ — , which can be fully explained by absence 



3°. Paradox of Stokes hypothesis 5*2 = S in arbitrary flows 

For infinite number of flows symmetric tangential stress in Stokes 
hypothesis (2) has zero values, i.e. in all points of the flow are equal 
to zero 

dv dv . 
7r ji(c) = Mi-^- + -=— ) = 0,i^ j, i = 1,2,3, j = 1,2,3 

J I 

Let's restrict ourselves with giving a short list of flows with 
velocity components, in which this fact take place 

1) u=F(sin kiX cos kiy ), v=F(-coskiX sink^y ), 

2) u= U(sinli2X cosli2y -cosli2X sinlt2y ) , 
v= U(sin &2-* cos ^2y -cos /^-^ sin A: 2 y ), 

3) u=W(-cos k^xsin kj y ),v= W(sin k^x cos & 3 y ), 

4) u=Q(sin k 4 x sin k^y ),v=Q(cos k A x cos k^ y ), 
5 )u=T(sin k$X sin k$y ),v=T(cosk$X cosk^y ), 
6) u x =M(sin kftX sin k^ y +cos k^ y cos k^x ), 
v=M(sin kfrX sin k^ y +cos k^ y cos k^x ), 

7)u=S(e ), v=-S(e ), 

for three-dimensional flows: 8)u=D((e % y -e 8Z )e gx J, 

v=D((e k * z .e k * x ) e k * y ), w =D((e k * x .e k * y )e k * z ), 
where coefficients kf = const, i = 1,2,3,4,5,6,7,8, are selected 
arbitrarily from infinite interval -co<ki < + 00 • Standing here 
differentiable functions F,U,W,Q,T,M,S,D are also arbitrary in 
choice. It is obvious that new any linear combinations of type 
u=F+U, v=F+U etc. may be formed from the indicated list. Velocity 
fields 1-7 correspond to plane flows and satisfy two-dimensional 

continuity equation 1 = , velocity fields 8) satisfy three- 

8x dy 
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minesional continuity equation 1 1 = . 

dx dy dz 

For all these flows, Stokes tangential stress is identically equal to 
zero in all flow points: 

dv. dv i 
^jnc) = Vi-^- + -r- ) = 0,i^j,i = 1,2,3, j = 1,2,3 , 

J ' 

For compressible gas, number of flows where zero strains of 
Stokes take place 7t^ c = 7T - c = 0,1 ^ J , 

infinitely increases due to presence of variable density yO in 
continuity equation. Thus, in flows with components of velocity of 
1,2,3,4,5,6,7,8 type, symmetrical strains (2) reduce to zero, 

^ij(c) ~ ^ , ' ' ^ h an d it comes so, that viscous fluid motion 
occurs without friction, which contradicts to fundamental Newton 
law 71 ns = /Jdv s I dn . It is easy to calculate by this formula that 
tangential stress in the indicated flows is not equal to zero 

dv; _. 

4°. Paradox of Stokes hypothesis $2 = S in Hagen-Poiseuille 

flow in a round pipe 

Viscous fluid laminar flow in a round pipe (ref./l/J in cylindrical 
coordinates has the velocities 

4ju dz dz 

where a —pipe radius, z - axial, r- radial coordinates. According 
to Stokes hypothesis (2) symmetrical tangential stress are equal to 
each other and are calculated by formula III 

,dV r SV 7X 1 dp 

'™-*-<»-'*-at + -* ) -2* r <7) 
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V z =-- ^(a 2 -r 2 ),V r =0,V^=0,^ = const, 



Let's consider flow in positive direction towards axis Z, it 

appears when pressure drops —±— = const < O • 

dz 

At this, by formula (7) longitudinal tangential stress are negative 

^rz(c) "^ ^ > by Stokes hypothesis (although flow in crosswise 

direction r is absent) due to symmetry (3) exist and are negative 
crosswise tangential stress 

n zr(c) =7T rz{c~)i n zr{c)<® 
Their directions are shown in Fig.l. 

<= -z +z^> 



pipe symmetry axis 



= 



n zr{n) 
dz 

<0 



<r~7T 



rz{n) 



pipe wall 



+ r 







<- 



z. 
^rzic} 



o 



t;r zr(c) <0 



+a 



w 



Fig.l 



According to fundamental Newton friction law, crosswise tangential 
stress are equal to zero 7Z zr(n} = JLldv r I dz = , because 

V r = , longitudinal tangential stress are negative 

7T rz{u) =/jdv z /dr = l/2(dp/dz)r<0, 
and both - tangential stress by Stokes hypothesis as well as tangential 
stress by Newton formula 7Z rz{c) <\j,7Z n(H) <\J have equal 
directions up the stream. With positive gradient dpi dz>0 fluid in 
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pipe flows in negative direction of z. axis (right to left), tangential 
stress change signs due to (7) 

Kfzlc-i >0 '^(O >0 



According to Newton friction law 



dv„ 



oz 
and flow will be directed " — z ", which is displayed on Fig.2. 

180° turn of Fig.2 results in Fig.3, which displays that flow in pipe 
is directed in the same direction as on Fig.l. 
<=-Z + Z => 




pipe symmetry axis 
y z 

— » 



+a 



= o 

^rzin) > 
pipe wall 



+ r U Fig.2 

Contradiction is in location of nonzero tangential stress by Stokes 

hypothesis ^ zr(c} , which in Fig. 3 is directed towards pipe 

wall, and in Fig.l is directed to pipe axis with equal flow direction, 

whereas strain ^"rz(c) in both cases is directed up the stream! 

On the contrary, according with Newton law tangential stress is equal 
to zero: 

n zr{n) — f^~^~ ~ u > 



dz 
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since V ,. — O because of absence of radial flow. 



+z 



-z 






dz 



\^. 



n zr{H) ~ ° 
dp 



O 



+ r 



II 



p/pe symmetry axis 



">->-> 



V. 



7T zr{c) > 



^ <- ^rz(c) > ° 



+ a pipe wall 

Fig.3 



It is obvious that for nonsymmetric tensor (6) contradiction 
confirmed by Fig. 1-3, does not arise, while for symmetric Stokes 
strain tensor (2) there is an indicated in Fig. 1,2,3 paradox. The same 
paradox with directions of symmetric Stokes tangential stress, 
obviously obtained in Poiseuille and Couette flows. 

§7. Nonsymmetric Newton strains tensor. Paradoxes of 

determination of viscous normal stress are related to hypotheses 

on pressure and to Pascal law for ideal liquids 

Law of linear dependence of strains on strain rates proposed by 
Stokes in 1845 as inadequate generalization of Newton friction law, 
can be related to Helmholtz formula or Taylor series, which is the 
same: 

vif + <$r,t) = v(r,t) + l/2[rotv,cT] + S(r , (1) 

equivalent representation (1): <9P = &S? + SSr . 
This formula needs to be presented detailed form 
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^ 1 3v. dv j 1 dv. dv, 

~T X 2 dx^ &. 2 ckj. Qx;. 

Referring (2) to infinitesimal period St , we obtain convectional 
acceleration component 

Sv. 3 i dv- dv t 1 dv. dv f Sx f 
ctf T^f 2 dxj dx t 2 Sx, &; dt 
On the basis of (2) Stokes suggested a hypothesis that strains are 

proportional to doubled deformation displacement 2.St5r . The 
hypothesis was formulated in the form of generalized Newton law 
with symmetrical strains tensor 

2 dv. dv, 

n m =-(p + -judiw)S ij +ju(— L + -^-),i,j =1,2,3, (3) 

J ' 

where S tj - Kronecker symbol. As was marked, appearance of 

multiplier "2" in 2 S Sr is connected with fitting tangential 

dv- dv ■ 
stress by Stokes (3) 7T- i(c) =//( — L H -),i ^ j to Newton friction 

dxj dx t 

dv. 
law (6) §6: n H(H) = ju — -, i =£ j . In this connection there 

dXj 

arises a question: if strains (3) are caused by deformation 

J^l dv. dv. 
displacements > — ( — '- -\ -)Sx , then the same strains cause 

j-}2 dxj dx t ] 

1 _ 'l Qv. dv. 
rotation displacements in {2) — [rotv , df\ => / — ( — '- -)Sx- ? 

2 "p*2 cbc^ 3x. 

Especially so because, rotational motion includes the same gradients 

dv i dv j . . 

, , that are also in 7t;: lr \, i ^ / , moreover 

ok ■ ox,- 

> — ( — ' -) — - is the constituent of convectional 

~~^2 dxj dx t St 

44 



acceleration ! 

Paradoxical neglect of Stokes of these forces was motivated by 
fitting of the hypothesis (3) to erroneous statement of symmetry of 
strains tensor. (Nonsymmetry of continuum strains tensor will be 
rigorously proven in the following paragraphs.) 

Answer to the question arisen above is the insistent need of 
considering strains 

ft* = -u( — )> 

11 r d Xj dx/ 

Corresponding to rotational displacements of medium, equally with 
strains 

« 1 ,dv t dv: 

ft a =~M( — + — -). 

2 dXj dx { 

arising in flow due to deformation displacements, i.e. it is impossible 
to neglect such important component of general motion as 

— \rotv,Sr] = &SP . In the result of this, total force is determined 

2 

in the following form 

ftjun) = ~(P + -/*livv)S ji + n\ + TV**, 

that after substitution 7Z-. , tt,-. leads to nonsymmetrical strains tensor 

1 dv. 

ftjHu) =-(p + -judivv)S ij +ju— -,i,j = l,2,3 (4) 

3 ax j 

This formula automatically corresponds to both - Taylor series 

4^dv j , dv 

ov = / ox. as well as Newton friction law tt = ju — - and 

~7 ox i on 

gives nonsymmetric values of tangential stress (6) §6: 

dv,. 

8X; 



71 'ji(s) = M — L ' i ^ J ' ■ m Cartesian coordinates, it will look as follows: 



./ 



du dv dw 

ftyx(H)=M — > ftxy(H)=M—, ftxz(H)=M — 
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_ du __ _ dv dw 

7r zx („)=M—,X zy(H) =/J—, 7r yz(H) =/J — 

Normal strains, due to equalities (7) §2 , are equal 
Sxx = e x > S yy = e y , S 7z = e z , 
e x = du/dx,e y = dv I dy,e z = dw/dz 

, 1 ,. -. du 1 _ 3v 

3 ox 3 <3y 

*"«(») =~(P + -VdlVV)+ M , 

5 oz 

where, as opposed to normal strains of Stokes hypothesis, there is 
1/3/jdivv , (effort to justify this formula was made by A. V.Lykov in 
/8/). 

Paradoxes of defining viscous normal stress 
1 st hypothesis on pressure states: arithmetical mean value of 
viscous normal stresses must be equal to pressure (refer /I/): 

KxxJH) +7r yy(H) +7r zz(H) _ ._. 

Due to this hypothesis, into normal stresses in Stokes law (3) 
artificially was included member 2/3/jdivv and into (4) by 
analogy included was member 1/3/jdivv . These artificial 
hypothetical inclusions provide for unconditional execution of 
pressure hypothesis (5). 

It is known that hypothesis (5) for ideal fluids p n = -pit is Pascal 

law, unconditionally executed for incompressible viscous fluids 

divv = , even when normal stresses in them are defined with 

arbitrary coefficient J2 : 

_ dv- 
n a =-p + jU —^ , i = 1,2,3 , (6) 

OX, 



and tangential stress are defined by Newton law n jt = jU — '- ,i ^ j . 



dv^ 
dx. 



j 
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When jTi=0 in (6), 7t u = —p, i = 1,2,3 , result in Pascal law. 

When value of Jx = fj, normal and tangential stress are determined 

by a single representation 

dv. 
n n =-P S ,j+i u ^> l > J = !' 2 ' 3 > ( 7 ) 

OXj 

i.e. there is canonical nonsymmetric tensor 

7t =—pE + /jS (8) 

If during definition of viscous normal strains, to proceed from 1 st 
hypothesis (5), there is an arbitrariness related to selection of 
coefficient /2 . (quite possible that in simulation of turbulent flows 

this circumstance may play an important role) 

2 st hypothesis on pressure states: that normal stresses in 
compressible gas divv ^ must satisfy equation 

7rxxiH)+7r ^" )+7r ^ ) =-p + J u'divv, (9) 

where arbitrary coefficient // (for its selection there are no any 

physical justifications) is declared a coefficient of bulk viscosity or 
second viscosity coefficient 

It is not difficult to calculate that 2 st hypothesis is observed for 
normal strains in the following form 

1 dv. 

^u(h) =-[p + (-ju-ju')divv'\ + ju— -,i = l,2,3, (10) 

3 dx t 

i.e. nonsymmetric strains tensor takes the following form 

1 dv. 

n n(n) =-ip + (-jU-ju')divv]S ij +ju— —,i,j =1,2,3, (11) 

3 ox j 

n H =-[p + (-jU-jU')divv)]E + juS (12) 

As already said, there is a definite arbitrariness in selection of 
multiplier fj . Disposing this arbitrariness, let's assume //'=—//, 

thanks to which normal strains (10) in compressible gases divv ^ 
will take the following form 
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dv 
tt u =-p + ju— -i, i = 1,2,3, (13) 

coinsiding with normal strains in incompressible divv = fluids (7) 

ju = /j, i.e. (12) will transform into (8) and there will be a uniform 

Theological law for both - incompressible fluids, as well as for 

compressible gases: 

dv. 
7r = -pE + juS ,7t„ =-pS ii +ju-^,i,j = \,2,3 (14) 

dXj 

Only hypothesis (9) will take other form 

= -p + -/jdivv , (15) 

realized for normal strains (13)! 

In the event of incompressible fluids with divv = hypothesis (15) 
automatically transforms into Pascal law. 

Thus it was shown: 1) that inclusion of complexes \l3judivv or 
2 / 3/jdivv is some absolutely obvious slanting of viscous normal 
strains to pressure hypothesis (5); 2) formulae of normal strains 
depend on form of pressure hypothesis. 

Refusal of artificial additions of types of 213/Mivv , l/3judivv , 

(2 / 3jU — ju' )divv, (1 / 3ju — //' )divv, considerably simplifies 
strains tensor, and, accordingly, vicous compressible gas 
dynamics equiations with divv ^ . 

Described above paradoxes with pressure hypotheses and Newton 
rheological law (14) allow making a conclusion: proportionality of 
viscous part of normal strains to diagonal elements of displacement 

tensor S , is justified by Newton friction law, due to which tangential 

stress are proportional to nondiagonal elements of tensor S . 

Newton rheological law (14) releases equations of viscous 
compressible gases dynamics in dissipative part from mixed 
derivatives, in connection with this, three-dimensional equations are 
reduced on 12 derivatives, and if to compare with equations of 
Navier-Stokes, then in total sum in new equations 18 derivatives less 
contained in comparison with equations of Navier-Stokes. Besides, 
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presence of mixed derivatives affects type of equations, i.e. for 
quasi-parabolicity creates specific difficulties in researching stability 
of difference schemes for their numerical solution. 
In §5 it was noted that in universal formula 

b-\ 

v(r +8r,t) = v(r,t)-\ [wtv,SP] + S b SP 

b 
when b=l, there is matrix Si = S , therefore transfer to 

displacement tensor S is performed when b=l 

v(r +<5F,r) = v(r,r) + S Sr , 

thus, it is proved that for formation of Newton strains tensor there is 
no special need in the first Helmholtz theorem, it is necessary to 
proceed directly from the Taylor series 

v(r + Sr,t) =v(r,t) + SSr. 

Division of Taylor series into symmetric S Sr and antisymmetric 
l/2[r<9?v,c5F] parts is directly connected with formation of 
symmetric continuum strains tensor and deformations theory (this 
issue is considered here in Chapter 9 ). 

L. D. Landau, E. M. Lifshitz in /10/ laid down the following 
conditions as additional requirement to rheological law of 
continuum: "... in viscous fluid rotating about axis like a solid body 
with uniform angular velocity all tangential stress must reduce to 
zero". Strains (6) meet this requirement (refer /8/). Indeed, in 

cylindric system r, s, z III, Newton strains (7) have the 
following forms: 

or r r os oz 

dv 7 _ dv E ju dv z 



71 r7 = U -,7t 7 „ = U , 7Z„ = 

or oz r os 

r A -w -n d V r 

n rr =-lP + (-JU-Ju )divv] + ju^—, 
3 or 

1 1 dv v 

n^ = -[p + (- ju - ju')divv] + ju( — - JL + - L -), 
3 r os r 
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7T = —\_p + (—// — // )dlVV\ + jU 



3" 5z 

,. - 1 <3(rv ) <3v 1 c?v_ 

<iz VV = — H H , 

r Or dz r ds 

for fluid divv = , rotating with uniform angular velocity cjl> , 

velocities are equal ( ex - angle between cb and T~ ): 
v,. = O, v z = O, v^ = |[<2>, F]| = 

= |<2>||F| sin ex = ctfrsin ex, 

Accordingly, all tangential stress are equal to zero 

n rz — 0, n zr — 0, n rs — 0, iig. — 0, n Z£ — 0, n^ — 0. 

§8. Backhground for fallacious conclusion on strain tensor 

symmetry 

Let's consider the system of material points with masses m„ 
i=l,2,...,N and with radius-vectors r„ i=l,2,...,N, upon which 

influence forces F i , i = 1, . . . ,N . Resultant (major) force is equal to 

F*=Z,F r 

i 

Moment of resultant (major) force is equal to M c = [r c , F ], 
where r c = / m^ I / m i radius vector of center of mass, denoted 

/ i 

further simply r . = r ; and resultant (major) moment is equal 

i 
Theorem 1. For the system of material points moment of 
resultant force is not equal in general case to resultant moment: 

M c ^M. 

Indeed, there is the obvious inequality 

besides, the simpliest example, confirming the theorem, is the pair of 
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opposite forces p — —F resultant force here is equal to zero: 

F* =F l +F 2 = 0, therefore M c = 0, 
whereas resultant (major) moment of couple is not equal to zero: 
M=[r 1 ,FJ + [r 2 ,F 2 ] = [r l -r 2 ,F 1 ]*0,(r l -r 2 )#F i (1) 
Theorem 2. For mutually given perpendicular vectors Q _L F there 

exists infinite number of vectors r —X I +y J+Z k, for which 

there is an equality [r ,F] = Q. If mutually given perpendicular 

vectors Q J. r , there is an infinite number of vectors F , for which 
this equality exists. 

Actually, this equation, solved in relation to x , y , z , has 
determinant equal to zero 

y F z -z F y =Q X , z F x -x F z =Q y , x F y -y F x =Q z ,(2) 
where 

F = FJ +F y J + F z k,Q = QJ + Qj + Q z k 

From the equality (r ,Q) = , comes expression xQ x +y Q +zQ=0, 
which is represented in the following form 

yQ y +zQ z =-xQ x 

and is solved along with the first equation of this system (2). In the 
result obtained is the infinite set of solutions 

* = Q x (Q z -xF y ) . = -Q x (xF z +Q y ) 

y F y Q y +F z Q z ' Z ' F y Q y +F z Q z 
The second part of the theorem is proved in the same way. 
Theorem 3. For the given vector Q it is possible to find infinite 

number of vectors r and F , for which there is equality 
[r',F] = Q. 

In this theorem system (2) is solved together with equations 
(r*,Q) = and (F,Q) = , 
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in the result, for 6 given values x ,y ,z ,F x ,F ,F obtained are 5 
equations, i.e. 5 unknowns will be expressed via the sixth, for 
instance as in theorem 2, via x . 

Resultant force is in the right part of dynamics equation (10) §3, 
where with the purpose shortness of strain will be denoted by inferior 
indexes 

71 x = 7t\ , TV y = 71 2 > 7Z z = ^3 > 
Multiplying (10) §3 vectorially by T , we'll find moments 

-, dv _ - ^ S7T ,. 

\r,p— -] = [r,pF + V— -+-] (3) 

at 7~7 ox j 

It is obvious that in (3) moment of major force is equal to 

M c =lr,pF + '£ i —^--\ (4 ) 

7=1 ^ X 7 

Firstly, symmetry of strains tensor in III is derived from theorem on 
moment of momentum change (angular momentum), represented in 
the form of integral 



dt 



j]J[r , vpSr] = j]J[r , FpSr ] + cg[r , ;?„<?] , 



with introduction of differentiation in the left part for integral 

\\\j[7,vpSt}= \\\[7,FpSt] + ^r,n n S\, (5) 

r r a 

which is unjustified, if volume r = x(t) depends on time. As for the 
left part of this expression, due to formula (7) §1, there is inequality 

-lp,vpSz]= \l\-[?,vpfc] + §r,vp\(v,n)S<y* jjj-[r,vpSr] 



From thus formulated theorem in III, in the same way, in other 
textbooks, derived was wrong in general case relation 

dv - 3 d 

[r,p—] = [r,pF-\ + Y—-[r,^j] (6) 

dt j~l ox j 
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Mismatch of expression (6) to fundamental theorem on change of 
moment of momentum of system of particles is considered below in 
§9. In the right part of (6) there is a resultant (major) moment 

- 3 d 

M =[r,pF] + ^-— [r,^] , (7) 

j= i ox j 

since (7) is derived from relation 

~7 j]J I' ^ T J j]Jl/' P^ St + $/» S ^ cos ("€*j W a . (8) 

where in the right part, there is a summary (in the form of triple 
integral) major moment of bulk forces 1 1 1 [/ , pF]Sr and major 

r 
3 

moment of surface forces jj^'^^j C0S(ftgX ; )]<x7, obviously, 



o 



j=l 



having different radius -vectors T~ . In fact, from (8) instead of 
formula (6), according to (9) §1 

derived is exact expression 

d _ _ 

— [r,vp] + d/v([r,vp]w)j +d/v([r,v/7]v)7 + 

+ J?v([r,v/?]w)fc + -g =[r,pF] + > -— [r,^-.], 

consequently, left part of (6) is not equal to left part of (9), and from 
(9) in no way results fact of symmetry of strain tensor in the general 
case. 

By differentiation (7) in III is represented as follows 

M=[r,pF + £-^] + £[--,^], (io) 
;=i 0x j i=i 0x i 
Due to (4) major moment (10) is represented in the following form 
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M=M c +Y t l--,Z J -\, (11) 

By theorem 1 these moments are not equal to each other: M & M . 

In III and other textbooks equality of these moments is 
traditionally assumed, which is possible only when the last member 
in (1 1) is equal to zero: 

t[f-,*,] = o (12) 

For the moment of resultant forces M be equal to resultant 

c 

moment ]VZ condition (12) must be observed in each point of 

continuum. 

Theorem 4. Condition (12) is observed in the following three 

situations'. 

1. Strains are parallel to coordinate axes: 

dr 
tf.-H-— ,; = 1,2,3 

dx j 

2. Strain components are symmetrical: 

n^n. .,/,; =1,2,3. d3) 

Thus, symmetry of strain tensor (13) and resulting from it 
symmetry of strains contradict to theorem 1, which is confirmed by 
given in §6 examples 1-8 and list of paradoxes. This implyies 
logical conclusion that symmetry of strain tensor by Stokes 
hypothesis is not physical property of continuum, but is the 
consequence of artificial equating of major moment M to the 
moment of major force J\/f , i.e. mathematical condition of 

c 

equality of these moments in the sense that if in some point of 
flow conditions of (13) are observed simultaneously, the major 
moment of force becomes equal to the major moment 

M c = M , if these conditions (13) are not observed, they are not 

equal to each other: M c ^ M . Consequently, theorem on change 

54 



of moment of momentum was incorrectly applied to arbitrary volume 
of continuum, which was stated in §3, because this equating occurred 
due to equating left parts of expressions (3) and (6). Strain tensor is 
not symmetric in the general case, nonsymmetry for individual flows 
was indicated in /3/. The following conclusion is just: for viscous 
media, Newton theological law with nonsymmtric tangential stress 
§7 is just, since Stokes hypothesis leads to the above-mentioned 
contradictions. 

§9. Momentum theorem does not mean symmetry of 
strain tensor. Strain tensor is nonsymmetric 

For a system of material points with masses m t , moving with 

velocities v. and affected by forces f i , momentum thorem has the 

following form 151: 

d N N 

-£[^m.v.]=£[^,/;.] (i) 

For taken as a parallelepiped volume fa = SxSySz °f individual 
particles of continuum, radius-vector r m moves towards the point of 

application of resultant bulk force Fdm , accordingly, radius- 
vectors of resultant surface forces are directed towards points of their 
action on each face of parallelipiped: 

Due to this, according to fundamental theorem (1) for fa — SxctySz 
major moment of forces will be presented in the following form 
M* = [r m ,SnF] + [r^, 7T x2 fy&] +[r^, ti_ xl fy&] + 

+ [>*> KyT.&Sx} + [^> 7T_ yX dzdx\ + 
+ [r t ^,7r z2 3cSy] + [r t ^,7r_ zl ScSy] (2) 

And resultant force is in the right part (3) §4 and is equal to 
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It is obvious that by theorem 1 major moment of force is not 

equal to major moment of forces: If, F pe3 \ ^ M . Let's 

consider in more details the integrand in the left part (8) §8 

— I' v pdt _, equal to |, Snv _ and representing moment of 
dt JJJ 

T 

momentum Snv of mass Sn = pSx , where OT is individual volume 
containing mass Sn = 2_ t m i , at this each particle has velocity V ( - , 

i.e. there is a conformity 7J — > m ; - — > V ; -,Z = 1,2,..., due to which the 
major moment of system of particles vectors in volume OT will be 
equal to /TJ" , m^. ] . According to theorem 1 it is not equal to 

moment of major momentum [r , Snv] : 

^[/;,m I .v 1 .]*[r,<Smr], 

i 

because 



[,<5mv = 






777 .. 



accordingly, equality between them is possible only when conditions 
of theorems 2 or 3 are observed. In passage to the limit r in the 

considered expression, [r,$wv] is a specific radius-vector of 

dv 

continuum point, meaning that by theorem 3 , will be the 

dt 

arbitrary vector, therefore equality between JjJ/>/0 — ]St , which 
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was obtained by integration in the total volume of continuum for 
both parts (3) §8 

and III — \,vpoT _ is impossible, thus the following inequality is 
justified 

T T 

since by formula (7) §1, left part of this expression is equal to 
~7 IW^^pSt] = JjJ— [r,vpSr] + j^[r,vp](v,n)Sa , 

T T <T 

and the right part transforms into expression 

f[[ — [r,vSn] = \\\{ — ,v$n] + [r,cm — ] + [r,v ]} = 

JJJ dt JJJ dt dt dt 

T T 

= §([^*4] +[ ?,Ai 

J f J dt dt 

Absolutely similar conclusions need to be made for moment of bulk 
forces [r,SmF], because F is the resultant of all forces F. in 
volume St , acting upon particles m t of continuum, however 

&n = J^m, , F = X l *,-^ / Z m .- = Ti m fi /dm 

i i i i 

By theorem 1 

[r,3nF]^^[r i ,m i F i ], 

i 

therefore, for a specific vector r by theorem 3, equality 

[r,SmF*]=Y J [r i ,m i F] 

i 

is possible only for such vector F , which won't be in general case 
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equal to true F : F ^ F , i.e. in the right part of (7) §8 there is 

in fact F , which is absolutely arbitrary vector. Classic formula 
from III 

— jjj[r, vSn] = j]J[r , FSn] + jfir, TT n 8a\ (3) 

due to reasons given above does not correspond to fundamental 
theorem on change of moment of momentum (1). 
Indeed, moment of momentum theorem 

d N N 

— Y l \r i ,m i v^ = Y J \r i ,f i ^ 

must be represented in individual volume St in the following form 



dt 



a 5l k 



accordingly, for volume T in the form of triple integral 

at ; ; h 

T ' T ' T <T Sz K 

According to theorems 2, 3, equality sign cannot be put between 
respective members of expressions (3) and (4), i.e. there are obvious 
inequalitites 

- \p,vsn] tJll-^MiipJw * H&inMi 

<$[r, n n ]Sa * j]}<0X , n nk o k ] , (5) 

emphasizing difference of expressions (3) and (4). We need to take 
into account the circumstance that moving coordinates, having 

beginning in the center of mass dm = 2_J n i > accordingly, moving 

i 

with acceleration dv/dt =£ , will be noninertial system. Here, in 
(5), as previously 7t n 3(7 = / 1 7t nk <7 k - is resultant surface force, 
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ing upon S<J = 2_.<J k , at this 0(7 E (7%. , <j£ —site of surface 
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S<J , which is occupied by particle T^Z^ , that is affected by strain 
ft n l , vectors 'ok are directed to particles ^ok , laying on surface 

G dt , QED. 

Symmetry of strain tensor does not result from (4), due to 
indicated above circumstances, i.e. 7t t - ^ 71 „, j ^i,, since in both 

parts (4) there are major moments, by theorem 1 not equal to 
moments of resultant vecors in (3). 

Let's once again show in detail fallacy of equality (6) §8, widely 
used in /l/,/2/,/3/,/4/ for proving of strain tensor symmetry. Let's 
represent it in the following form 

[r,p^-]-[r,pF]-y^[r,7fj] = (6) 

at j~[ oxj 

Symmetry of continuum strain tensor is obtained theoretically from 
this expression and, as a result, slanting of symmetry of viscous fluid 
tangential stress in Stokes hypothesis III. 

Let's prove that in fact left part of (6) is not equal to zero: 

[r,/A-[r f /^]-Y^[r f ;?,]*() 

dt ~~^ dxj 

With this purpose we'll engage precise formulation of moment of 
momentum change theorem in elementary volume St of continuum: 



d_ 



■2[^w / v i ] = ^[^m, J F;.] + c| , ^[r o ,,^crJ, (7) 



c. 



Sz 



Between vectors T in (86) and r t , r A in (7) there is an association 

n =r+r i \r ak =r+r ak \ 
wherefrom, an association between velocities results 
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dr. dr dr.' _ _ 

— = — + — , v. =v +v,.' 

J? d? J? 

Substituting them into (7), we'll find 

dtt i 

+ ^[F + ^',^o-J (8) 

Let's make necessary transformations 

— y[r+r i :,m i (v+v i ')] = — [r,Ym i v] + — [r,V m i v,'] + 
<#T dt T dt i 

Here, dm = 2_, m i = P^? was taken into account in volume £r . 

i 

Having in mind = Sxdivv , let's make the following 

dt 

transformations 



d ^ ^ r Jr Jv ^ ddm 

— [r,omv\ = [ — ,omv\ + [r,an — J + L*% vj = 

dt dt dt dt 

r dr <iv r _ d(pSr) ^ n Jv 1 

= [— ,pSzv] + [r,pSr—-] + [r,—^ v] = [r,pSr— ] + 

dt dt dt dt 

+ [r,(p-— + £r-Rv] = [r,p£r— ] + 
a? a? a? 

+ [f , (pSrdivv + St — )v] = [r , /^r — ] 
J? dt 
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dr _ _ 

Here equality [ , porv ] =0 was used and also equality to zero of 

dt 

continuity equation: pSrdivv + St = . 

dt 

Finally we have 

— V | + jj.,in | (v+v/ =[r,pSr— ] + 

dti ~ dt 

+— [?,^m,v.'] + -^[^,m,v] + -^[^,m i v ! '] (9) 

dt i dti dtj 

Further on, in the right part (8) 

^[r + r i \m i FJ = [r,J^m i F i ] + J^[r i ',m i F i ] = 

i i i 

= [r,F3n] + Y J [r i \m i F i ] = [r,FpST] + Y J [r i \m i F i ], 

i i 

o* r k 



Oft k °Sr k 



+ 



°Sz k 



where by theorem of Ostrovsky-Gauss and theorem on mean for 
elementary volume St there is 



<?Sz k <?Sz 






In the result of these transformations (8) will take the following 
form 

d-V d _ •cp J-r -1 - - fi?\-i- 

[r,p£r— ] +— [r,> m,.v,. ] + — > [r i ,m i v] + — > [^ ,w.v ; ] = 
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J J— 1 OX : I 

Be dividing this expression by St and having done regrouping, 

we'll obtain 

[r,p—]-[r,pF]-^—[r,Jtj] = 

dt ~T X ox. 

Right part of (10) is not equal to zero 

- — {— [r,^m.v.'] + — £[*}', m,.v] + — ^',m.v/]} + 
dr at i dtj at; 

quod erat demonstrandum. 

Accordingly, equation (6), which is used in textbooks of 
Loitsyansky III and Sedov 111 as theorem on change of moment of 
momentum is fallacious. Obviously from (10), that in formula (6) 
Loitsyansky III and Sedov 121 there can be no equality sign. In fact, 
due to (10) there is inequality: 

[r,p—]*[r,pF] + Y—- [r,7iA 

dt j-fdXj J 

It is widely known that theoretical conclusion about symmetry of 
strain tensor is established on the basis of equality (6). Since (6) does 
not apply due to (10), there is no symmetry of strains tensor in 
continuous medium in the general case, quod erat 
demonstrandum. 
Note. In equilibrium, velocities of all particles of fluid are equal to 
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zero v. =0 Vi , v = 0, — = , continuum equilibrium 

dt 

equation 

pF + ^~^ = 0, (12) 

M dxj 

is equality to zero of major force (see §8), besides, the expression 
(11) will be equal to zero 

^It^E^I+tIR'-M+tI l',"vVj=o (13) 

or dt i dti dti 

By fundamental determination of moment of momentum (7) will be 
equal to zero in the right part (10) and expression 

^(T J [r i ',m i F i ] + ^Y J [r ak ',;r nk a k ]) = (14) 

due to equalities (13) and (14) in equilibrium, there will be equiality 
of major moment of foreces to zero 

M=[r,pF] + f j ^[r,x j ] = (15) 

y=l OX j 

From (12) results equality to zero of moment of major force: 

M c =[r,pF + Y—^] = (16) 

By theorem 1 §8 in the general case they are not equal to each other 

by definition: M c =£ M . Their equality M = M takes place 

under conditions of theorem 4: 

/. Strains are parallel to coordinate axes: 

dr 

dXj 

2. Strains components are symmetrical: 

71.. =71.., i,j = 1,2,3. (17) 

This implies logical derivation that symmetry of strain tensor of 
continuum even in equilibrium is not a physical property of 
medium, but is the result of artificial equating of the major 
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^■11— j = 1,2,3, 



moment M to moment of major force ]\/[ , i.e. mathematical 

c 

condition of equality of these moments in the sense that if in some 
point of flow conditions (17) are observed simultaneously, the 
moment of major force becomes equal to major moment 

M c = M , if these conditions (17) are not observed, they are not 
equal to each other: M c ^ M . 

§10. Paradoxal application of theorems on angular momentum 

(momentum) variation. Fallacy of the Navier-Stokes equiations 

Let's remember the following logical order of deriving moment 

of momentum change theorem. It is very important here. In 

particular, that from law of conservation of momentum of individual 

material points (see 111) 

d(m,.Vj) 



dt 



= F,+ Y J F lk , i=l,...,N ( i) 



k=\k*i 

results, firstly, law of conservation of material points system 
momentum (see 111) 

WK-v,) v v v 



dt 



=£ F < + £IX' *=u,ff, (2) 



i=l Ul i=l i=\ k=lk*i 

where sum of internal forces is equal to zero 

N N 

therefore, there is theorem on change of momentum for the system of 
material points 

d N N - 

-t5/W=£f„ i=l,...,N, (3 ) 

dt j=i ;=i 

secondly, (1) results in moment of momentum for each individual 
material point 

at 



k=\Ui 
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summing up of which results for the system of material points in 
equality 

^£[^m ( vJ = £[^] + £ fff^l i = l,...,N, 

where sum of moments of internal forces is equal to zero 

A' A' 

;=1 k=lk*i 

Due to this theorem on change of moment of momentum in closed 
system of material points above is formulated as follows 

d N N 

"' (=1 i=\ 

N 

where / ,[^,^] - sum of moments of external forces. 

!=1 

Basing on (5), precise formulation of theorem on change of 
moment of momentum in lementary volume St of continuum was 
given above: 

d N N _ - cc 

— Z^' m ^] = Z^' m ^] + §Z^'^^] ( 6 ) 
at , =1 , =1 asi k 

So, (1) results in (4) and then (5), but not vice versa, i.e. (5) does not 
result in (1) explicitly. Let's apply this sequence to deriving theorem 
on change of moment of momentum in continuous medium in other 
formulations. Let's represent continuum dynamics equation under 
strains 

dv - dn x dn dn 7 

p = pF + - + - + ±- (7) 

dt dx dy dz 

As was shown in inductive method §4, (this equation is the prototype 
of equation (3)). Multiplying (7) vectorially by radius-vector r , 
we'll derive a prototype of moments theorem (5): 

at ox dy oz 

From equivalent transformation (8) 
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— [r,pv]-[r,v-^-] = [r,pF] + [r,— - + — - + — ■ -] (9) 

at at ox Oy oz 

results theorem on change of moment of momentum pv , due to a 
unit of volume: 

— [r,pv] = [r,v-^-] + [r,pF] + [r,— - + — - + — -] (10) 
at dt ox oy oz 

For incompressible continuous media p = const and theorem on 

change of moment of momentum has the following form 

-[r,pv] = [r,pF]+[r,— ± + — L + — 1 -], (11) 

at ox oy oz 

fully coinsiding with moment definition (5). Moment of individual 

volume St is equal to pStv = Snv , therefore moment of 

momentum change theorem from (8) is derived for elementary 

volume St in the following form 

„ _ dv -; 3ir x Sir dn 7 ^.„ 

[r, /x?r — - ] = [r, p£rF] + [r , — ^ + — -^ + —HSz 
dt ox oy oz 

or 

_ dv n r ^ - ^ d;r x Sir <3tz- -.- 

For compressible media St = <5r(0 and = Svlivv , 

dt 

dSm 

= and this expression transforms into theorem on change of 

dt 

moment of momentum for individual volume St for compressible 

media 

d r ^ r ^ - r ^ 5^ x d;f dn 

— [r,anv\ = [r,anF] + [r, — - + — - + — -]dr (12) 

dt dx dy dz 

Integration (12) for the total volume r results in expression 
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T T 

f ff r - d7T x V7T d7T 7 , _ 

J ~ J ox ay oz 
in principle different from familiar formula 

Jib" ^' ^ pST ^ = Jlf^' ^*^ + tf F ' ^" S<7 ^ ' ( 1 4) 

r r a 

given in books of Loitsyansky, Sedov and other. Obviously, 
conclusion of symmetry of strain tensor in no way comes from 
correctly formulated theorem on moment of momentum change 

(12). Fallacy of (14) and formulae derived from it was discussed in 
previous paragraphs. 

Well, list above paradoxes and evidences dissymmetry of strain 
tensor of continuum confirm mistaken of Stokes hypothesis. As 
exactly by this hypothesis derived equations of Navier-Stohes, 
really, equation of Stokes 

dv i dp „ ^ 8 r ,dv i dVj 2 5. 
at ox. i^ox. ox, ox. 3 ox. 

1 J J J ' * 

dt 2 ~7~7 ax j dXj 

- pdivv — // 4li vv ^ 

are mistaken. 

Equations of ductile incompressible fluid, set Navier by 
dissymmetry of strain tensor of Newton 71 H , 

dv 
p[ — + (v,V)v] + Vp = juAv +pF, 
8t 

(V,v) = 

are rectilineal. 
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§11. Analogue of Stokes hypothesis. Antisymmetric tensor strains 

§7 set absolutely logical issue: if, according to Stokes hypothesis, 
symmetric strains 

dv. dv,. 
n m = ~(P + 2/3judiw)S fi + //(-^ + —±), i, j = 1,2,3 (l) 

dxj dx i 

J^ 1 dv. dv ■ 
are caused by deformative displacements > — ( — L h )Sx , 

j^2 d Xj dx/ J 

what strains create in Taylor series (in convective acceleration) 

Jn 1 3v. dv,. 1 dv. dv,. 
"~I 2 ca^. fix. 2 Sx y . ox i 

rotational displacements — \rotv,Sr = V— ( — ! -)&•? 

2 ^2 fix,, fix,. 

An alternative hypothesis is put forward to answer this question by 
analogy with Stokes hypothsis: let tangential stress be proportional 
to components of rotational displacement, i.e. there are 

1 dv. dv. 
antisymmetric (skew-symmetric) strains n ji =—fl{ — ), 

J l 

accordingly, alternative tensor will be equal to 

x m) =-p E + 2 v£> 

dv- dv, 

^m = ~P s ij + M^r - ^r)'*' J = l > 2 > 3 ( 3 ) 

J * 

Substituting (3) into equiations of dynamics in strains, we obtain 
alternative equations of viscous flud dynamics: 

dv. dp ^ d dv. dv,. 

a? ax, l^i dx. ax,. ax 

Normal strains 7T ii{Alt) = —p,i = 1,2,3 correspond to Pascal law. The 
following fact are of interest. With constant viscosity jil = const 
equations (4) are represented as follows 
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dv dp F ^ ^d 2 v t a ^ dv, 

at ox i ~[ dXj ox t ~[ ox j 

from which for incompressible fluid without sources and due to 

continuity equation divv = /, — = obtained are Navier 
equations for incompressible fluid 

P"T L + ^ 1 = P^+//2.^f' 1 = 1 ' 2 ' 3 ( fi ) 

The same Navier equations are also derived from Stokes equations 
with symmetric strains (1): 

dv, dp „ ^ 5 r ,dv, oV, 2 5 , ,. _ . ,„„ _. 

«f OX; j=l OX J OX } OX, 3 OX; 

and the same equations (6) are obtained from equations of dynamics 
with nosymmetric Newton strains: 

7T H =-\_p + (jU/3-jU')divv)~\E + juS , (8) 

p dv L + ^ = pF+ j^ JL^^L)- JL[( M /3- M ')divv-],i=l,2,3 

dt dx t H-V x j v x j u x i 

But Stokes hypothesis (1) leads to the above-given paradoxes, 
because, as was proved, continuum strans tensor in the general case 
is nonsymmetric, and in alternative hypothesis equations (4) in the 
general cases of variable viscosity jU ^ const loose property of 

ellipticity. These two circumstances confirm adequacy of 
nonsymmetric Newton strain tensor (8) to flow of viscous fluids. 



§12. New equations of vicous fluids dynamics with Newton 
nonsymmetric strain tensor 

1 ..,.-,„ dv. 

; — '- 

OX, 



Xjhu) = -iP + (~ju- /Li)divv]d fi + /u^,i, j = 1,2,3 



j 

Substituting components of Newton tensor given in §7 in Cartesian 
coordinates, into the general equation in strains and projecting to 
axes of coordinates we derive new equations of viscous fluid 
dynamics (Jakupov K.B. /8/): 
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,du du du du. dp 
p{ — + u — + v — + w — ) + — = 

dt dx dy dz dx 

d , du. d . du. d . du. d ,.1 ,. ,. _ „ 

= —(M—) + — (M—) + — (M—)-—[(-M-M)diw] + pF x , 

dx dx dy dy dz dz dx 3 

,dv dv dv 9v s dp 
p( — + u — +v — + w — ) + — = 

dt dx dy dz dy 

ddvddvddvdl 

= ^^—)+—(M—)+—(M—)-—i(-M-M')divv]+pF 

dx dx dy dy dz dz dy 3 

,dw dw dw dw. dp 
p( — + u — + v — + w — ) + — = 

dt dx dy dz dz 

d dw. d , dw. d , dw. d 1 _ 

= —(//—) + — (//—) + — (//—)- — [(-//-// )<fcw] + pF z 

& & qy qy oz <9z <9z 3 

from the total equation of energy balances thermal conductivity 
equation is derived 

dT dT dT d7\ 

PC V ( \-U hV \-W ) = 

dt dx dy dz 

d dT d dT d dT 1 

= — (A — ) + — U — ) + — (A — ) - pdivv -(-ju-ju' ){divvf + 

dx dx dy dy dz dz 3 

r ,&,2 ,3w. 2 ,du. 2 ,5V^2 
+ //[(—) + (— ) + (— ) + (— ) + 

ox dy dz ox 

dv i dv o dw -> dw ■> dw i 
+ (— ) 2 + (— ) 2 + (— ) 2 + (— ) 2 + (— ) 2 ] 
ay <3z dx dy dz 

In incompressible fluid p = const , divv = : 

,<9« du du du. dp 
p( — + u — + v — + w — ) + — = 

9? ck ay 9z ck 

d . du^ d , du. d , du. „ 
= — (//— ) + — (//— ) + — (//— ) + yC^, 

& & qy oy az <9z 

,9v <3v <3v dv. dp 
p( \-u hv \-w — ) + — = 

dt dx dy dz dy 



(2) 
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d dv d dv d dv 

ox ox dy dy dz dz 

.dw dw dw dw. dp 
p( — + u — + v — + w — ) + — = 

dt dx dy dz dz 

d , dw w d , dw. d , dw. „ 
= — O — ) + — (u — ) + — (M — ) + P F ?> 

ox ox oy oy oz oz 

ffl dT 8T d7\ d ,,dT. d ,,dT. d ,,&T S 

pc v (—+u—+v— + w—) = —(X—) +—(X—) +—(X—) + 

dt ox oy oz ox ox oy oy oz oz 

t ,du. 2 .du. 2 ,du. 2 .dv. 2 

+ //[(—)+(—) +(— ) +(— ) + 
ox dy oz ox 

,dv. 2 ,dv. 2 ,dw. 2 ,dw. 2 ,5w. 2 , ,.. 

+ (— ) +(— ) +(— ) +(— ) +(—)] (4) 

oy oz ox oy oz 

Derivation of Navier equations for incompressible fluids 

^ + fU^. +/ £|^ = l,2,3 ,5, 

at ox i ~r x ox j 

from Stokes equations in constant viscosity 

a? ox ~r ax , ax , ox 

as is known, uses continuity equation tffc'vv = 0/1/- If the flow 
contains discretely located sources or drains with specific power J , 
dix>X> = J will be a discontinuous nondifferentiable function, 
therefore derivation of the indicated equations in the form (5) is 
impossible, while (5) from new equations (3) us derived without the 
use of continuity equation divv = O • 

As opposed to Navier-Stokes equations, equations (1) contain 9 
derivatives less in three-dimensional case and 4 less in two- 
dimensional. Besides, equations (1) are easily represented in 
divergent form: 

-^ + div(pv iV ) + ^ = 

dt dx i 
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d 1 

= pF i + div( / ugradv i ) [(—ju-ju )divv\,i =1,2,3 (6) 

dx t 3 

If in the right part of this equation clivv = Ois included, derived 
is analogous divergent expression of equations (3) of incompressible 
fluid. 

Brief expression of energy balances equation for nonsymmetric 
tensor 

AT 1 

pc v — = div(AgradT) + 
dt 

3 3 dv 1 

+ jliV'Y j (— ^f -pdivv -{-ju- ju'Xdivv) 2 (7) 

in divergent representation is similar with dynamics equation (6) 

c v [ 1- div(pTv)] = div(AgradT) + 

dt 
3 3 q v i 

+ M^Zli-^ 1 -) 2 ~ pdivv - (- ju - ju')(divv) 2 

In new equations of dynamics, attention is attracted by dissipative 
members div(jUgradv .) , i=l,2,3, which in their structure are 
similar with dissipative member in energy balance equation 
div(AgradT) , derived from basic Fourier conduction law. 

Prandtl number rr = connects viscosity JU with heat 

A 

conduction of medium X , which emphasizes uniform essence of 
molecular transport of substances, which in one case is the 
temperature T, and in other cases - components of velocity 
v t , i = 1,2,3 or concentration C m according to Fick law. 
Note. Due to arbitrariness of choise of second coefficient of 

viscosity, if ju'=—ju is included, the given equations will be 

reduced by 12 derivatives. 
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§13. New equations of vicous fluids dynamics with Newton 

nonsymmetric strain tensor in cylindrical coordinates 

n H =—[/? + (1/ 3// — //' )divv]E + juS 

Deriving similar equations in cylindrical and spherical 

coordinate systems does not result in any difficulties, the only thing 

here is to take into account nonsymmtry of strains in these systems 

(Jakupov K.B. /8/). 

In cylindrical system for constant viscosity ^/ = const and 

density, dynamics equations coinside with known Navier-Stokes 
equations 131, in the case of variable viscosity and density, they are 
derived by substituting with given in §7 nonsymmetric strains 



re r \ ~ ' ' er —■>&•** 

or r r as oz 

dv 7 dv F u dv 7 

or oz r os 

r A .w -n ^ V - 

n„ = ~VP + {- M - M )divv] + ju 



3 dr 

&* =-\.P + (-M-M) dlvv ] + tf. — r J - + — ]' 
3 r as r 

3 oz 

in continuum dynamics equation (Lykov A. V. 131): 

dv dv v dv dv v 1 d(rit ) 
p(—+v r — + ^ —+v z — — s -) = - — + 
dt dr r ds dz r r dr 

1 Stt d7V n^ 

r ds dz r 

,dv F dv r v di> dv.. vv 

p( — - + v r — - + — — - + v, — - + - J - L ) = 
dt dr r ds dz r 
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r dr r ds dz 

.dv 7 dv 7 v dv 7 dv 7 . 

p( — - + v r — - + — — - + v z — -) = 

dt dr r ds dz 

1 d(m ) 1 dn„ dn 22 

= — - — — + - + — - + pF z 

r dr r ds dz 

energy balance equation takes the following form 

dT dT v cT dT id gr. 

pc v { + v r + -^ + v 7 ) = (Xr — ) + 

dt dr r ds dz r dr dr 



1 8 ..3T. d ,,57\ 

— — (X — ) +— (X — ) 
r ds ds dz dz 



^ ^^ f ^2 A dV F V r . 2 ,dV T . 2 ,dV r . 2 

® v = MM + (" T 1 + -) + (^) + (^-) + 

dr r ds r dz dz 

Sv 1 l dv -> dv v i 
+ (— ) + ( — — ) + (— - — ) + 
dr r ds dr r 

dv o 1 dv 9 _ 1 ^ o 

+ (—^) 2 +(--^) 2 ]-pdivv-(- J Li-M')(divv) 2 
dz r as 3 

Note. Due to arbitrariness of selection of the second coefficient of 

viscosity, if //'=/// 3 is included, the given equations will be 

reduced not less than by 12 derivatives. 

§14. New equations of vicous fluids dynamics with Newton's 
nonsymmetric strain tensor in spherical coordinates 

n H = —[p + (1 / 3// - //' )divv ~\E + juS 

Components of nonsymmetric strain tensor according to Newton 
law have the following form (Jakupov K.B./& /): 

1 dv E sine 5 v, dv E v 

rsms d(p r ds sin s dr r 

Li dv 1 dv 

7t„=-—^,7i rr =-[p + (-pi-pi)divv] + pi—^, 

r ds 3 dr 
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r A in j- -n A dv„ V r . 

Kse =-[p + (-M-M )divv] + jU( — — ^ + — ), 
3 r os r 

1 dv r d v, 

^ = ;" — — .^ =Mr — ( —), 

rsin ^ <90 or r 

, 1 ^ v r v £ c?g£ 

+ M— — ~^T + — + — )' 

rsin 6" c>0 r r 

,. - 1 d(r 2 v r ) 1 5v> 1 d(v_sins) 

dlVV = — — H H , 

r dr rsin^ d(p rsin^ ds 

are substituted into equations of dynamics in strains (Lykov A.V. I'M): 

2 2 

,dv r dv r v F dv r v, dv r v e +v ^ 
dt dr r ds rsms d<fi r 

1 d(r 2 7T rr ) 1 SO-,, sin ^) 



r 2 dr rsin s ds 

1 5 ^ r Kse + ^ r 

+ — + pF r , 

rsins d(j) r 

.dv E dv s v s dv s v, dv s v r v e vjctgs 

p( — £■ + v r — - + — — - + — + - L - L - — ) = 

dt dr r ds rsms d(j) r r 

1 d(r 2 n rF ) 1 d(7t„srn.s) 



re' j ■ bi. '_ I 



r dr rsins ds 

1 d7Tji„ 71 + 71 ctss 

+ — ^f + ^ --—*„ + PF S > 

rsms dtp r r 

/?( — -+v r — v - + -^- — v - + — f- - + ^— + — -ctgs) - 

dt dr r ds rsins d<f> r r 
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1 d(r n H ) | 1 g^ | 1 g^ | 



r 2 3r r 3s - rsins - 3^ 

1*^+7*4, IctgS 

r r 

in energy balance equation the following is included for 
nonsymmetric tensor 

3r r ds r rsing- d<f) r r 

9r r r 9e r 5f sine r sin e 30 

+ (^—-j) 2 +[r — (^)f}-pdivv-(- J u- J u')(divv) 2 
r sin £• o<p or r 3 

Afote. Due to arbitrariness of selection of the second coefficient of 
viscosity, if ju'=—ju is included, the given equations will be 
reduced not less than by 12 derivatives. 
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Chapter 2. FALLACY OF PRANDTL EQUATIONS IN 
BOUNDARY-LAYER THEORY 

Fallacy of theory proposed by Prandtl in 1904, states that in 
system of equations of boundary-layer derived from reductions III, 
121 carried by him, second Newton law and law of conservation of 
mass are not observed. 

§1. Velocity vector components are determined with the help of 
Newton's second law or momentum theorem 

This is the truth that became trivial, but here appeared necessity in 
its application for proving statements given in the chapter title. Let's 
use the following formula of Second Newton law for a material point 

with mass m=const, that is under resultant force / 



dv - 
m — = j 
dt 
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For force / , not depending on velocity v = ui + vj + wk , 
computation of all velocity components is performed with the help of 
integration of second law on time slice [? , t] 

t 
v\,-v \ t + \fdtlm, 

which gives in projections 

/ t t 

u\ t -u\ t + \f x dt/m, v| f =vl f + \f y dt/m, w\ t -w\ t + \f x dt/m, 

h h 'o 

In case of force / , depending on v , integration becomes more 

complicated, which in general takes place while computation of 
flows of viscous fluid. 
Let's consider Navier equation of viscous incompressible fluid 

p — = -Vp + //Av + pF, 
dt 

having multiplied it by elementary volume St , we'll come to 

formulation of Newton second law for mass dm = pSx : 

dm — = / , 
dt 

where / = (— Vp + juAv + pF)5m is the major (resultant) force. 

It is obvious that according to Newton second law, velocity vector 
v , and, consequently, all its components must be computed only 

from the major equation of dynamics dm — = / , and not from 

dt 

equation of law of conservation of mass divv = , as is the case 
in Prandtl equations. 

The purpose, for which law of conservation of mass is used in the 
form of continuity equation, is considerd below. 

§2. Action of law of mass conservation upon flow is performed 
via pressure only 

Let's consider equation of viscous incompressible fluid dynamics 
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dv _ _ 1 _ - 

— + (v,V)v+-V/? = vAv+F, (1) 

along with continuity equation 

(V,v) = (2) 

In §1 it was showed that computation of velocity components must 
be carried out only through solving of Navier equations-(l) in 
projections, accordingly, logically, pressure p must be computed 
only from equation of law of conservation of mass, in other words, 
from continuity equation (2), which, properly speaking, is realized 
in grid solution methods of initially -boundary value problems for the 
system (1), (2) (see 131). 

Technique of deriving equation determining pressure is known for 
long and nalitically is in the following. Differentiating continuity 
equation in time t, nedded expression is derived 

dv 
(V,— )=0 (3) 

at 

In the result of substituting into (3) a time derivative 

dv ~ 1 

— = vAv + F - (v, V)v Vp 

dt p 

Found is desired analytical expression for pressure, realizing 
fundamental connection between pressure and velocity 

Ap + (V,(v,V)v)-(V,pF)-(V,/zAv) = (4) 

It is obtained in natural manner from law of conservation of mass. 

Equation (4) is apparently an elliptical equation in relation to 
pressure, which thus emphasizes diffusive character of pressure 
distribution in continuous medium. In Cartesian system, due to 
equation (2) it is significantly reduced becase equality to zero of the 
last member in the left part (4) (V, juAv) = juA(V, v) = 
accordingly, there is a more compact equation 

AP + (V,(v,V)v)-(V,pF) = 0, (5) 

where there are no forces of medium viscosity, i.e. this equation is 
equally just for both - viscous and ideal fluids. 

Similar association is present in compressible media as well. 
From divergent dynamics equations /3/ 
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fWi + div{pv.v) + -2- = (6) 

dt dx i 

d 1 
= pF i + div{jugradv i ) [(— // — //' )dzvv ] , i = 1,2,3 , 

ax. 3 

c v [ h g?/v(/?Tv)] = div(AgradT) + 

5? 



3 3 5v 1 

;_1 ; = i OX; J 



are singled out from (6) time derivatives 

a i 

+ div{jugradv . ) [(- //-//' ) J/vv] , i = 1,2,3 (7) 

a^,. 3 

and are substituted into differentiated by ? equations of law of 
conservation of mass 

and because of which hyperbolic-type equation is derived for 
pressure 

— r(- J —) = > T— F + — 5-[(-//-//)Aw] + 

a* 2 /?r y V dxf 3^ 

5 
H [div(pv i v ) - div(/ugradv . ) - pF. ] , (9) 

dx t 

i.e. in gases, pressure distribution has a wave-like behavior. In 
numerical methods of solution of initially -boundary value problems 
for systems (1),(2) or (6) computed pressures, as a rule, are carried 
out by difference analogues of equation (4) or equation (9), derived 
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at* 

= (Vp,n)| =({-p[— + (v,V)v] + J uAv+pF} 

a dt 



,n), 



from difference analogues of continuity equations (See 131). At this, 
immediate approximation of equation (4) by homogenous difference 
schemes with the use in the capacity of boundary conditions on 
boundary (j of dynamics equations (1): 

dp 

dn 

as was shown by numerical experiments, is not an effective method 
for solution of initially-boundary value problems for system (1), (2) 
(iteration algorithms for computation of pressure will be vary). 
Absolutely similar situation takes place also for gas dynamics 
equations (6). Therefore, in numerical solution of systems of 
equations (1), (2) or (6) used is a bit different computational 
technique of fundamental associations relatiozations (4) or (9). At 
this substantial is that depending on net domain nodes location, 
difference analogues of associations (4) or (9) are derived in the 
form of closed system of inhomogeneous difference equations for 
pressure, which are easily realized by suitable iteration algorithm, for 
instance, simple but convergent parametric iteration method /3/. 

§3. Fallacy of Prandtl equations is in incomplete observance of 

second Newton law and due to absence of action of law of 

conservation of mass upon flow via pressure 

Exact equations, describing viscous flows, including those in 
boundary layer are undoubtfully equations (1),(2) for incompressible 
fluids and equations (6) for gases. 
Ludwig von Prandtl supposed that in boundary layer near rigid 

surface there are estimates y « vv , v « v v , u « 1, x « 1 , basing 
on which he discarded in equations (1) for plane flows of viscous 
fluid 

du du du I dp ,d u d u. 

— + u — + v — + — - = v(— T + — T ), (10) 

dt dx dy p 8x dx dy 



dv dv dv 1 dp / d 2 v d 2 v 
— + u — + v — + — — = v( — j + — 
dt dx dy p dy dx dy 



+ u — + v — + —^ = v( irY + —), (11) 
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du dv n „_ 

— + — = 0, (12) 

ox ay 

relatively small members and introduced simplified system 111 

du du du 1 dp d u 

— + u — + v— + -— = v — -, (13) 

ot dx dy p ox 5y 

-^ = 0, (14) 

poy 

? + ^=° < 15 > 

ox oy 
The main error of Prandtl system of equations (13), (14), (15) 
according to §1, is in neglect of projection of second Newton law to 
y-axis, which is confirmed by equality to zero in (14), i.e. second 
Newton law is not observed in full, preserved only projection (13) of 
the second Newton law to axis x . 

Obviously, for equations (13), (14), (15) analogue of fundamental 
association (4) will be equation for pressure in the following form 

d p d v d , du du^ d u 

dx oydt 8x dx oy dxdy 
According to above-stated concept of solution of initially-boundary 
value problems for equations (1), (2) or (6), during realization 
considered is system (13), (14), (16). (14) results in dependence of 
pressure on x and t only: p=p(x,t). Bearing this in mind, system (13), 
(14), (16) may be represented in the form of two equations 

du du du I dp d u 

— + u — + v — + — = v — -, (17) 

dt dx dy p dx dy 

d p d v d , du du^ d u 

— f + p— (u — + v— )~ju ^ = (18) 

dx oydt dx dx dy dxdy 

Thus, for 3 unknown functions u,v,p derived are two equations (17), 
(18) only. Obviously, one more equation for computation of lateral 
velocity v is missing. 

Continuity equation (15) may not be used for defining v, since it is 
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derived from equation (18), at this, at the initial time t=0 there must 
be an equality 



,du <9v\ 

(^T + — ) 
ox dy 



= 0. 

t=0 

Accoring to §1 transverse component of velocity v must be 
computed from projection of the second Newton law to axis y (11). 
But in this equation, Prandtl neglected all members associated with 
velocity component v 

dv dv dv ,d 2 v <3 2 v N 

+ U +V V( r-H r-) =0, 

9? dx dy dx dy 

and left equality to zero (14), which was said above. 

Vital error of Prandtl is that according to his concept in the form 
of equations (13),(14) major law of dynamics, in this case second 
Newton law, act in the flow in boundary layer only in 
longitudinal direction along axis x , in crosswise direction along 
axis y do not act, accordingly, in the direction ofy -axis there must 
be no motion according to major law of dynamics , and if there is no 
motion in crosswise direction, this means that there is no 
apprpopriate velocity copponent, i.e. v = . 

If so, Prandtl equations (13),(14),(15) are reduced to form 

du du 1 dp d 2 u I dp _ du 

— + u — + - = v — -, - = 0, — = 

dt dx p dx dy p dy dx 

and final form of Prandtl equations will be as follows 

du 1 dp d 2 u I dp 
— + — = v — =-. — = 

dt p dx dy p dy 

The last equation indicates that pressure will be a function onle x and 
t: p = p(x,t), and also u — u(x,t) . 

In the result, for two unknown functions only one equation. Obtained 
contradiction, is apparently connected with the fact that Prandtl 
system neglects projection of second Newton law to y axis. 

^ dv dv dv dv . . 
Equating to zero of acceleration \-u hv \-w — = U in 

dt dx dy dz 

steady flows is justified by the following considerations: that spatial 
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dv dv dv 
members of transfer u 1- v h u — are negligibly small. Of 

dx dy dz 

course, there can be stationarity, but products in 

dv dv dv 
u h v h w — may be pretty large. Let s snow it on analogue 

dx dy dz 

of only one of products of this expression of transfer, for example, on 

dv ■ , 

v — .Let V = £"Sin CXy, where a = const »1 is a very large 

dy 
but finite number, e — parameter providing for infinitesimality of 

dv 

function v , its derivative is equal to — = aecosory . Therefore 

dy 

dv 2 ■ 

product v — = as smaycosay will be an infinite nonzero 

dy 

number, it is enough to take a = const * £~ . The same will happen 

for the whole class of functions with removable discontinuities. 

Thus, equating to zero of acceleration 

dv dv dv dv _ 

Yu hv \-w — = 

dt dx dy dz 

is not always justified. Unjustified first of all because of 

contradiction with the major law of dynamics with second 

Newton law, which states: effect of force causes acceleration 

which is nonzero. 

Logical demand of strict observance of the second Newton law 
gives instead of fallacious Prandtl system (13),(14),(15) the 
following approximate system: 



du du du 1 dp d u 
+ u — + v — h — = v 



2 ' 



dt dx dy p dx Qy 

dv dv dv 1 dp d v 

— + u — + v — H — =V j 

dt dx dy p dy dy 
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dv 


*tJv~, 


dv 
u — ? 


» 4^, 


dv 
v — 


dt 




dx 




dy 



V, l / TT SJ% ' V ' 



du dv 
— + — = 

dx dy 

This system is also obtained from Navier equations, based on carried 
estimations, according to which all members of the second equation 
of this system of one infinitesimal order: 

dy : 

and neglected transverse pressure gradient of a higher infinitesimal 

1 dp 
order « v , therefore it may be neglected if desired. 

pdy 

In short, neglected was only this pressure gradient, and all the rest 
members of equation (11) are preserved as values of equal 
infinitesimal order, thus, second Newton law is observed partially, 
i.e. in both directions, as opposed to fallacious system of Prandtl 
(13),(14),(15). 

Obvisouly for these equations, analogue of fundamental 
association (4) will be equation for pressure in the following form 

d 2 p d 2 p d . du du^ d 3 u 

—Y + —Y + P — (u — + v — )-ju j + 

dx dy dx dx dy dxdy 

d . dv <3v x d 3 v 
+ p — (u hv — )-ju — - = 

dy dx dy dy 

This equation may be reduced to the following form 

d 2 p d 2 p d . du du^ d , dv dv . _ 

y + Y + yO (U + V ) + p (U + V )=0, 

dx dy dx dx dy dy dx dy 

having in view continuity equation. 

In our system for three unknown functions u,v,p there are 
corresponding three equations, at this second Newton law and law of 
conservation of mass are observed, and also satisfied is fundamental 
association (4), as opposed to Prandtl equations (13) (14) (15). 

Steady-state Prandtl equations 

du du I dp d u 

M ^ + v ^T + ~/ = v ^T' (19) 

dx dy p dx dy 
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I* = 0, (20, 

pdy 

e l + ^=° < 21 > 

dx dy 
have the same deficiencies that are in (13) (14) (15). 

Due to (20) p = p(x) and dynamics equation (19) is 
represented as follows 

du du 1 dp d u 

u — + v — + — = v — - (22) 

dx dy p dx dy 

By analogy with (13) (14) (15) here also for 3 unknown functions 
u,v,p exist only 2 equations (21) and (22) 
equation for computation of lateral v is missing, since according to 
concept (4) in a steady-state case from (22) we can single out 

du I dp d u. .du 

u=(-v- -J- + V—)!— (23) 

dy p dx dy dx 

and, substituting into equation of law of conservation of mass (21), 
hence obtain equation for pressure in the following form 

d , du 1 dp d 2 u „ ,du n dv 

— [(-v - + v— )/— ] + — = (24) 

dx dy p dx dy dx dy 

Equation (24) forms system along with dynamics equation (22) and 
again for 3 unknown functions u,v,p there are only 2 equations (22) 
and (24), at this continuity equation (21) is the result of equations 
(22) and (24). 

As is known 111,121, Prandtl, for computation of pressure p, does 
not use this equation (24) or (18), thus broken is the fundamental 
association between pressure and velocity. Bearing in mind (20) or 
(14), i.e. dependency of pressure only on x and t: p=p(x,t), Prandtl 
considers equation (17) in external incoming flow, considering it as 
moving with velocity 

u= U x (x,t),v=0, 
and obtains for pressure gradient relation 
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1 dp dU x TT dU^ 

y = ^--U^ (25) 

p dx dt dx 

(This Prandtl technique, resulted in (25), resembles action of baron 
Munchhausen, who pulled himself and his horse from bog pulling 
up hair on his head). 

Obviously approach (25) of Prandtl contradicts to equation of 
fundamental association (18). In steady-state flows, 
assuming homogeneity of external flow 

u= C/ 00 =const, v=0, 
Prandtl, from (25) derives equality to zero of pressure gradient 

dp 

— = , consequently, according to Prandtl, in boundary layer 
dx 

pressure turns out to be constant throughout the flow p=const and 

equations (13) (14) (15) are reduced to limit 

du du d u ,„^ 

u — + v— = v — -, (26) 

dx dy dy 

du dv n _„ 

— + — = (27) 

ox dy 

In this system (26),(27) longitudinal velocity u is determined from 
second Newton law (26), transverse velocity v is determined from 
continuity equation (27), which is fully in contradiction to given in 
§1, §2, §3 proofs. 

Fallacy of Prandtl equations (26), (27) was confirmed by 
numerical solution of Navier equations (10), (11), (12) by semi- 
implicit scheme 131, used for computation of longitudinal flow for a 
plate with length L in the area [0 < x < L,0 < y < S] , transferred 



into dimensionless area [0 < x'< 1,0 < y'< 8/ vRe] . Results of 

computations, carried out with number Re = U^L/v = 10 on grid 
A^ x = 20, A^ >; = 20 , showed that in boundary layer of the disc 
pressure, there is function of both variables x,y, pressure gradients 

dp dp 
differ from zero — ^ 0, — ¥= . Let's give table of pressure values 
dx dy 

87 



in three cross sections of the disc l)x'=0.05, 2) x'= 1 — 0.1, 

3)jt'=l-0.05,0< j'<8/VRe: 

1)0.982 0.982 0.982 0.982 0.981 0.980 0.979 0.978 0.978 0.979 0.980 
0.981 0.983 0.985 0.988 0.990 0.994 0.997 0.998 1.000 

2) 0.923 0.923 0.924 0.924 0.925 0.926 0.927 0.928 0.929 0.930 0.932 
0.933 0.935 0.937 0.938 0.940 0.942 0.944 0.946 0.947 

3)0.903 0.903 0.903 0.903 0.903 0.903 0.903 0.903 0.903 0.903 0.903 
0.903 0.904 0.904 0.904 0.904 0.904 0.904 0.904 0.904 

Distribution of longitudinal component of velocity u in cross section 
x'=l, 0<y'<8/Jjte 

0.000 0.090 0.176 0.259 0.339 0.417 0.493 0.566 0.636 0.700 0.760 0.811 
0.855 0.892 0.922 0.946 0.964 0.977 0.986 0.991 1.000 

differs from values obtained from solution of Prandtl equations with 
U K = const 

0.000 0.108 0.216 0.324 0.429 0.530 0.625 0.711 0.784 0.845 0.893 0.929 
0.954 0.971 0.983 0.990 0.994 0.996 0.998 0.999 1.000 

In flows of incompressible fluid due to deceleration of particles on 
solid streamlined surface on external edge of boundary layer the flow 

will not be homogenious U ^ {/ w , V ^ , transverse component of 
velocity will differ from zero but extremum condition will be just 



du 
dy 



= 0. 



y=S 



dv 



= 0. 



y=S 



So, approach of Prandtl to computation of pressure by formula (25), 
and steady-state system itself (26), (27) contradict with the equation 
of fundamental association (16) between pressure and velocity, and 
this means that law of conservation of mass is not observed. 

§4. Laminar boundary layers in compressible flow 

All the above-mentioned paradoxes are present in equations of 
laminar boundary layer in compressible flow III: 
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. du du w dp d . du . 
p(u — + v — ) + — = — (jU ) , 

ax dy dx dy dy 

J_3p =0 d{pu) | d{pv) = Q 

p dy dx dy 

dT dT d ..dT dp ,du. 2 
pc (u — + v — ) = — (A ) + u — + ju( — ) 

dx dy dy dy dx ' dy 

Correct, i.e. not violating the second Newton law and law of 
conservation of mass, will be the following equations 

. du du w dp d . du . 
P(u — + v— ) + ^- = — (//-— ), 

dx dy dx dy dy 

. dv dv^ dp d . dv . d(pu) d(pv) 
p(u — + v — ) + — = — (jU ), ^ + H =0, 

dx dy dy dy dy dx dy 

. dT dT d..dT dp ,du. 2 
pc p (u — + v—) = —(A—) + u^ L + ju(—) 

ox dy dy dy dx dy 

§5. Three-dimensional boundary layers 

Obviously equations of three-dimensional boundary layer III 

du du du 1 dp d 2 u 

u hv \-w 1 = v — -, 

dx dy dz p dx dy 

dw dw dw 1 dp d 2 w 

u hv \-w 1 = v — — , 

dx dy dz p dz dy 

I dp _ du dv dw n 
— — = 0, — + — + — =0 

p dy dx dy dz 

have the same disadvantages that are peculiar for (13), (14), (15). 
These equations must be replaced with the following ones 

du du du 1 dp d u 
u hv \-w 1 = v- 



dx dy dz p dx dy 1 

dv dv dv I dp d 2 v 

hv hw h = v — - 

dx dy dz p dy dy 



u — + v — + w — + — ^ = V TTT' ( 28 ) 
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dw dw dw 1 dp d 2 w du dv dw _ 

u hv \-w h — = v — ^-, 1 h — = 0, 

dx dy dz p dz dy dx dy dz 

but the most objective will be use of full equations of viscous fluid 
dynamics 121. 

§6. Substitution of Blasius equation 

Thus, pressure gradients are not equal to zero and in 
apprpoximation of Prandtl even in steady-state flow, equations must 
have the form (19) 

du du 1 dp d u __ 

u — + v — + -^ = v^, (29) 

dx dy p dx dy 

-^ = 0, (30) 

pdy 

^ + ^ = (31) 

dx dy 

Differentiated (29) by y 

d u du du d u dv du \ d p d u 

+ + v — T + + — = v — - (32) 



dydx dy dx dy 2 dy dy p dydx dy 
In equation (32) due to Prandtl assumption (30) 

^ = ^ = 0, (33) 

dydx dxdy 

therefore, from (32) derived is third order equation 

d u du du d u dv du d u 

u + + v — - + = v — - (34) 

dydx dy dx dy dy dy dy 

Introduction of flow function 

u = dy//dy,v = -dy/jdx 

allows precise integration of continuity equation (31) and from (34) 

results equation of forth order for y/ : 

dy/ d 3 y/ d 2 y/ d 2 y/ dy/ d 3 y/ d 2 y/ d 2 y/ d 4 y/ 
dy dy 2 dx dy 2 dxdy dx dy 3 dydx dy 2 dy 4 
Let's use self-simulated Blasius variables 
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' U ~ , = f(?])^U^,u = UJ'(rj),v = Urf'-f) p^- 

i w 2 V x 

As the result we arrive to ordinary differential equation of 4 th order 

2f""+ff'"+f'f"=0 (36) 

As opposed to third-order Blasius equation 2f'"+ff"=0. 
Weakness of Blasius equation is absence of pressure gradient: 
dp/dx = 0. In principle, from (36) results equality 

2f'"+ff"= — , but here dp/dx * ! For equation (36) needed 

p dx 

are four boundary conditions, which corresponds to equation (35), 
while for Blasius equation there were three such conditions, which 
means loss of an important physical condition of flow. 
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Chapter 3. ABOUT TURBULENT FLOWS SIMULATION 

Problem of doing Raynolds equations is widely known. Various 
semiempirical theories, major of which are given in monograph III 
are used for practical purposes. These models are very remotely 
related to original equations for averaged functions and contain a 
number of functional inclusions and constants with indefinite 
physical properties and doubtful universality. Many semiempirical 
theories, like mixing length theory (Mischungweg) of Prandtl til, 
simulate averaged turbulent flows in boundary layer by equations 
which lack pressure gradient, thus law of conservation of mass and 
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second Newton law are not observed, which was discussed in detail 
in chapter 2. 

Of course, appropriate results from such theories may be obtained 
only if slanting coefficients. 

§1. Paradox of Reynolds averaging formula 

In Raynolds approach to derivation of equations for averaged in 
time values, actual values of hydrothermodynamic functions / are 

represented in the form of averaged / and fluctuating components 



(+— t+— 

i L, . . -.. i 2 



/ = / + /',/= — jf(x,y,z,T)dT,f'=— jf'(x,y,z,T)dT = 



(1) 



2 2 

While averaging nonlinear equations of Navier-Stokes, energy 
balance, diffusion etc. appear additional, so called Reynolds or 
fluctuating strains, heat and diffusion flows of the following form 



H 

i 2 

MfW } = — [A{f'(x, y, z, r)if/\x, y, z, r) }dr 
t i 



(2) 



t — 

2 



( A — differentiation operator), definition of which creates artificial 
problem of closing Raynolds equations. 
Paradoxicality averaging formulae (1), (2) is in limits of integration 

[t — t" 1 2, t + 1" 1 2] , where t° - averaging time, and t - current time. 
According to Raynolds (1) it comes that preliminary it is necessary to 
know values of functions for future moments in time [t,t + t" /2] 
and only then carry out averagings by (1), which is not real. In 
computing experiments there is always possibility to carry out 
averaging to current moment of time t , using collected in memory 
values of function in previous moments of time: 

1 ' 1 ' 

f = f+f,f=- \f(x,y,z,T)dr,f=- \f\x,y,z,T)dz = Q (3) 

(-(° t-f 

According to averaging (3) there is 
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1 t 

A{/V'} = — \A{f(x,y,z,T)y/'(x,y,z,T)}dr (4) 

t , .o 



§2. Approximate calculation of averaging integrals 

2.1. There are several simple formulae for approximate calculation 

of definite integrals. For example, in averaging by Raynolds (2) 

integral may be approximately calculated by average rectangular 

formula 



;+ — 
i 2 

A{/V'}=— \h{f'{x,y,z,r)y/\x,y,z,r)}dT = 



(5) 



( — 

2 



,o\2 -.2 



= A{f'(x,y,z,tW(x,y,z,t)} + { —?-- Y Mf'(x,y,z,fy(x,y,z,f)}, 

24 at 



t° t° 



t e[t- — ,t + —] 

2 2 

By Lagrange formula 

—A{f'(x,y,z,t*)ij/'(x,y,z,t*)} = —rA{f'(x,y,z,t)ij/'(x,y,z,t)} + 
at at 

a 3 

+ (t' -t)—A{f'(x,y,z,tW(x,y,z,t *)} (6) 

r**e[?,f*] or f"e [**,*] 
Substituting (6) into (5) we have 

A{fY}=A{f'(x,y,z,tW(x,y,z,t)}+ (7) 

+^-|yA{/'(x,y,z,ry(x,y,z,0} + O(^-ir-H) 

Taking into account that period of averaging t° - one of the most 
important characteristics of turbulent flow and is determined in real 
flows as quantity comparatively small in relation to characteristic for 
this process time scale, in formula (7) residual member 
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+ ^^^YMf'(x,y,z,t)iy'(x,y,z,t)}, (8) 



(t°) 2 

0( \t —t\) may be neglected as relatively small quantity of 

24 

third order, having in mind inequality \t —t\<t°. 

In the result from (7) derived is approximation formula 

A{/V'}= A{f\x,y,z,t)iy\x,y,z,t)} + 

(t") 2 d 

24 dt 
having error 

(t°) 2 * (t°) 3 

0(^-\t -t\) = 0(^-) 

24 24 

The same we do for averaged third and higher moments 
A{f'y/\.xp'} = A{f\x,y,z,t)if/\x,y,z,t)..xp\x,y,z,t)} + 

+ ^rzrvMf(x,y,z,t)y/\x,y,z,t)..j(p'(x,y,z,t)} (9) 

24 dt 

2.2. In real averaging (4) for approximate calculation of integral it is 
possible to apply rectangular formula 

i t 

A{/V'} =— \h{f\x,y,z,T)y/\x,y,z,T)}dT = 

t J „ 

t-t 

= A{f'(x,y,z,t)if/'(x,y,z,t)}- 

t° d Af ... ., ,. ... (io) 

-——Mf(x,y,z,t )y/(x,y,z,t )}, 
2 dt 

where t e [t—t°,t]. 

By Lagrange formula 

d * d 

—A{f(x,y,z,t )i//'(x,y,z,t )} = — A{f(x,y,z,t)y/'(x,y,z,t)} + 

dt at 

+ (t t -t) — A{r{x,y,zX t )\ir\x,y,zX t )} (11) 

dt 

t e.[t,t ] wiu t e|7 ,t] 
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Substituting (11) into (10) , we'll find 

i i 

A{/V'} =— \h[f\x,y,z,T)y/\x,y,z,T)}dT = 

t „ 

t-t 

= A{f'(x, y, z,i)y/\x, y, z,t)}- 

t° 3 M, (12) 

- — — A{f'(x,y,z,t)y/'(x,y,z,t)} + 0(—\t -t\) 
2 ot 2 

From (12) approximate formula is derived 

A{/V'} =A{f'(x,y,z,t)iy'(x,y,z,t)}- 

t° d 
- ——A{f'(x,y,z,t)y/'(x,y,z,t)}, (13) 

2 ot 

t° * t° 2 
having error 0{ — \t —t\) = 0{ ) due to inequality 

\t* -t\<t°, since t*e[t-t",t]. 
The same we do for averaged third and higher moments 

A{f'i//'..xp'} = A{f'(x,y,z,t)i//'(x,y,z,t)...<p'(x,y,z,t)}- 

t" d 
-— A{f{x,y,z,t)yj\x,y,z,t)..jp\x,y,z,t)} (14) 

2 at 

Formulae (8), (9), (13), (14) allow closing systems of averaged in 
accordance with Raynolds equations of hydrodynamics with 
involving pulsations equations or equations for second moments 

/>', />', y/>' etc. 

§3. Application of pulsations equiation for closing Reynolds 

equations 

In equations of viscous incompressible fluids 

p[^- + div(v,v^ + ^- = pF i+ juAv t , i = 1,2,3, 

at ox i 

pc p [— + div(Tv)] = AAT,Y-^ = 0, 

Ot f OX: 



95 



dC„ 



\-div(C m v) = D m AC m , m = l,... (15) 



V; 



dt 
= (p.,T\ s =q,C m s =q m , Vi \ t=0 =v°,T 



t=0 



= T\C 



t=0 



=Cn 



let's show use of formula (8). 

Thanks to averaging by Raynolds (1), from (10) derived are 
equations for averaged functions 

dt j dXj j dXj ox t 



pc( — +> v. ) + c D > J — = AAT,y — - = 0, 

p Dt Y dx, P Y dx, iSx t 



^'-E^^ + S^=^c., m =i,..., (W) 



& 



ck, 



dx, 



; =^jL=^,c M L=^ m ,vl =0 =v.°,r 



/=o 



=r,c 



r=0 



=c, 



When subtracting equations (16) in appropriate manner from (15) 
obtained are pulsations equations 

p — -+\p( v i — +v i — + v i ) - X — — + = 

5? y 1 ck^ 9^. ck^ ^ a* ; . 5x. 

dv' 
= pF' i+J uAv' l , i= 1,2,3, S^ = ' ( 17 ) 



pel +> (v, + v' ; + v' ; )]-c„> — = AAT, 

^" dt y J dxj ' dxj ' dx/ pA t 8k j 



«\ +Z(?/ __W,__V, x. 



dt 



dx . 



j 



,-Y-_A--L = ,> iC , 

' a, ' a/ Y a ; " '" 



' ,«' T 1 '! ' /"" ' ' I '0 T 1 ' T'lO /-ft 



1=0 



= c 



, 



Further on, included into (16) and (17) pulsation expressions 
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z 



dpv'j v\ 



• dx j 



c > — , > — ^ are substituted by 



; & y 



; Sx 7 



formula (13), in the result of which closed system of equations is 
obtained 

p(— +y v,— L +y — — — y — j — )+-*-= 

dt y ' dx y dx j 2 j dtdx j fy 



pF i + juAv i , i = 1,2,3 ,%-+■ = 0, 



(18) 



.ar v _ d T ^dv'jT r d 2 v'T 

pc( +> v. h> — > — )=AAT , 

p dt j J dxj Y dx j 2 j dtdx j 



5Q 
dt 



■ + 



dv',C 



v,.c 



£-i i 3,.. Z-( ^„ o ^-i 3*2 



& 



;' ■> 



3jc. 



2 Y 9f3xj 



= D m AC m , 



m m ■ 



/? — -+} (v, — L + v',— ) + — > — + — = 

a? y ax . ; ax . 2 y dtd Xj ax, 

rV 
^pF'^juAV,, * = 1,2,3,^-^ = 0, (19) 

ar ~ _ dv , sf. f v 3Vr 

pcj +> (v ; . + v\ ) + — > —] = AAT\ 

^ p dt j J dxj J dx/ 2Y dtdxj 



dt 



+ Y(v.^ + v'.^) + t -Y dVjCm =D AC , 



dx 



dXj 2 j dtdxj 



<Pn T \s = <l> C m s = <l m > v i L Q =v- jL q = t\c 



V : 






r=0 



(=()" 



=c, 



r 1 

■ ^ 1 



Sum of equations (18) and (19) is Navier, heat conduction and 
diffusion equations. 
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Averaging time t" , organically included into closed system (18) 
and (19) along with such such physical properties of flow as the 
following coefficients - viscosity // , heat conductivity X , heat 

capacity c , diffusion D m and density p , and also initial and 

boundary conditions for averaged and pulsations must be defined for 
each calculated flow. Numerical solution of closed intertwined 
system (18) and (19) for given // , X, c , p,D m , 

(p { , q>\ ,q,q' ,q m ,q' m and initial time conditions does not represent 

difficulties in kind and may be carried out by difference schemes, 
developed and theoretically justified in 141, 151. 

Equations for second moments v' . v' ■ , v' 7" , v' ■ C may be 

composed with the use of Keller- Friedman method III. Let's limit 
ourselves with consideration of equations for second moments 

T.. = v'. v' • : 

y i j 



dr it ^_ dr. r . dv t dv f dp' , dp' 

ot k dx k t ox k ox k ox i ox j 



k dx k t dx, dx k 



dp' dp' 
Expressing in system (20) second moments v' ,v\ 

dx, dx, 



^dv'.dv' dv\v'v' k 
/ by formula (13), and third moments > by 

k dx k dx k ' k dx k 

formula (14) we'll obtain closed system along with equations (16) 
for v. and equations (17) for v\ , having set appropriate boundary 

T it o and initial conditions r :i ,_ n for r„ : 

IJ \ IJ \ l—\) 'J 

Ot j OX j j OX j ox i 
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Z f8x i 
p[ — L + > (v. — L + v' / — L + v,— )-> — -] + ■ 



a* 



dx ; 



J ax ; 



5x. 



; 5x. dx ; 



^ P F\+juAv\, i= 1,2,3, Z^ = ' 

r 3T g v _ 5r s v av. 5v ^ f 8 . , dp\ 

8t k 8x h - !> -"■ ~- ° - ! ' 






5x 4 5x t 



ax,. 2 a? 7 &,. ' 

^2 



+ v) (v' ^-) = -pY ( ' J * '-^-^-) + juAr iF 

dx. 2 dt ax- k dx k 2 dtdx k 

Sv', dv' t " 8 ,8v\ dv'j 



t " 5x 4 5x 



o ..i 



; = ^ v ; ,=o =v «> v 






2 5? 5xj. cbc t 



(=0 V i» 7 j/|s ^lij^ij 



t=Q T ij 



In transition from (18) and (19) to equations in dimensionless 
variables, having taken as scales of - velocity V* , linear dimensions 

L, time t, =L/V t , pressure p t = pV t , temperature T t , obtained 

are known similarity criteria: numbers of: Strouhal - Sh, Euler - Eu, 
Raynolds - Re, Froude - Fr, Prandtl - Pr, Schmidt - Sc. Along with 



them one more dimensionless complex appears here - Dg= — or 

fV, 

Dg= , defining relation of averaging time to typical time scale. 



§4. Turbulent flows simulation with account to friction caused by 

normal strains 

As experiments show, turbulent flow has pulsating, irregular 
character with intense mixing of layers, which results in multiple 
increase of fluid flow resistance forces III. Basing on these facts, in 
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semiempirical theories, coefficient of turbulent viscosity jU, T is 
introduced, which by several orders exceeds coefficient of molecular 
viscosity // , thus emphasizing role of increase of friction in 

developed turbulent flows. For friction force in laminar flows it is 
enough to determine tangential stress by Newton law with coefficient 
of molecular viscosity // , and neglecting action of normal strains for 
friction (refer to Chapter 1). In turbulent flow friction increases. In 
this connection the following hypothesis is suggested (Jakupov K.B. 
IAI): increase of friction forces in turbulent flows is connected with 
effect of action of normal strains upon tangential stress as a result 
of fluctuating mixing of layers that are moving irregularly, i.e. 
arises necessity to use known formula of friction of two contacting 
surfaces of continuums F mp =k mp F n . 

In laminar flows, in this formula, dimensionless friction coefficient 
"k" if equaled to zero: "k=0". 
This formula 

F =k F 

mp mp n ' 

where F n - force acting normally to friction surfaces, in transition to 
tangential stress is transformed as follows. 

Let strain 7l l —7l ii l-\-7l--]-\-7l & K act upon elemental area 
8x j 8x k , lying on coordinate plane ox j x k , j ^ k ^i, . According to 
above hypothesis, with the purpose of taking account of action of 

normal stress A a in tangential stress ^y > ^n along with friction by 
Newton law, the indicated strains are represented in the form of sum 

n ij ~ %h) + %T) . Xik = ^ik(u) + Xik(T) » i * h * * k, (21) 

where ^ij{n) ' ^ik(u) are tangential stress of molecular 

transport defined by nonsymmetric formulae of Newton friction law 
in Chapter 1, and tangential stress 

fly(T)i fljk(T) are contributions of action of normal stress and are 
calculated by friction formula F mp = k mp F n as follows. Let friction 

100 



force F mp on area &c j dx k has projections ^ mpij ' * 'mpik for 
directions of orts j ,k : 

F = F 1 + F k 

mp mpij J mpik ' 

having divided by elemental area Sx -8x k , we'll find tangential stress 
F F 

mpij mpik 

%n = ^T^T ' n ik(T) = ~ ~ i * j, i * k (22) 



3c j 3c k 



3c j 3c k 



Meaning that in this case, force acting normally is equal to 

F n = \^ii^j^k\ = M&j&k , 
friction formula is represented in the following form 



nip 



= kF n =k\7r ii \Sx j 3c k 



Component of velocity vector 

v =v i i+v j j+v k k 

in area SXjSx k is vector Vj k = VjJ ^V k K . 

Let vector "ft make angles a and a k accordingly with orts j,k , 
then cosines of these angles are obviously equal 



C0S<2 ; =Vj I 



jk 



cosa k =v kl 



jk 



(23) 



V 



jk 



I 2 2 £; 

~\ V j k ■ But friction force /*„,„ is directed in 



mp 



opposite direction from velocity " jk : " jk T4- F , therefore, 
there are formulae for this force components 



F = 

mpij 



F = 

mpik 



F_ 



mp 



F_ 



mp 



cos^ y . = -k mp \tt u \3cj3c k cosaj , (24 ) 
cosa k = -k mp \tt u \SxjSx k cosa k (25 ) 
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Substituting (24) and (25) into (22) and reducing Sx -8x k , we'ss find 






7T ;; 



7T U 



COS CCj, 

COSCXr 



(26) 



Substituting cosines (23) into (26), represented via velocity 
components 

_i_ 



x iK T ) =-k mp \x ii \-v j (v j +v k ) \ 



(27) 



_l 

, 2 \~2 



^it(T) = -K P Wu | • v k (Vj +v k ) 2 ,i*j,i* k, 

Thus, in formulae (21) tangential stress ^h(t)» ^ik(T) are defined in 

the form (27), and ^ij( H ) ''^ik(H) are determined by 

nonsymmetric formulae of Newton friction law in Chapter 1, which 
in Cartesian coordinates are equal to: 



7T 2 2 — Kyyi 



^32(h) — M 



dv 



7Z\2(h) ~ M 



' ^33 — ^"zz ' 

dv 



dx 



ft 21(h) ~ t^~ 



du 



dy 
dw du dw 

ft 23(h) ~/^^~' ft31(H)~M^T> ftl3(H)~M^^, 

dy dz ox 

where normal strains are preserved in the initial form 

r A -w -n du 

™ xx =-\-P + (-M~M )divv] + ju — , 
3 ox 



K yy =-[p + (-M- M' )divv] + jLl 



dv 

dy ' 

A .x j- -n 3™ 

n zz =-ip + {-ju-ju )divv] + ju—- 

3 dz 



(28) 



(29) 
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For pretty small values of molecular viscosity in (29) when 
Raynolds numbers are great, modules of normal strains may be 
assumed equal 

ii 1 dv. 

Ffi ^ =1 -[/? + (- //-//)<^vv] + //—^l= (30) 
3 ox. 



\ ..,.-. dv. . 

= p + (-ju-ju )divv-ju — -,i =1,2,3 

3 ck. 

For simulating turbulent flows, tangential stress 

^ 'yx ' ^ xy ' 71 ' zx ' T^xz ' ^"yz ' ^"zy m equations of continuum 

dynamics 

.du du du du^ ^ dn xx dn dn 7X 

p( — + u — + v — + w — ) = pF x + — — + — — + — =- , 

dt dx dy dz ' dx dy dz 

,dv dv dv dv N „ dn dn dn 
p( — + u — + v — + w — ) = pF + - + — — + — -, (31) 

dt dx dy dz dx dy dz 

,dw dw dw 3vv\ „ dn I7 dn dn„ 
p( — + u — + v — + w — ) = pF z + — — + — — + — ^ 

dt dx dy dz dx dy dz 

are determined by formulae (21) with account to expressions 

(27),(28),(29),(30) 



7Z\. 



u(u + w ) 2 , 



du , 

n = u k 

yx dy 

dV 1 I 1/2 2x"9 

7r * y=M dx~ k m P F xx \-v(v +w ) 2 , 

d U 7 I I / 2 , 2x~2 

dz 



mp \ ZZ\ 

dw 

~dx~~ K ""' 1 



^ xz =M : ir L -k mp \?rJ-w(v 2 +w 2 ) 2 , 
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Q\? II 2 2 

*"zy =M—-K P WJ\-V(U +V ) 2 , (32) 

<5z 

7T = Li K 

dv mp 



71 



yy 



w(u + w ) 



^ 2 



Substituting expressions (30) into (31), we obtain sought for 

equations for simulating both, laminar (k=0), as well as turbulent 

flows (k *0): 

du 9 , du 1 d . du 

p— = pF x + — {ju- — [/? + (- ju-ju)diw] } + — {{* — - 

at ox ox 3 dy dy 



7 r A in J- - ^1/2 2^1, 5 . 9w 

-k[p + {-ju-ju)divv-ju — ]w(w +vv ) 2 }H {// 

3 dy dz dz 



1 



ftv, 



k[p + (—jU-ju')divv-jU :i:i -~\u(u 2 +v 2 ) 2 }, 
3 9z 



dv 9 , 9v 1 ^ 9m 

p— = pF y +—{ju- — k[p + (-ju-jU )divv -//—]• 

rf? 9x 9x 3 9x 



9 r 9v 



1 



v(v z +w z ) 2 } + — {//- — [/> + (-//-//')<ffvv]} + 

9y 9y 3 



3, * ,. .1 ,. ,. _ 9w> . 2 2n _ 9i 

+ — {ju — -k[p + (-jU-ju)divv-jU— ]v(u +v ) 2 }, 

9z 9z 3 dz 



dw 8 , 9w , 1 inj- - d« n 

p— = pF z + — {ju- — fc [/> + (-ju-jU )divv -ju — ]• 

a? 9x 9x 3 9x 



9 , dw 



1 



9v, 



xw(v +vv ) 2 }h {// k[p + (—jU-ju')divv-jU — ] 



dy dy 



dy 



2 , 2 ,~ , - S r 9W r ,1 



w(u +w ) 2 } + — {// — -[p + (-ju-ju')divv]}; 

dz dz 3 
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d 8 8 8 8 ,„„, 

— = — + u — + v — + w — (33) 

dt 8t 8x 8y 8z 

continuity equation joins these equations as well 

-^- + div^JVj-0 (34) 

dt 

and energy balance equation (49) from Chapter 1 

p — (E + \v\ 2 /2)= (35) 

dt 



= p(F,v) + — (x x ,v) + — (7T y ,v) + — (7r z ,v)-divq + pQ, 



d ,_ _ d _ _. 9 ^ ^ 
— (^,v) + — (^ v) +— i 
ox qy oz 

in which strains may be taken in the form of (29),(32) or transform 
(35) into the below form 

,dT DT DT d7\ _ 6v N 
pc v (— + u — + v— + w— ) = Or,,— ) + 
dt ox oy oz ox 

dv dv 

+ (7i , — ) + (7T Z , —) + div(AgradT) + pQ (36) 
oy oz 

In incompressible fluids p = const, if viscosity is constant, 
jU = const . For many laminar flows friction coefficient k = . 

In turbulent flows friction is caused because of mutual penetration 
of particles from neighbouring layers, which is eventually connected 
with velocity fluctuations in Euler approach. Therefore, in the 
system (33) coefficient "k" must be connected with the degree of 
turbulence in the given point of flow and, most probably, was 
selected in the form 

i 

k=k Q (u' 2 +v' 2 +w' 2 ) 2 1 w Macm , k Q = const > 0, (37) 
where u' =u —u, V = v — v,w' = w — w -fluctuations, u,v,w- 
average in the given point values of velocity components. First 
alternative - as u,v,w may be used averaged for the previous time 
interval values. Second option - use of net averaging at the given 
moment of time in a definite net node with indexes i,j,m, for 
instance, arithmetical mean value of function 
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rn _ , rn rn rn rn rn rn if" \IH (T.S.) 

J ijm w i—ijm J i+ijm J ijm J ij-lm J ij+lm J ijn—l J ijn+l' ' ^ ' 

where n - number of time layer t n (Of course, for averaging (38) it is 
possible to involve great number of surrounding nodes). Then 
fluctuations at the given moment of time t n are easily counted 

tn n — n tn n — n tn n — n /'3Q\ 

Imn Imn Imn ' Imn Imn Imn ' Imn Imn Imn ' ^ / 

where u" mn , v" mn , w" mn - actual values of velocity components, 
counted by difference scheme on computer. Third option - is to 
assume constancy k = const and varying it, to study obtained 
numerical results. Fourth alternative - is to determine values 
k = const , comparing results of numerical computations of flows 
for longitudinal flow around disc or in pipe with experimental data. 

Among these three alternatives, the most hardly executable is the 
first one, where fluctuations u',v',w' must be calculated from 
appropriate system of equations, methodology of obtaining of which 
is developed in §3. Obviously in second and third options, 
involvement of fluctuations equations is not required. 

While setting initial and boundary conditions for computation of 
turbulent flows, random disturbances must be included 121. 

Numerical solution of initially-boundary value problem for 
equations (33) and other is performed on the basis of difference 
schemes /4/. 

Note. In those points of flow (in difference methods - in net nodes), 
where velocity components are equal to 0: u=0,v=0,w=0 in (32) and 
in equations (33) cosines 

I l^ l^ 

u(u +w ) 2 ,v(v +w ) 2 ,u(u +v ) 2 , 

II 1 
w(v + w ) 2 , v(u + v ) 2 , w(u + w ) 2 , 


will become a equivocation — . In this case, basing on chaotic 

nature of turbulent flows, these cosines may be taken with the help of 
sensor (built-in procedure) of random numbers as equal to random 

106 



value from interval [-1,1], if such necessity is confirmed by 

L'Hospital rule. 

Note. Due to arbitrariness of selection of the second viscosity 

coefficient, if to include ju'= —ju , the given equations will reduce 

by not less than 40 derivatives. 
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Chapter 4. FALLACY OF APPLYING "DARCY law" IN 
FILTERING THEORY 

Polubarinova-Kochina P.Ya. in III certify that engineer Darcy in 
1856 published in 121 formula of linear dependence of filtration 

velocity vHvl from piezometric gradient s: 

dh 
v = -Z—, (1) 

as 

P 
referred to in III the Darcy law. In (1) h = \- z is declared a 

PS 

piezometric head, axis z is directed straight up, direction s makes 
angle of piezometric gradient with axis z, Z ~ filtration factor (for 
sand x ls measured in interval 0.01- 0.0001 m/sec). Important is that 
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in (1) v = Jv 2 x + v 2 + v 2 =1 v I , since due to porosity of medium, 
filtration is essential spatial motion, occurring with velocity 
v =vj +v j +vk , due to which, most probably, Darcy 
measured in his experiments not only one velocity component, but in 

/ 2 2 2 

particular full speed value v = Jv x +v +V 7 , therefore, formula (1) 

must be represented in details as follows: 

dh 



4 



v 2 x +v 2 y +v 2 =- X '— (2) 

as 



; 



In the event when direction of gradient is vertical, i.e. s=z, from (2) 
results, as a consequence, the following formula: 

v 2 x +v 2 y +v 2 z =-Z^r(— + z) (3) 

dz pg 

Again porosity of medium, fluid particles going around solid 

components of porous medium, makes the flow three-dimensional: 

V2 2 2 

v x + v + v z , which is accounted for in (3). 

In the same book III values of filtration factor % for various media 
are given. It can be expected that values of filtration factor jp will be 
different, if following second Newton law, according to wchi force 
causes acceleration, i.e. instead of Darcy formula, to use the main 
law of dynamics 

dv s _dh 

-r = -z—> v, =(v„s), is 1=1, 

at as 

and thus avoid contradictions with laws of physics. In principle it is 
necessary to base on the fact that in building mathematical models of 
this or that phenomenon, where forces causing motion act, to 
indisputably follow Newton laws, laws of conservation of matter 
and energy. Models, violating these laws must be rejected, as 
having no physical sense. 

§1. How "Darcy law" was created for multidimensional filtration 

Darcy formula (2), transformed to formula (3) of vertical 
gradient, further was reduced to formula 
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d , P x 
w = -%—( — + z), 

dz pg 
introduced only for one vertical component of velocity w, 

generalized for spatial filtering with velocity v = ui + vj + wk and 
represented in vector form as "Darcy law" 131: 

v=-xgrad(— + z) (4) 

P§ 

It is necessary to note here that the basis for (1) Darcy laid a number 

of experimental data, therefore use of the above-presented formula in 

engineering practice for rough estimates of vertical fluid flow rate in 

soil is quite acceptable, whereas generalization (4) conflicts laws of 

physics and encounters a number of problems of setting boundary 

layers adequate to the process of filtration. 

Later on an effort was made to give justification to "Darcy law" 

(4), using vicous media dynamics equation 

dv , - 

p — = -gradp + pF + pAv , (5) 

dt 

in which acceleration is equaled to zero 

dv 

— = 0, (6) 

dt 

and, taking into account essential resistance of medium to fluid 

motion (filtration), was made an assumption, that in soil consisting 

of solid particles it is possible to use formula of friction in solid 

bodies slipping (see Polubarinova-Kochina P.Ya. Ill) 

K P = *v , 

according to which, forces of viscous friction are assumed to be 
proportional to velocity 

ju/SS? = — kv 
In a result, "Darcy law" for multidimensional filtration is applied in 
the form of a system 

kv = —gradp + pF , (7) 

divv = O , (8) 
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F - density of bulk forces. 

Equation (5) is the analogue for the second Newton law, falling on 
a unit of volume 

.dv dv dv dv. , - 

p{ \-u hv \-w — ) = —gradp + pt + /jAv , 

dt dx dy dz 

according to which in model (7), (8) it ocurrs that during filtration of 
fluid particles acceleration 

dv dv dv dv dv 

— = — + u — + v — + w — = (8) 

dt dt dx dy dz 

are equal to zero, accordingly, their velocities V = const are 
constants and according to the first Newton law fluid particles must 
move along rectilinear trajectories. 

§2. Problem of setting boundary conditions 

Complexity of setting boundary conditions adequate to simulated 
physical process of filtration is already seen from one-dimensional 
equations of "Darcy law". Indeed, system of equations (7)-(8) 

kv = —gradp + pF , divv = O 

in one-dimensional flow in parallel to axis z with velocity 

v =0i + Oj + wk in projection is reduced to equations 

dp dw 

kw=-—-pg, — = 

dz dz 

From the last equation results constancy of velocity w = const . 
For defining this const, it is necessary and enough to set velocity in 
any one point w(z ) = w () = const . Then, for computation of 
pressure from the first equation 

dp 

kw Q =-—-pg 
dz 

it is also enough to set its value in one point and solution will 

have the following form 

p(z) = p(z Q ) - (kw + pg)(z - z ) 

This circumstance must be taken into account in the theory of 
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multidimensional filtration while formulating boundary 
conditions for components of velocity and pressure. 

§3. On setting boundary conditions for unknown functions u, v, 

w,p 

Here arises a number of problems while setting boundary 
conditions for equations in partial derivatives of "Darcy law", related 
to fact that in (7), (8) included are only first derivatives from 
unknown functions u, v, w, p. Indeed, in Cartesian system these 
equations have the following form 

ku = -^- + pF x ,kv = -% + pF y ,kw=-^- + pF z (9) 

ox dy oz 

du Dv dw 

1 1 = O , (10) 

Sx dy Sz 

accordingly, for u and p there must be set one boundary condition 

on lines parallel to coordinate axis x, for v and p - one boundary 

condition on lines parallel to coordinate axis y, for w and p - one 

boundary condition on lines parallel to coordinate axis z , i-e. only in 

separate sites of boundary of flow area under study, and not on the 

whole boundary as a whole. 

According to theorem of Ostrogradsky-Gauss 

iUdivvdr = i\(v,n)d<j from continuity equation (8) results a 

condition for boundary values of velocity vector components 

\(v,n)dcr = (11) 



# 



§4. Boundary conditions otDirichlet and von Neumann 

Problems of setting boundary conditions arises for two- 
dimensional equations 

ku = -^- + pF x ,kw = -^- + pF z (12) 

ox oz 

du dw _ 

ox oz 
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For simplicity of reasoning, we'll select boundary er of plain 
domain r in the form of a rectangle. According to above-said, if on 

section of boundary <j l set was boundary condition u \ CT = cp x , and 



on section of boundary c 2 set was W 



(Pi , then due to 



dp dp 
contained in equations (12) first derivatives of pressure — , — there 

dx dz 

must be set boundary conditions on opposite sections of boundary 

<7 3 , <7 A in the form p= X and p=q. 

Used until now methodology of solving equations (7), (8) uses 
elliptic equation 

div(-(gmdp - pF)) = , (14) 

k 

which is the result of substituting (7) into (8). Integral for area r 

from both parts of (14) leads in the general case to equality 

J|j<iz'v(-(graafp - pF))dT = <&(-(gradp - pF),n)d(7 = (14') 

Thus, for two-dimensional equation (14) due to association (12) 
boundary conditions will be 



cr=^P 



= (J 



(15) 



Having solved boundary value problem (14), (15) in relation to 
pressure/?, velocity components are calculated from formulae (12): 

u = (-^ + pF x )/k,w = (-^ + pF z )/k (16) 
ox oz 

throughout the area T , including sections cr 3 , (7 A of the boundary. 
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Fig.l 



The question arises: do values of velocity components obtained in 
the result of boundary value problem (14),(15) solution, satisfy 
condition (11), which for plain domain r will take a form 

(v ,n^)dcr = O 



* 



(17) 



in this case <T - curve limiting T . 



§5. Boundary conditions of von Neumann type 

On the same two-dimensional problem, let's consider setting 

boundary conditions of von Neumann on sections <7 3 and <r 4 of the 
boundary, on cr, and a 2 they are already set: 
dp dp 

( ~ac Fx)c7 ' =k(P l'(~Q +pF J°2 =k<P2 > (18 ' } 



ftp s\ , S P^\ 



(18) 



According to Darcy law (7), boundary conditions (18) result in the 
equality of velocity boundaries on sections <7 3 and c 4 



M 



"a 



= {-(p,+pF x )lk, w a4 = (-<p 4 +pF z )/k (19) 
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Thus, boundary values 

(p x ,(p 2 ,(-(Pi +pF x )/k,(-<p 4 +pF z )/k 

must be satisfying condition (17), and in three-dimensional 
problem, condition (11). There must be satisfied solvability 
condition (14') with account to boundary conditions (18), (18'). 

Paradox of application of "Darcy law" in form of 
multidimensional problem (7), (8) is in the fact that due to continuity 

equation (13) for u one boundary condition U a — (p x will be 

enough, and because of (19) there arises the second boundary 

condition U a = (— #> 3 + pF x ) I k , which is needless, since 

equation (13) of the first order by x, similar situation ocurrs also for 
w: two boundary conditions 



W 



= <P: 



2 ' 



w 



a =(-<p 4+ pF z )/k, 



dw 



because in continuity equation (13) there is first derivative — , and 

dz 
the same situation of excess boundary conditions for velocity 
components arises in three-domensional problems. 
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It is clear that boundary condition of von Neumann 

d P\ < A -\\ 

— \a=(Sfadp,n)\ a =i// 

on ' 

must satisfy according to (14') on area boundary 

cpf(— {gradp - pF),n)da = <n(— ty — {pF, n))d<j = 
JJ k JJ k k 

solvability criterion cpf (— y/)d<7 = <+] (— (pF, n))da 



§6. Paradoxes of "Darcy law" (7), (8) 
First paradox: "Darcy law" (7), (8) in conservative field of force 
describeees vortex-free, i.e. potential flow. 

Let's compute circulation of velocity vector cl^Jr by arbitrary 

r 
shape r, which in filtering theory by "Darcy law" (7), (8) is equal to: 

c[<, dr j= j((-gradp + pF) /k,dr) (20) 

r r 

Let the flow ake place under action of conservative gravity 

F = -gradn , (21) 

where 77 - potential energy. 

In this case, integrand in (20) is equal to 

{{-gradp + pF) I k, dr) = {{-gradp I p - gradll)plk, dr) = 

= -{plk)d{pl p + n) 

Due to this for p = const , k=const integral in the right part (20) 
will be equal to zero 

ji,dryj{{-gradp + pF)/k,dr) =-{plk)jd(pl p + 77) = (23) 

r r r 

As the result from (20) for potential forces obtained is equality to 
zero for circulation of velocity vector within any closed shape 

cj"<, <ir 2=0 (24) 

r 
According to Stokes theorem, there is transition to surface integral 
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jH,drj= ^(rotv,n)ds=0 (25) 

r s 

Here S - arbitrary surface, stretched to shape 7". By Du Bois - 
Raymond lemma, from (25) results equality to zero of velocity vortex 
(vorticity) 

rotv = (26) 

throughout the area of flow. 

Another proof of equality (26) results from identity 

rotv = rotgrad (p + pll) = 

k 

Flows, in which velocity vortex is everywhere equal to zero 
rotv = , are, as is known, potential floes of ideal fluids 13/. 
Velocity potential for two-dimensional flows was considered by 
Polubarinova - Kochina P.Ya. in III. In the same book III it is 
admitted that resistance forces, which affect fluid particle in a pore, 
depend on internal friction of fluid. Internal friction - is the 
integral feature of viscous fluids, in viscous fluids potential flows 
are impossible, vorticity is different from zero rotv ^ . Thus, 
filtration flow by "Darcy law" is initially irrotational flow of ideal 
fluid, in ideal fluids there are no friction force, no 
moleculartransport, which conflicts reality, in particular when extra- 
heavy oil motion is studied. 

Second paradox: "Darcy law" (7), (8) contradicts to law of 
energy conservation. 

Velocity potential is introduced by formula 

v = gradO (27) 

Substituting (21) and (27) into equation (7) of Darcy law, we'll find 
gradikO + p + pTT) = , (28) 

which results in integral 

k0 + p + pIJ = const , (29) 

where const is the same for the whole flow area. On the other part, 
for steady-state irrotational flows of ideal fluids there is the integral 
of Lagrange -Cauchy 131. 
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|-|2 
V 

p J — L - + p + pTJ = const (30) 

Obviously, (29) does not result in integral of Lagrange-Cauchy. 
Equating const in (29) and (30), we arrive to paradoxical equality for 
"Darcy law " 

\grad@\ 
k& = p l - L_ ( 31) 

Paradoxicality of (31) is connected with the fact that velocity 
potential & is determined accurate to arbitrary constant 121, and 
gradient neglects this constant. 

It is easy to see that integral of Lagrange-Cauchy (30) expresses 
law of conservation of mechanical energy in the unit of continuum 
volume. 

And in flows by "Darcy law" appears integral (29). Obviously, 
integral (29) of "Darcy law" conflicts with law of energy 
conservation (31). 

Third paradox: "Darcy law" (7),(8) conflicts with law of variation 
of momentum or second Newton law. 

Major dynamics law = / i f k states that any force V^ f k 

dt t V 

causes variation of momentum mv , therefore is formulated as law 

of momentum variation. Having differentiated, it is possible to 

present it in the form 

dv _ dm 
— + v — 
dt dt k 

Any motion of material medium in Newton mechanics is 
subordinated to this law, accordingly, motion in porous medium 
arranged in the form of "Darcy law" 

kv = —gradp + pF 

must follow this law. Simple comparison of these two laws, one of 
which is pseudo "Darcy law" leads to the following fact 



m — + v = ^A 
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y^ \f k = — gradp + pF , v = kv , m — =0 

^ dt dt 

From the last equality it results that in pseudo- "Darcy law" 

dv 
acceleration is equal to zero — = . 

dt 

But his means that according to 1 st Newton law, particles must 
accomplish straight line motion with constant velocities, or rest, 
which is possible only in one-dimensional flow towards axis "z". 

In hydrodynamics, descriptive illustration of such flows are: 
Couette flow, channel Poiseuille flows between parallel planes, 
Hagen-Poiseuille flow in circular pipe etc. 131. 

Due to this application of "Darcy law" (7), (8) for computation of 
two-dimensional and three-dimensional flows at once loses sense 
because of contradiction with 1 st Newton law, as soon as curved 
streamlines appear in the flow. There can be no rectilinearity of 
streamlines due to above-given paradoxes with boundary conditions. 

^ dv dv dv dv _ . 
Equating to zero of acceleration Yu hv Yw — =0 is 

dt dx dy dz 

justified by the following considerations. First: in steady-state flow, 

dv n 
of course, partial time derivative — = , but for nonsteady flow 

dt 

dv 
local acceleration is not equal to zero — ^ , and undertaken in 

dt 

monograph III effort to equate it to zero ca not stand stricture. 

dv dv dv 
Second: spatial transposition members u hv Yu — are 

dx dy dz 

negligibly small. Steady-state of filtration, does not probably always 
take place due to action of external conditions, and products in 

dv dv dv 

u hv Yw — may be not small at all. Let's show this on 

dx dy dz 

analogue of only one product from thus expression of transition, for 

example, on w — . Let w = sF(cosaz), where a = const »1 is 
dz 
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very large but finite number, s - parameter, providing 

infinitesimality of unction w, F- limited with derivative differentiable 

function of its argument, then derivative will be equal to 

dw dw 2 

— = -sF'asmaz . Therefore, product w — = -as FF'smaz 

dz ' dz 

is the finite number, it is enough to include a = const * £~ . 

Sufficient number of similar examples may be given, for example, 

functions with weak discontinuity, in order to prove that equating to 

zero of acceleration 

dv dv dv dv _ 

\-u hv \-w — = 

dt dx dy dz 

is unjustified. Unjustified, first of all, because of contradiction 
with the major dynamics law - in particular, with 2 nd Newton 
law, which says: action of force results in acceleration, not equal 
to zero: 

dt 
(Paradox - is it real that fountains of oil and water from underground 
occur under negligibly small accelerations?) 

As was noted in (14), widely spread method of solving equations 
(7), (8) of "Darcy law" is reduction of this system to one elliptical 
equation 

div(— gradp) = div{— pF) 
k " k 

in above-indicated boundary conditions. 

It is necessary to admit that the most important criterion of 

accuracy of obtained solution of equation (14) is that for velocity 

vector v , defined by formula (7) "Darcy law", not only condition 

(1 1) must be exercised, but first of all equality to zero of acceleration 

dv dv dv dv n 
\-u hv \-w — = , 

dt dx dy dz 

i.e. pressure in steady-state filtration must satisfy equation resulted 
from the below equality: 
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(- J + pF x ) ^ ((-gnafc + pF)/k) + 
ox ox 



dy y dy 



+ (-^ + pF y ) — ({-gradp + pF) Ik) + 



+ (-^ + pF 7 )^-((-gradp+pF)/k) = 
oz ' oz 

in all points of investigated area, which is not an automatic 
consequence of solution of equations (7), (8), as this occurs in flows 
of Couette, Poiseuille between parallel planes, Hagen-Poiseuille in 
circular pipe, where acceleration is identically zero 131. 
Obtained is 4 th paradox: pressure in "Darcy law"(7),(8) is defined 
through solution of equation (14), on the one part, according to the 
first and second Newton laws, pressure must be also solution of the 
system of 3 equations: 

ox ox ox oy oy ox 

oz oz ox 

ox ox oy oy oy oy 

+ (pF.-*)|-[(pF,-|)/*]=a 

oz oz oy 

(pf, 4>U(pf, -%ww, -fjiy, 4)/« + 

<3x <3x 3z oy oy oz 

oz oz oz 

Alltogether, for one pressure function p arise 4 equations which once 
again shows irrationality of applying "Darcy law" (7),(8) in the 
theory of multidimensional filtration. 

Due to the above-mentioned reasons, application of so-called 
"Darcy law" in the form of system of equations (7),(8) in the 
filtration theory is erroneous. 
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§7. In «Darcy law» only in one-dimensional steady-state flow, all 

three laws of physics are observed, if equation of continuity 

accounts for mass source 

Let's consider one-dimensional steady-state flow of ideal fluid 

with velocity components u = 0, v = 0, w ^ 0, parallel to direction of 

gravity g : 

dw dp ,„_ 

pw— = --f—pg (32) 

dz dz 

This equation is projection to axis "z" ofEuler equation 

p — = —gradp + pF 
dt 

Multiplying both parts of this equation by elementary volume 

St , we'll arrive to the second Newton law 

dm — = / , 
dt 

where / = {—gradp + pF)Sr - resultant force, dm — pSr - mass. 

Further on, taking into consideration filtration in soil, in the law of 
mass conservation let's introduce strength of source or drain in the 
following form 

Pi 1 1 ) = £ (33) 

Qx 3y dz 

This equation in case of one-dimensional flow u = 0, v = 0, w ^ 0, is 

reduced to the form 

dw 
p— = k (34) 

dz 
Projection of acceleration 

dv dv dv dv 

YU hV YW 

dt dx dy dz 
for steady-state flow towards z will be equal to 

dw k 

w — = w — 

dz p 
Substituting (34) into (32) 
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dp 

kw = pg 

d z (35) 

we'll obtain formula (3). In conservative field of gravity force 
potential energy is equal to 77 = gz + const , due to which law of 
energy conservation (30) (integral of Lagrange-Cauchy) for one- 
dimensional flow takes the form: 

2 
W 

p \- p + pgz = const , (36) 

having differentiated which on z we obtain dynamics equation (32). 
On the other part equation (32) is brought to equivalent form 

—(^— + P + pgz) = 

dz 2 
from which obviously results law of conservation of mechanical 
energy (36) in unit of volume of medium. 

Clear is the following simple circumstance - according to formula 
(1) of engineer Darcy, first problem of dynamics is solved: knowing 

velocity of fluid w(z Q ) and pressure p Q = p(z Q ) on the given level, 
it is possible to compute force affecting surface at the distance of 
(z — Z ) from the selected level. Let's demonstrate it on one- 
dimensional flow. From (35) it is easy to derive 

Mz ) = -y(^^ + pg), 

k z-z 
wherefrom pressure on level z is computed: 

P = Po -(z-Z Q )[kw(z Q ) + pg] , 

where in the right part, there are measured quantities. Preliminary, 
for calculation of coefficient k formula (34) is applied. 

Integrating (34) by height of filtration, it is easy to establish 

p(w(z 2 )-w(z l )) = k(z 2 -z l ) 
wherefrom the coefficient is calculated 

p(w(z 2 ) - w(zj) 
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where, obviously, levels z 1 and z 2 are defined in experimental 

soil sample to be so, that on these levels velocities w(z 2 ), w(z x ) 

are known. 

Thus, founders of so-called "Darcy law" instead of the first 
problem of dynamics are trying to simultaneously solve the first and 
second problems of dynamics, which in the result led to above- 
mentioned paradoxes. 

§8. Forzheimer model of filtration in saturated porous media 

As opposed to "Darcy law" more suitable (from the point of view 
of observance of laws of physics) for describing filtration in 
saturated porous media is the Forzheimer model 151 



Sv ,._,, s v ., 



i 



p f — = -(p(V(p + p f gz) + p f —v+p f c f K 2 I v I v) , (37) 
ot K 

divv = O (38) 

These equations take into account local acceleration, equation (37) 
is represented as major dynamics equation, thus, observed is second 

Newton law, as opposed to "Darcy law". In (37) *P - coefficient of 

medium porosity, K - permeability coefficient, f - fluid density, 
c f - dimensionless Forzheimer friction coefficient. Accounted is 

also fluid viscosity, since kinematic viscosity coefficient V is 
included 

Regretfully, for the system of Forzheimer arises similar to the 
above-mentioned problem of setting boundary conditions for 
pressure and velocity vector, and adequacy of obtained solution to 
simulated physical process depends on boundary conditions. 

§9. Forzheimer model correction 

The indicated problem is easily removed, if Forzheimer model is 
altered, making it close to equations of viscous fluid, at the same 
time replacing local acceleration for an individual one 



123 



I- 1 

P/ — = -^[V(77 + p / gz) + p / — v+p f c f K 2 \v\v] + juAv, 
at K 

divv = O (39) 

Numerical solution of initially-boundary value problem for system 
(39) is carried out by methods developed for viscous fluids equations 
161. 

§10. Numerov model of filtration in saturated porous media 

Numerov S.N. in 1968 turned attention to necessity of taking 
account of inertia forces in major equations of filtering theory and 
proposed to use in /8/ the following system of equations: 

1 ^T* 1 

+ -(v,V)v + gradh + f (v)v = 0, divv =0 

go dt go- 
Thus, Numerov S.N. for the first time drew attention of researchers to 
inadmissibility of violating Newton laws of dynamics, which state 
that: "mass by acceleration is equal to force", if acceleration is equal 
to zero, a body performs straight line motion with constant velocity 
v = const or is resting. 

§11. Interest to "Darcy law" is related to applications of complex 

variable theory 

Indeed, in plane problems (7)-(8) elliptic equation is obtained for 

pressure, if pF = const : 

d 2 p d 2 p _ 

dx 2 dy 2 

for finding solution of which effectively used was theory of complex 
variable functions, which is in sufficient details described in /l/,/4/ 
etc. At present, development of computing machinery and 
mathematics allow obtaining of numerical solution of equation (40) 
with the given boundary conditions in a matter of seconds, therefore 
equation (40) is of interest more as learning material 111, than 
scientific and research, due to above-mentioned reasons. 
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Chapter 5. FALLACY OF THE "FICTITIOUS DOMAIN 
METHOD" 

Recently, there appeared many works devoted to fictitious domain 
method (f.d.m.), idea of which was in due time proposed by V.K. 

Saulyev. 

Fallacy of "f.d.m." is that solutions of two independent initially- 
boundary value problems for differential equations in partial 
derivatives are connected in arbitrary manner. There are especially 
many works on "f.d.m." for Navier-Stokes equations. 

Idea of "f.d.m." is as follows. Let in "inconvenient" for 
introducing computational meshes of nonrectangular physical region 
D with boundary S, initially-boundary value problem for Navier- 
Stokes equations is solved 
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dv - - 

f{— + (v,V)v] + Vp = ^v + pF,(V,v) = 0,v\ t _ =d,v\ s =(p(i) 
at 

In "f.d.m." to regions D added is fictitious domain D £ with 
boundary *J e , in which initially -boundary value problem is solved, 
for example, for case tp e — vJ like this: 

p[^+(v;,v)v;]+v Pe =Mv;+^-^,(v,v;)=o,v;u=rf>;i s =^(2) 

at s 

It is supposed that region D U D E is convenient for introducing 
net domain while solving with the help of finite-difference methods. 

Here, it is necessary to pay special attention to additive term to 

V 

equation (2) — - , in denominator of which, artificially introduced 

s 

parameter £ tents to zero: £ —> ! 

It would seem that a single look at equation (2) is enough to reject 
f.d.m. Obviously, £ - equation (2) absolutely does not coinside with 

Navier-Stokes equation (1) with small parameter S — >■ , but tents 
to it under conversely large parameter value S — > 00 . And the 
following questions arise at once. 

Question 1. What must be specific numerical value of dimensional 
parameter £ (dimensionality [ £ ]=kg/( M c) )? What are criteria for 
its selection? 

It cannot be equal to zero. Assuming that £ =0, fictitious problem 
(2) transforms into a nugatory equality 



v =0 



(3) 



throughout the whole fictitious domain D , which means there will 
be no necessity in initially -boundary £ - problem (2). At the same 
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time it is obvious that continuation of solution of initial real problem 
(1) into fictitious domain D will be non-zero 



K 



V 



D *° , (4) 



vector v 

s 



i.e. there is a contradiction with (3). 

Question 2. If S ^ , what boundary conditions must meet velocity 

s = cp s on boundary S £ , i.e. how must be selected vector 

q> s so that solution of fictitious S -problem would correspond to 
continuation of solution of initial real problem (1)? 

If there is no precise match of solutions of problems (1) and (2) in 
the domain D £ , use of "f.d.m." in finite-difference methods with 

introduction of common net domain loses its sense. In addition to 
these questions, the below-given contrary instances will be enough 
for understanding irrationality of "fictitious domains method". 

Contrary instance 1. Let's consider problem of "f.d.m." from the 
point of view of the second Newton law. For this, we will represent 
equation of vicious fluid dynamics (1) with the help of substantial 
time derivative 

p^ = -Vp + juAv+pF, (5) 

dt 

multiplying which by lementary volume St we arrive to formulation 
of the second Newton law 

Sm^- = f, (6) 

dt 

where dm = pSr — mass, acceleration, 

dt 

f = (— Vp + juAv + pF)8r - resultant (major) force. 
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Multiplying dynamics equation in system (2) by St , by analogy 
with expression (6) we obtain second Newton law in the form 

Sm^ = f-^Sr , (7) 

dt s 

v 
where the last member — - St determines reactive force, which tents 

s 

to infinitely large values since f — > . It occurs that by virtue of (7), 

in fictitious domain the particles of fluid are affected by infinitely 

large proportional to velocity forces dependent on parameter S . 

Contrary instance 2. Let's start from a simple ordinary differential 
first-order equation 

y'=y,y(0) = l, (8) 

Let's condition that solution of Cauchy problem (1) is ought in real 
domain < x < <x> . Exact unique solution of the problem is obvious: 

y(x) = e x (9) 

Let domain — co < x < according to idea of "f.d.m." be fictitious 
domain. By analogy with (2) this problem (8) is connected to £- 
problem 

y £ =j £ --j £ ^ £ (o)=i do) 

£ 

Solution of Cauchy problem (10) has the following form 

(l--)x 

y e = e £ (11) 

Let's note, that initial condition j^(0) = l in (10) is satisfied by 

solution (11) with any £ ^ . When S > U solution (11) in point 
x=0 will have uncertainty in the index 

o 
y s (Q) = e~ T *° (IT) 

128 



Thus, in "f.d.m." appears a problem of setting boundary conditions 
when £ — ^ . 

Theorem 1. When £ — ^ solution (11) of Cauchy fictitious 
problem (10) may by arbitrary large value differ from extension 
into fictitious domain of the real solution (9) of the initial problem 
(8). 

Proving. Let's assume x = — b e (— co,0] ,b ^ . For this arbitrary 
point from fictitious domain, solution (11) of "fictitious" problem 
will be equal to 

y e (-b 2 ) = e * X =e x e £ =e~ b2 e £ =e~ 6 V =y(-b 2 )e £ , (12) 

where y(— b ) - extension of real solution (9) in fictitious domain 
(—oo,0] . Let's make a relation 

y -^^ = e^ , (13) 

y(-b 2 ) 

from which it can be seen, that when £ ^ U solution (11) 
y £ (— Z? 2 ) is different from extention of real solution (9) y(— b ) into 

b]_ 
fictitious domain by arbitrary large value, since e s — > oo tends to 
infinity, QED. 

Let's consider their difference in the real domain. Let 
x = b 2 e[0,co),b ^0. Then, equalities 

(I--)X -- 2 "— -— 

j^(Z? 2 ) = e £ =e x e £ =e e £ = )>(£> )£ £ result that in the real 

domain as well difference between solutions (9) and (11) when 

_b]_ 
£ ^ " also becomes arbitrarily large, because e f — »0 and 
fictitious solution tents to zero y s (b ) — » , while real solution is 

nonzero y(b ) ^ 0. In this example, solution (9) of the initial Cauchy 
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problem (8) as well as solution (11) of "f.d.m." problem (10) clearly 
satisfies boundary condition y £ (0) = 1 , y(0) = 1 . 

Contrary instance 3. Let's consider another example with 
nonhomogeneous equation 

y'=y+l,y(P)=0, (14) 

Exact solution of which is the following 

y(x) = e x -\ 
(15) 

Let's condition that solution of Cauchy problem (14) should be 
sought in the real domain < x < co . Let domain — oo < x < 
according to "f.d.m." be a, fictitious domain. By analogy with (2) this 
problem (14) corresponds to £ -problem 

? e = y e --y B +hy B (.o) = o (16) 

£ 

Solution of Cauchy problem (16) has the following form 

(l--)x 1 

y =e E -1/(1 ) (17) 

£ 

Solution of fictitious problem (17) approximately satisfies initial 
condition for £ ^ 

y e (0) =1-1/(1--) (18) 

£ 

When & — ^ just as in (11'), there arises uncertainty in index 

-^- 1 
y E {Q) = e '-"-1/(1 -) (18') 

Let's assume x = — b e (— oo,0] . Solution (17) in this point will be 
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*1 1 *1 61 i 

y (-fc 2 ) = «T*V -1/(1 — ) = (e" fe2 -l)e^ +e £ -l/(l ) = 

£• fi- 

fe^ fel . 

= y(-Z? 2 )e £ +e £ -l/(l--), (19) 

£ 

here y(— Z? ) there is an extension into the fictitious domain of the 
real solution (15). From (19) is seen the infinitely large difference 
between the fictitious solution y £ (—b 2 ) and actual solution (15) 

y(— b ) , arising when £ — > . 

In contrary instances 2 and 3, it is shown that solutions "f.d.m." 

when S — > may differ by infinitely large value from extension 
into fictitious domain of the real solution because of (13) and (19), 
and these particular solutions in fictitious domain are used in mesh 
methods as additional conditions for initial problem in the real 
domain. 

Contrary instance 4. Let's show inconsistency of "f.d.m." for 
flow of viscous incompressible fluid in a channel with parallel walls 
(Poiseuille flow), that is for boundary-value problem. For the given 
laminal flow, Navier-Stokes equations will take the following form 
111: 



d u dp 
dy 2 dx 



J u -r^ = -^- ( 2 °) 



Exact solution of this equation: 

M = __L.4P(^_/) t ^- = const (0 , (21) 

2jU dx dx 

uQy) = ui—b^ = O (22) 

In fictitious domain D £ by "f.d.m." £ - equation is solved 

dp d 2 u c 1 

JU-ty- u c ( 23 ) 



dx dy s 
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Exact solution of this equation in physical domain D with 
boundary conditions (22) will be: 



i I 

(— ) 2 -y -y(— ) 2 

dp ri e £// + <? £// , 



(24) 



Difference of solutions (22) and (24) 



|«(y)-«,(y)| = -— l^-(6 2 -y 2 ) + ea- 



y(— ) 2 



->'(-) 2 



g * +e * 



Jj 2// 



«±->i 



-H^ 



")l (25) 



g * +e * 



shows that limlw(y) — uAy)\ ^ 0, i.e. for S^O 

e-->0' ' 

w^ (y) does not tend to w(}0 • 
Now, let's consider equation (23) in fictitious domain D £ . 

According to "f.d.m." (2) at the boundary o e of the fictitious 
domain D introduced is homogeneous boundary condition 

U s (—y H ) = \) , besides, on the boundary of physical domain given 
is U s y—b) — 0. Under these boundary conditions, exact solution of 
equation (23) in fictitious domain D £ have the following form: 

l 



U £ =S 



dp 
dx 



K2.,...b .y H 



t.^wl^ 



(_)2.( y+ +^£L) _(,,+ +^L).(_) 

e m 2 2 +e 2 2 * 



,yn > w l,i 



( ^4)-(-) 2 



e 2 2 w -\- e 2 2 %l 



(26) 



Similarly to (25) difference of solutions (21) and (26) in fictitious 
domain D does not tent to zero when s —> . This implyies, that 
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solution (26) in fictitious domain does not add up to solving Navier- 
Stokes equation (21), therefore, in non-physical fictitious domain, 

D e S — quation may not be applied. 

Obviously, exact solution of Navier-Stokes equations (21) on the 

boundary " — y H " e S of the fictitious domain is nonzero 

u(-y H ) = ~- d ^(b 2 -y 2 H )^0, (27) 

2ju ax 

therefore, homogeneous boundary condition U E \— J# ) = U 

"f.d.m." (2) is fully in contradiction with exact solution (21), 

according to which u(— y H ) ^ 0, which can be seen from Fig. 

So, in this contrary instance, solution (26) in fictitious domain D e by 

infinitely large value differs from extension of solution (21) into this 
domain, accordingly, the given inequality of boundary values 

u s (-y H )=Q, u(-y H )*0 (28) 

proves senselessness of applying "f.d.m." (2) in problems of 
hydrodynamics. 

This hydrodynamics contrary instance 3 convincingly shows the 
problem of setting adequate to extension into fictitious (trans- 
boundary) domain to solution (27) respective boundary conditions 
for fictitious problem (23), because homogeneous boundary 

condition U £ (—y H ) = "f.d.m." (2) turned out to be fallacious. In 

this example a simple domain was considered, at this analytical 
solution (21) of the real problem is known in advance, thanks to 
which exact bounadary condition (27) is also known. 

In solving three-dimensional problems, where analytical 
solution of type (21), is not known in advance, complexity, 
frankly speaking, insolubility of the problem of setting boundary 

conditions on fictitious boundary b E , not matching the real 

boundary S , in three-dimentional problems of hydrodynamics, 
is indisputably obvious. 
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Given above contrary instances prove that true is 

Theorem 2. When £ — ^ solution of fictitious problem in 
fictitious domain will not match with extension of solution into 
fictitious domain of a real problem. 

Main purpose of "f.d.m." is in the use of solution of fictitious 
problem in fictitious domain as additional boundary conditions in 
mesh solution techniques for initial primal problem. Contrary 
instances and the theorem prove fallacy and uselessness of "f.d.m." 
for solving hydrodynamics equations in nonstandard domains. 
Theorem 3. When s — > oo solution of fictitious problem 

/$■ + (v„ V)v e ] + Vp e = iM, + P F e - ^£-,(V,v.) = 0,v, L= dj s l,= <j> t 
at s 

in fictitious domain will not match with extension of solution into 

fictitious domain of a real problem: 

8v - 

fA— + (v,V)v] + Vp = ^Av + pF,(V,v) = 0,v\ t _ =d,v\ s =(p 
at 

Theorem proving is supported by previous results and common 
sense, which is in absolute impossibility to set such initial and 
boundary conditions 

v e \,*,= d e ,v e \ s =p e 

in a fictitious problem, under which solution of a fictitious problem 
would exactlymatch extension of the initial real problem into 
fictitious domain (see (28)). It is clear that if there is no such exact 
fit, solution of a fictitious problem, engaged as an additional 
condition, will injure solution of a real problem, which is 
unacceptable. 

Suggestion. In problems of flows in domains with solid curved 
boundary more suitable is extrapolation to mesh points covered by 
fictitious domain, values of obtained in points of real domain, spline 
functions formulae can be used for extrapolation, and pressure must 
be computed from continuity equation. It is much easier and more 
efficient than solution of a lengthy, not corresponding to laws of 
physics system of sophisticated equations with ungrounded initially - 
boundary conditions, in fictitious domain. 
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Chapter 6. CONSTRUCTIONS OF MONOTONOUS SIMILAR 

SCHEMES 

§1. Monotonous schemes constructing technology 

For building monotonous similar schemes second and more upper 
order approximation of younger terms effective is next technology 
(see Jakupov K.B. /If). 

Consider initial-boundary task for one-dimensional equation with 
convection term: 



,dco dco. 8 co 

+ u ) = M — 2 

dt dx dx 



p(—+u—) = J u— r + f,a>\ t=Q =d, (1) 



<3)(0, t) = cp x (t),a>(a, t) = <p 2 (t),0 < x < a, t e [0, T] 
Introduce net: spatial 

n h ={x i ,i = 0,...,N x ;y j J = 0X...,N y ;z k ,k = 0X...,N z } 

and by time Q r ={t n =nr,n = 0,1, . . . ,N T } with steps 

h xi = x i - x._! > 0, i = 1,2, ...,N X , h yj = y j - y ._, > 0, j = 1,2, . . . ,N y , 
h zk =z k -z k _ l >0,k=l,2,...,N z ,r n =t n -t n _ 1 >0. 
Using designations of net functions: 

fijk = / \ x i ^y jiZ- k ,t n ), ji+ijk = / \ x i±\ ^yj^ k ,t n ), 
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Jij±\k - j \ x i ' y,i±i iZ-k^nh hjk±\ - J \ x i > yj > Zfcu > f „)» /,$ - / (•*,■ » y,- , Zj , t n+l ) , 

Introduce difference derivatives: 



/■ n rn rn rn rn rn j . 

, _ Mjk ~ Jijk rn _ J ijk ~ Ji-ljk rn _ Ji+ljk ~ Ji-\jk ^ _ h x j+\ + <l x% 
Jx — , >Jx — , >Jx — . >"> 



h. 



xi+\ 



h, 



Km+K 



f - f. 

n _ J ijk Ji-ljk rn _ *■ / 
' ' J XX ~ 1 v 



/■" = 

J X 



h, 



^ rn rn rn rn rn rn 

1 /Ji+ljk J ijk J ijk Ji-ljk \ rn J ij+lk J i 

")tjy = 



ijk 



h* K+\ 



/■n r n r n r n 

ijk J ij-lk rn J ij+lk J ij- 



fj" = 



K 



/n 
v = 






K+i+K 



' n yj r, 'Jy - 



/•n rn 

ijk ~ J ij-lk 



fi„ 



1 rn rn rn rn rn rn 

/n 1 /J ij+lk J ijk J ijk J ij-lk \ rn J ijk+1 J i, 
yy = 7 ^ i i ' 'Jz = 



ijk 



n yj K* 



h. 



f" = 



yj 

ijk J ijk- 1 rn J ijk+l J ijk— I 



h 



/■ n r n 

ijk J i 



h 



/■ft rn 
ijk+1 Jijk—i ,_ 
z = '« 



afc+i 

Kk + i+Kk 



zk 



Kk +i + h zk 



'"■zk 



/•n rn ^ rn rn rn rn 

ijk J ijk— I rn ^ sJ ijk+1 J ijk J ijk J ijk — 1 

z = I ' J zz = 7 V 



h. 



h zk Kk+l 



h 



) 



Technology III concludes in using decision equation (1) by next 
way. Take place scheme 



dco_ I < I +u[ 

dt + 2 



col'-coL K d 2 ® 
2 dx 2 



J i " / >'-l. + .« 



h. 



+ 0(h 2 J 



u]-\u] 



<i-< K+\$ a + {h 2 ) 



h 



'.«'+! 



2 dx 



_ d co - r _ r . ,. 
jU—Y + f,jU=-,f 



+ (2) 



M i f 



dx 



P P 



In (2) is made replacement from (1) ^ , 

a>" -co", h. 1 \dco dco 



d co 1 \dco dco j] 

T = —\ + u f\- 

dx ju [dt dx J 



dco I u" I +w" 

— + — - 

dt 2 



• '-i _ + •«' 



_. +M--/+0(^) 

2 ju [dt dx J 



+ (3) 
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+ - 



n i n i 

U; -\U; I 



c LZ ^_^ ±L i|^ +M ^_ / j +0(/j2+i) " 

L 2 fi{dt dx J 






after that in (3) using approximate 

dm _ coT l -® n i 
dt TV, 



+ 0(r„ +1 ) = <+0(r n+1 ), 



5© I w" I +M 

W =— ' 

dx 2 



4-dU 



+ 0(hj 



+- 



n i n 

U; ~ I U; 



dU-4 



xi+\ 



+ 0(h xi+l ) 



,(4) 



M 



d co 2ju t G>M-e>" a>*-a>\ 



dx h xM +h xi 



h._ : 



i -) + 0(Mh xM -h xi )) + 0(m 2 xi ) 



Reduction similar terms prove clear scheme 2 order mistake 
approximation of convection term 



I "li n fl ^-.n n I M I ,-H ^-.n 



/' 



. p \U { \+U, , \U { l-M, , , /z ■,, 

1 + -^- (— ! ^ h xi + — ' L h xM ) * ,+1 

2// 2 " 2 x,+1 






K: ft. 



In (4) is approximation 1 order. By the same technology, use in (3) 
instead of (4) approximation 2 order 



do I m" I +u" 



I I 



dx 



col-coU ^h xi tfco 

h; 2 dx' 



+ 0(h 2 J 



+ 



n i n 

«, - I U: 

+ ~ ~ 



®M-®! Km^O +0{ j2 ) 



'.Vl'+l 



2 cT 



(6) 



further, state 



3 o> 1 ,3o> 3» 



& 2 // ' dt dx 
term, prove scheme 



= — { Yu /} in (6) and put similar 
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n n 



* „ \U) \+U l CO- -co-, u.-\u t \co i ^-co i -> 
2 ^ 2 ^ +1 



jU 



-Jl rn 

i=o 2// 2 2 

^° = d„i = 1,...^ -1,a£ =^ n ,< =^,Vn,M= 0,1,2,..., 

which by ffl"=(ffl" + -co") It are clear schemes, convergent and 
stable in network form C by execution of conditions 

Ip( — -+— — -) + 

Vi 2 ^ 2/2 .« + i 
+ ^ -(— + — ) JO,V(/,n), 

L [ V~ ( T H * + O h *M )] 

M) 2 // 2 2 

by <y" = (<y" — <y"~ )/r are implicit schemes, of course convergent 
and absolutely stable in the same form C. 

Error of scheme (7) depends on M>0. For M = scheme (7) has 
"approximated ductility" and 1 order of exactness and well-known as 
Buleev-Petriscthiv scheme. For M = 1 has 2 order of error of 
approximation of convection term III: 

Vi =0(T n+l ) + 0(\u:\(h 2 xi+ h^ M ) + 0(ju(h xi+1 -hj)+ 0(<uh 2 xi ), 

For M = 2 on analytical grid h xi = h xi+l =h xi = h x ,\/i , new scheme 
(7) takes type 
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2 A. 2 A 






w )+1 x.n/ ; -r u^_[ 



2/u 2/u 

For error of scheme z" =a>" -co" , where <y=<y(jc,£) — 
accurate decision of task (1), take place estimation 
llz"ll<rmaxll//ll,n = l,2,..., 

0<fc<n 

whence flow convergence of data of monotonous schemes (7), but 

also stability of schemes by starting condition d and by free term f, 

that is for disturbance of scheme S" —CO n —a>" take place estimation 

\\S n M\s\\+Tmax\\^ k \\,n = l,2,.... 

0<k<n 

Generalization of this technology on equations in cylindrical, 
spherical and other coordinate system, but also on multidimentional 
equations don't arouse special difficulties. 

Necessary notice, that in 131 for equation with younger derivative 
(ku')'+r(x)u'—q(x)u = —f(x) , 
monotonous scheme was built by Samarckyi A.A. technology 
X(ay-) x +b + a (+l) y x +b~ay x -dy = -<p, Z = U(l + R),R = 0,5 1 r I hlk, 
with error (h 2 ) on analytical grid. 

§2. Monotonous similar schemes 2 order of approximation of 

convection terms in equations 

with constant coefficient of molecular transport 

Above-mentioned technology III for multidimentional equations is 
applied by similar way. 
For two-dimensional equation of parabolic type 

,da> dco dco. ,d 2 co d 2 co, 
1- u h v ) = ju( — y H 2 

dt dx dy dx dy" 

similar (5) clear scheme takes bype 



p(— + u — + v—) = J u(— r + — Y ) + f,co\ t=0 =d,co 



= <P> 
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«+l n i n i , n n \ n \ \ n i , n n \ n \ 

CO;; -CO:; I U:: I +U :: «;; ~ I U, ,■ I I V,-, I +V,-; V, ~ I V,, I 

p(-* L +^ L co" + ^ '-to" +^ '-to" +^ '-a") = 

T n+l 2 A 2 x 2 ' 2 v 

= -^- (fi , : _^ ) + J^ K _^ ) + / » 

Ki + i+Ki h yM +h yJ 

where coefficients by dissipative terms are equal 

i n i . n i n i n i n \ , n i»i n 

I v" I +v" I v" I -V* 

A =[^ + p(-JL_JL^ + _L_JL^ )] / fi> 

I n I . « I n I n 

1 W„. I +«,.. I «,,. I -«„■ 

4 = [// + rt- 1 ^ K + -^-*- ^ )] / q 

On analytical grid h J+1 = h yj = h ,\/j,h xM = h xi = h x ,\/i these 
coefficient are simplified 

\v"\h \u"\h x \v"\h 

K =(//+p^-^/[i+/?Hr^+^^)], 

2 2/i 2/i 

\u"\h x \u"\h x \v"\h 

2 2// 2// 

The scheme is stable and converge in form C by execution of 
conditions 

I n i . n i«i n \ n \ . n \ n \ n 

Vi 2/l *.- 2A„. +1 2/j y . 2fc y/fI 

2 ^x , 1 l ^ 2 K' 1 1 N T> n w ,. • s 

+ — — ( + — ) + — ( + — ) _f0, V(i,j,n) 

Km+K Ki + 1 K Kj + l+ h yj Kj + l Kj 

It is not difficult to build by this methodology absolutely 
convergent and stable schemes of Krank - Nicholson type or 
schemes of method of fractional steps. Exactly in work III similar 
approximation is applied for Helmholtz equations in scheme of 
variable directions. 

For three dimensional equation of parabolic type 
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,dco dco dco dco. ,d co d co d co. . 

p( — + u — + v — + w — ) = ju( — Y + — T + — 2") + /' 
dt dx By dz dx dy dz 

monotonous similar clear scheme has type 

©! +1 -col \u\\ +u" u"-\u"\ I v" I +v" v" - 1 v" I 

yC( + C0- x + fflj + C0 y + ffl + 



+ (0-- + w, ) = \(o r -0=)+ (©,, -<&) + 

2 - 2 z^+c i -^ ^+v ' >' 

r\ UV 

+ , ^ (a n z -a) n z ) + f;,i = \,...,N x -lj = l,..,N-l,k=\,...,N z -\, 
where coefficients by dissipative terms are equal 

+ 4 ^ + 4 ^ +l) ' 
,r r„i J v » l+V » h , lv »'" v ^ 'V +H k ,, , lw »'- w # , v 1/n 

A = L"+ A " «»■ + " V + 4 "* + 4 «* + i)J'G. 

I n i . n i n i h i n \ , n \ n \ n 

I n I . n I « I H I n I . r I « I n 

\U~, \+U, \U-, \-U-i V., +V. V V., -V-i 

,,"" r,, i n ( lk ,k l i '> k ,k u i >J k lk u i 9* lk u \im 

m z =[m+p( — - — K + — - — Km + — - — K + — - — h yi+ i)VQ 

On analytical grid h xM = h xi = h x ,h y]+l = h yj = h y ,h zk+l = K k = K 
these coefficients are simplified 

Q=i+-(-^K+-*-hy+-*-\), M 7=lM+P(-*-h,+-*-\)VQ> 

/.I 2 2 2 2 2 

V?=[M + P(-^h z +-*-h x WQ,M?=[M + P(-l-K+-l-h y )yQ 
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Condition of stable and convergence in form C: 

I n i . n \ n \ n \ n \ , n i«i n 

T n +l 2 K 2 Km 2 K 2h yj+l 

'<!+< '<!-< , 2//7 . 1 1 . 

H — H — ) H — ( 1 ) + 

2/7 2/7 h -\-h h h 

^ ri zk ^ rl zk+l n xi+\ ^ n xi n xi+l n xi 

+ 2^_ 1 + 1 + 2^_ ( 1 + 1 } j a V (/,M,n), 

§3. Monotonous similar schemes 2 order of approximation 

of convective terms in equations with variable coefficient of 

molecular transport 

Consider similar equation with variable coefficient of molecular 
transport // = ju(x,t) > : 

.dco dco. d , dco. r 
p{ + u ) = — (// ) + /, 

dt dx dx dx 
co\ t=Q = d ,co(0,t) = q> x {t),co{a,t) = <p 2 (t),0<x<a,t e[0,T] 
For differentiable function / u = / u(x,t) this equation may be 
transformed to type 

dco , du. dco d co r 

P \-\pu — — ) — = ju — 2" + /, 

dt dx dx dx 

and to build for it monotonous scheme §2. 

In case of disruptive non-differentiable coefficient // = //(x, t) , 
apparently, this transformation is not apply. There is technology of 
building of monotonous schemes 2 order of approximation of 
convection terms to apply to equivalent equation 

,dco dco. d r/ . dco^ d 2 co 

P(— + u — ) = — [O - // min ) — ] + // min — T + / , 
dt dx dx dx dx 

where indicated in // min minimal value of function 
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fi(x,t) > fi mia = const . By technology §1 proved monotonous schemes 
of equivalent equation for k>l: 

I n i . n n n n \ n \ n n 

. „ \U i \+U i 0) i -CD., U i -\U i \(O M -00 i . 

p(co, +— '- - — '- '— + -* '■ - -) = 

2 h xi 2 h xM 

= 4 ^ min u . ^ 

* M i n i . n i n \ n A^min _J 

tt 2// mm 2 2 

« 1 ( <,-a>r <-< 1 }l 

Xi xj+1 xi 

| l ( Kl+A"<l-^" A"+/Cl^"-<l )|y; n t 

fc* 2 fc» + l 2 fc* 

with approximation 2 order and above of convective term. 
For two-dimensional equation of parabolic type 
,da> dm doo^ 8 , doo^ 8 , doo^ 

pt^+^+v^r^^r^— )+—(/*— )+fM^=d,co s =<p 

ot ox oy ox ox oy oy 
monotonous schemes are similar by construction: 

n+l n i n i . n n \ n \ i n i . n n \ n \ 

a> 9 -a> 9 \u,\+u u r \ Uij \ lv ff l+v, Vff-lv ff I , 

/A + 1 % + 1 ^ + 1 ^y + 1 ®y) = 



^n+l 



— (0)" - COl) + (/C« - Anin ) ^ 



= ( Al - AnJ — K - ® " ) + (<„ - Anin ) ~ K ~ «" ) + 
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n . " 



+ — ( -co x -cl>v) + 

h xi 2*2* 

Yi ft ft ft 

ft,. 2 2 

where coefficients by dissipative terms are equal 
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I n i . n i n i n \ n \ . n \ n \ n 

O I U„ I +Ufi I U„ I ~U a I V;; I +V„ I V;; I -V;,. 

* u 4 xi 4 xi+l 4 * 4 yj+l) ' 

lv-1+v- | v »|_ v « 

^ = [a* + *- JL ^ JL K + - JL ^ L v. w & 

For three-dimensional equation of parabolic type 

.dco dco dco dco. d , dco. d , dco. d , dco. . 
p{ — + u — +V — + w — )=—(// — ) +—(m — )+—(M — ) + /> 

dt dx dy dz dx dx dy dy dz dz 

monotonous schemes have type: 

n+l n i n i . n n \ n \ \ n \ . n n \ n \ 

A— + — 2 — x — 2 — — 2 — ' — 2 — ' 

+ — ^ — -on + —i g> z ) = (a - ^,)r^ - °h ) + 
+« -// min )^-K -<)+(//; v - J u m j-*-( (0 n z - ( o n z )+ 

n n n n n n n n 

1 ,Mi+ijk + Mtjk „ Mijk + Mi-\jk „s 1 r Mij+ik + Mi,i „ Mtjk +A;-u „. 
H ( ft),. ft). H ( ft),, (Dt.) + 

K 2 2 *„ 2 ' 2 ' v 

*»■ 2 2 

coefficients by dissipative terms are equal 

I n i . n I " I n \ n i , n i n i » 

= 1 + - J — (— h, + — /?,, , + — h v! + — h M + 

i, A " A A y] A yj 

/"mill ^ ^ ^ ^ 

I n I . n I n i « 

, ' Wg ' +W ijk , , ' gV I -w» 
+ 4 ft* + 4 ^ +l) ' 
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I n I . n I " I n I " I i i I ft I ft 

\W-, \+W: h \W-, \-w.. t \U., 1+U-, \U- h \-U- l 

\u",\+u", \u n ,\-u n , lv" l+v" I v" I -v" 

4 4 4 ■' 4 

Condition of stability of clear scheme 

_£__ | p( '^ l+ ^ h l %lz t V | 'V I+V ^ , 'V 1 -^ , 

Vl 2 ^rf ' 2k xM 2h yj T 2 ^+l 
\n\,n\n\n ^ ^ -. 
H H )+(A* _ Anin) ( 1 ) + 

2/* 2/z h h h 

^"-zk ^ n zk+\ n xi n xi+\ n xi 

+ (My" ~ A»ta ) — (" + — ) + (M7 ~ Anin ) ~ (} + T" ) + 

»yj h yj + l h yj H zk Kk + l h zk 

A ,Mi+ljk + Mijk Miik + A'-l ft >. A / Ay+l*r + A» Mijk + An-U s 
+ ( — + — ' — ) + (— ' — + — ' " ) + 

h xi 2h xi+l 2h xi h yj 2h yj+l 2h yj 

+ j_ ( a; + i+a; + a;+a;- 1 ) ^ a WJXn) 

h zk 2 K + l 2 K 
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Chapter 7. PARADOXES OF DEFORMATIONS THEORY 

AND ELASTICITY THEORY EQUATIONS OF NAVIER- 

CAUCHY-LAME 

§1. Paradoxes of deformations theory 

On page 64 in «Mechanics of continua, V 1» by Sedov L.I., in 
paragraph «On dependency of attendant system basis vector on 
time» it is stated: «... Indeed, in motion of deformed body, distances 
between its points M and M' change. Coordinate curves of attendant 
coordinates system are deformed, and basis ^vectors change in 

time so that their values and angles between them change as 
well...». Let's remind that in the «Basic course of theoretical 
mechanics, part 1» of Buchgolz N.N. for varying vector 

a = aj +a j +a z k , with basis changing in time i,j,k in 

moving body-fixed attendant coordinates Oxyz, time derivative is 
defined as: 

da da r r da _. da 7 - di dj dk 

— = — -i H -j + — -k +a x — + a — + a z — , 

dt dt dt " dt dt dt ' dt 

di _, r dj _ - dk _, r* 
further by Euler formula — = [co, i], — = [<D, j ] , — = [a>, k ] , 

dt dt dt 

where 3 - angular rotational velocity. 

Reducing by dt, we have 

da =dai +dai +da k +adi +adi +a dk (1) 

Therefore, on page 65 in the same section «On dependency of 
attendant system basis vector on time» of the cited book by Sedov 
III expressions 

d? = d£,% and dr'= dE,% (2) 

conflict with formula (1), and, accordingly, must be replaced with 
corresponding to formula (1) expressions 

cff = d%% + gd§ { , dr'= d£% + ?d3\ (3) 

with all resulting consequences in the deformations theory. Strictly 
speaking, doubtfulness of formula (2) from the point of view of 
formula (1) is confirmed by statement that «... and basis $\ vectors 
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are changing in time so, that their values and angles between them 
change... » because 3/ , i= 1,2,3 - basis vectors at the initial moment 

of time t , p ( , i= 1,2,3 - basis vectors at the current moment of time 

t, and also r and r' and, as already said, variables. 
According to Sedov III for representations (2) it is assumed 

\dr\=ds = J$^d£ I >&»=$■§> 

I dr' \=ds'= Jg~d?d?, g' tj = 3; ■ 3) (4) 

obviously, for correct representations (3) instead of (4) obtained are 
absolutely other expressions, since 

I df 1= (dE,% + ?<$ , dE,% + gdf) 2 , 

i 
I dr'\= (d?3; + gds: , d?3; + gd5: Y (5) 

Further on, based on (4) derived is formula III: 

£ tJ =\ ( y^ ] +i7 A -v>,v.>vj (6) 

Product of derivatives V-hvV-w in (6) is commonly supposed to be 

a small value and it is neglected, thus (6) is simplified down to the 
following known expression of motion tensor 

1 dw f dw. 
e„ = -(— + ^- L ), i, j = 12,3 , (7) 

" 2 dx, dxj 

If Vj-W^V w is not small, (7) does not take place. Obviously, due to 

(5) instead of (6) and (7) there will be absolutely different relations. 
(From the point of view of common sense, formula (6) if absurd, if 
Wj and Xj would have different dimensionalities, for example, [ vv y ] 

=m/c, and [ X- ]=m). The same circumstance is present in the book by 

G. Mase «Theory and problems of continua mechanics)), if to bring 
Mase's symbols into compliance with symbols of Sedov: dx=dr , 

dX=dr\ £ =x v ^ 2 =x 2 , ^ =x 3 , s tj =E. j ,w i =u t , w=u etc. 

Formula (6) in 111 has the following form 
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1 du i duj du k du, 
E n =-( — '- + — — -) 

2 dx ; dx : dx : dx : 



(8) 



j i i j 

Thus, remarks (2) and (3), (5) are related to derivation of formula 
(8). It is considered that in (8) gradients are small in comparison with 
one and they can be neglected. Strictly speaking, and it is evident, 

k d u k d u t 
that dropping of gradients V-w fe V w and is connected 

1 dx, dx, 

with slanting (7) to false symmetry of continuum strains tensor. 
Continuum strains tensor in the general case is not symmetric (refer 
to Chapter 1). 

§2. Alternative presentation of relative movement du 

In symbols of G.Mase /2/ relative movement vector du (in the 
book by Sedov L.I.Ill - denoted as dw) is introduced as difference of 
vectos of oint Q Q initial position and final position of the same point 



P a : du=u 



(So). u W 



, W 



Expansion of u in Taylor series into 



neighborhood of point P Q gives du=Kdx, or in componentwise 
representation (inexact differential) 



du, , 
du. = dx. , j'=l,2,3 

dx i 

(by index j summing up from 1 to 3 is performed). 

Mark that tensor K is nonsymmetric: 



(1) 



K= 



<5i*~i <5i*~ x c2i*~ x 
c?je r c?jc^ c>je 3 

G>M,~, G>1*~~, <3z^.- 



~-\ 



c?jc x iDje^. iDje^ 
C>14.^, c?i^^. c?i*.^. 



c?jCi c)je^ ctje^ 



J 



(2) 



In /1/-/3/ expansion (1) is traditionally reduced to the form 
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A du. duj l Qu. du, 
du i =[-(—^ + — L ) + -(—^ L )]dxj, (3) 

2 dx ; dx ; 2 dx ; dx t 

J ' j l 

or in symbols of Sedov III 

Jw. = [- ( — '- + — -) + - ( — '- L )]dx . 

2 5x ; ck,- 2 3x y 3x; 

First sum in (3) is represented by Euler tensor of linear 
deformation 



r 



E = 



9«, 1 „ Qu, <9i/ 7 _ 1 . 9m, 

L _ ( L _, ?_) _ ( L 

Sxrj 2 Sx 2 Qx x 2 cbc 3 



<3w 



\ 



1 9i/ x 

2 9-3C 2 

Qu x 



Qu 2 Qu 2 1 



) 



ebc 2 2 



2 Sx, 



9x x 

<3w 3 1 Qu 3 

Dx 1 2 9x 2 



Qu. 

QXr, 



Qx x 
Qu 



3 ) 

2 \ 



3x 



<9w. 



3 

Qu, 



Qxr, Qx 



(4) 



y ^ u^ 3 tv^v- 1 ^ «--'- A -2 tv-A.3 t^-A.3 y 

This is a symmetric matrix; second part in (3) are components of 
Euler linear rotation tensor 



1 du 

®s =-( 



du 



L ) 



2 dx , ex 
At the same time they are components of antisymmetric tensor 



(5) 



r 



R= 



V 



2 cbc 2 

1 ,Qu^ 

— ( - 

2 9xj 

1 Du 3 

2 cbq 



oV/., . 1 , C?M, 

— ) — ( - 

Qx 1 2 dx 3 



2 <9x 3 



Qx 2 
Qu x 

9x. 2 Qx~ 



Du 3 

Qx x 
Du 3 

Qx~, 



\ 



etc, 



)o 



y 



Obviously Taylor series is presented in two ways in the form (1) 
and (3): 
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du=Kdx h <iu=E<ixH — [rotu,Jx] 
or with the help of antisymmetric matrix as follows 

du=Edx+ Rdx] 



(6) 



(7) 



In addition to (6), (7) there is the infinite number of various forms 
of representing Taylor series (1). With this purpose we'll introduce a 
family of one-parametric matrices 



c = 



<9wj 

dx l 
f 



b — \ du, 1 du. 



v b dx t 



■ + - 



b dx 



2 J 



b — \ du, 1 du. 



v b dx l 



b dx. 



3 J 



b — \ du. 1 du 1 



b dx, 



y b dx 2 



- + - 



i J 



du 2 

3x, 



b — \ du, 1 du n 



b dx. 



v b dx 2 



■ + - 



3 y 



Z? — 1 du. 1 3m, 



Z? 5x 



v v b dx 3 



■ + - 



i J 



b — \ duj 1 du^ 



b dx, b dx. 



2 J 



du 3 

dx 



3 J 



(8) 



Taylor series (1) with the use of matrices (8) has the infinite number 
of alternative or universal representations 



du = C, dx+ 



b-\ 



[rotu,dx] , b ^ 0, 1 b \< co 



(9) 



When b=l, it comes out that C , =K, when b=2 - C 2 =E etc. 



§3. Paradoxes of Navier-Cauchy-Lame hypothesis 

For applying Hooke law in the theory of deformable solids 
hypothesis of Navier-Cauchy-Lame /!/ is used, stating that in Taylor 
series (or inexact differential) 



1 du t du i du, duj 

dU, = [- ( L + -) + - ( L ~)]dX : 

2 dx, dx, 2 dx, dx, 

J ' J l 



(1) 



for defining components of strain tensor first half of this series is 
enough 
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1 du du : 
du i =-(—*- + -±)dx 1 , (2) 

2 a*, ax,. J 

second part of the series is neglected, assuming that 

1 du du j 

-(-L--2.) = 0, (3) 

2 or, a*,- 

although both these expressions consist of the same gradients. 

As a result, according to Navier-Cauchy-Lame hypothesis, Hooke 
law was defined and is widely used in the following form 

1 du. dU: 
tt m = AS }J divu + 2/je M , e fi =-(—^ + —±) ,i,j = 1,2,3 , (4) 

[o,/* j, 

where S..=< , A, jU—Lame coefficients, ;r, 7 = ^». 

[1, / = j, 

Thus, according to (4) it is stated that forces deforming a body 
create only linear deformation Edx. According to Navier-Cauchy- 
Lame hypothesis Euler linear rotation is equated to zero 

1 

— [rotu,<ix]=Rcix=0, which is specifically expressed because of (3) 

as equality of displacement curl to zero (see Lurie /3/): 

rotu=0 (5) 

(in particular, solvability condition of A.N.Konovalov has the form 

/A/ [rotudM = 0) 

D 

Accordingly, formula (1) is cut and takes a shortened form (2), far 
from Taylor series: 

du=Edx (6) 

Of course, the indicated paradox is connected with slanting (4) to 
formula of Euler tensor of finite deformations 

E = l(fy , dH J du k du k ^ 

'' 2 dx, dx. dx. dx. 

J ' ' 3 
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in which, according to small deformations theory products of 

du k du k 

, that can be quite large are dropped, and suggested is 

dx, dx, 

1 J 

1 du. du, 
doubtful formula E =-( — L + — L ) . 

2 dx, dx i 

J l 

So, if to follow Navier-Cauchy-Lame formula, according to which 
rotu=0, than from representation of Taylor series in universal form 

b-\ 

du = Cj dx+ [rotu,dx], b ^ 0,1 b \< qo (7) 

b 

dropped is the second addend and, based on expression 

du = Cj dx, (8) 

logically, Hooke law must be defined in the form 

7c n = X8 n divu + bus, ,s H = + , b * 0,1 b \< oo (9) 

11 v ^ ll J ' bdxj b dx t 

In formulae (4) Navier-Cauchy-Lame strains tensor is symmetric, in 
formulae (7) there are infinite many nonsymmetric strain tensors 
when b ^ 2 . Formula (7) is the Taylor series with any values of 
b ^ , including for b = 1 . Hooke law (9) was formulated for 
representation (8), obtained from (7) when rotu=0, and for b = 1 
takes the following form 

OH 

n .. = A5 tl divu + US ..,£.. = — L , i, 7=1,2,3, (10) 

1 ' J ] dx. 

And since when b = 1 , C , =K, then Hooke law (10) corresponds to 

full Taylor series du=Kdx\ 

The main paradox of Navier-Cauchy-Lame hypothsis is as 
follows. According to this hypothesis there are equivalent equalities 
(3) or (5), which result in equality to zero of curl components 

da, d u i n . . , ,,„ du, du- 

— ~ = 0,i, 7=1,2,3 or — '- = — -,/,j=l,2,3 (11) 

OX; OX; OX; OX; 

J ' J l 

Strains (4), built by Navier-Cauchy-Lame hypothesis, can be 
represented as follows: 
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7T H = 18-divu + 2U£ H , s .. = - ( — L + — ^) = — - — ( — '- L ), 

11 y h J ' J ' 2 fa j 8k/ dxj 2 8k j 8k/ 

which, due to equalities (11) are transformed into formulae (10): 

on 
n H = XS n divu + jus n , e H =— L , i, j= 1 ,2,3 (12) 

dXj 

and are symmetric due to second equalities (11). 

Again result is Hooke law in the form of (10)! Main paradox is 
also in the fact that solutions of Navier-Cauchy-Lame equations 

d 2 u 
p — — = p F + (X + u)graddivu + uAu (13) 

dr 

must satisfy equalities (11). Let's study this issue. With this purpose 
we'll apply operation rot to equation (13): 

d 2 rotu - 

p — = p rotF + juArotu , (14) 

dt 

since rotgraddivii = . 

Theorem 1. In nonstationary problems of elasticity theory rotu ^ , 
accordingly, Navier-Cauchy-Lame hypothesis (4) is not correct, 
strain tensor may not be symmetric. 

It is obvious, that for wtF ^0 equation (14) has a nonzero 
solution rotu ^ 0, therefore Navier-Cauchy-Lame hypothesis is not 
correct, which means that strain tensor needs to be defined by 

formula (10) with nonsymmetric strain tensor. Let rotF = 0. Then 
equation (14) will take the following form 

d 2 rotu 

~d? 

This wave equation has general nonzero solution 

rotu = (sin 7vc l sin 7vc 2 sin 7ix 3 sin y/SjUM)^ + 

+ (sin 7vc l sin 7vc 2 sin tzx 3 sin j3juM)i 2 + (16) 



p — — ^ — = fjhrotu (15) 



+ (sin 7vc x sin tvc 2 sin 7vc 3 sin -j3jUM)i 3 ^ 
accordingly, Navier-Cauchy-Lame hypothesis is incorrect again! 
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Let's consider equation of elastic equilibrium, obtained by Navier- 
Cauchy-Lame hypothesis: 

p F + (A, + ju)graddivu+ ju/Su = (17) 

Curl operation (17) gives an elliptic equation relative to rotu : 

p G rotF + /jArotu =0 (18) 

Theorem 2. If wtF ^ , equation (18) has a nonzero solution 
rotu ■£ , accordingly, in problems of elastic equilibrium, 
hypothesis of Navier-Cauchy-Lame is not correct, strain tensor 
may not be symmetric. 
The proof is straightforward. 

Now let bulk forces are such that rotF = . Then equation (18) 
will transform into homogenous elliptical equation 

Arotu = , 
which has a zero solution 

rotii=0 (19) 

Theorem 3. In problems of elastic equilibrium, hypothesis of 

Navier-Cauchy-Lame is correct, if rotF = , i.e. rotu =0, 
accordingly, strain tensor may be symmetrical. Displacement 
vector has a potential. 

Proof is already given in the form (19). Potential is introduced by 
formula 

u = grad<& , (20) 

because rotu = rotgrad® = . And in this case strain tensor is 

defined as (10), symmetry takes place because of (19), which results 

in equalities (11). 

So, according to theorems 1,2,3 strain tensor must have the 

ou 
following form (10): n n = AS ij divu + fis- t ,e H = — '-, i, 7=1,2,3. 

dxj 

According to Navier-Cauchy-Lame hypothesis, forces acting upon 
a body cause only linear deformation Edx, in the result a distorted 
expression is obtained, as said above, 

1 du. du: 
du : = — ( — '- H ) dx ,. , du=Edx, 



2 dxj dx, 
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having nothing in common with Taylor series (1) du=Kdx. Rejection 

of Navier-Cauchy-Lame hypothesis means that if to base on formula 

of Taylor series 

1 
du=Edx+ — [rotu,dx] , 

it is necessary to consider also forces of strains causing Euler linear 
rotation. With this purpose let's denote strains proportional to Euler 

linear deformations by jr.. -, strains proportional to Euler linear 

** 
rotations by 7T- : 

/. a iL+ & t n - E^Jjh.) (21) 
" 2 dx, ex, ' 2 &, at, 

J l J l 

By Hooke law total force will be 

7t ■■ = Xdjdivw + 7V „ + /T*. 

j i y J 1 j i 

Following substituting n- and n- obtained is 

(j\A 

7T n = ASjjdivil + ju — '- , i, 7=1,2,3, (22) 

dx. 

where <5- =•< , A,jU— Lame coefficients. 

[!»' = 7. 

Strain tensor in (22) is nonsymmetric 7r„ ^ 7T jn i ^ j , because 

here equality to zero (11) component of vector rotu is not suggested 
initially. So, strain tensor (22) again matches tensors (12) and (10). 

Nonsymmetry of strain tensor in continuum dynamics is proved in 
Chapter 1. Rigorous proof given there is easily transferred to 
deformed bodies in theory of elasticity as well. In monograph of 
Lurie A.I. /3/ et al. symmetry of strain tensor in theory of elasticity is 
established for the state of elastic equilibrium 

d7T C)7T; firr -, 

^ + ^L + ^ + pF=0 ( 23) 

dx dy dz 

Similarly to Sedov and Loytsynasky, considered is moment of forces, 

acting upon body in volume r with surface o , 
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t a 

from which, due to arbitrariness of r derived is relationship 

dx dy dz 

wherefrom, based on equation (23) Lurie et al. make a conclusion of 
symmetry of strain tensor in elastic medium. Criticism of such 
deductive method was given in §8 Chapter 1. Therefore, it shall be 
assumed, based on the above-mentioned facts, that tangential strains 
n ij ^ n ' ii 'J^i ( m Cartesian coordinates n ^ 7T yx , 7T zy ^ n , 
n xz ^ 71 , ) in the general case are not equal to each other. Equality 
of tangential strains - symmetry of strain tensor - at any point of 

du. duj . . 
continuum, i.e. execution of equalities (11) = ,i, 7=1,2,3, 

dxj dx j 

just means the fact that according to results of §8 Chapter 1, at the 
given point, moment of major (resultant) force matched with the 
major moment of forces, i.e. with resultant of moments of forces. 

Thus, it is natural that there is a long-lasting need of revising 
equations of dynamics of deformed solids, in other words, equations 
of Navier-Cauchy-Lame. 

§4. Paradoxes of Hooke law for symmetric strain tensor Navier- 
Cauchy-Lame 

Let's, for example, consider equations of elastic equilibrium (17) 

§3 for a two-dimensional case, when force F = Oij + 0z 2 + F 3 i 3 is 

perpendicular to plane (x 1 ,x 2 ), and will represent (17) in 

projections 

d 
(A + ju) divu + /jAu k = 0, k = 1,2,3 , (24) 

dx k 

or in more details for a two-dimensional problem 

, d .du, 3w 2 , ^ d ,du k du t . nj 10 ,„„ 

A ( — - + — -) + ju\ — ( — - + — L ) = 0,k = l,2 (25) 

dx k dx 1 dx 2 t^i dx t dx i dx k 
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For infinite number of displacements, i.e. solutions of system of 
equations (25), symmetric tangential strains of Navier-Cauchy- 
Lame hypothesis have zero values, i.e. in all points of the body are 
zero 

du. du j 

Ha = TTji = M^r + -^) = o,ij = 1,2, 

while normal strains are different from zero n u = 2jus u ^ . 

Let's limit to giving a small list of displacements, in which this 
fact is actual (for shortness, denoted as 
u x =u,u 2 = v, x x = x, x 2 = y,u 3 =w,x 3 = z) : 

1) u=F(sin kiX cos kiy ), v=F(-cosk\X sink\y ), 

2) u= U(sink,2X cosk,2y -coskjX sink,2y ) , 
v= U(sin kjX cos A: 2 y -cos k.2% sin k^ y ), 

3) u=W(-cos k^xsin £3 y ),v=W(sin k^xcos k^y ), 

4) u = Q(sink^X sink^y ),v=Q(cosk^X cosk^y), 

5) u = T(sink^Xsink^y ),v=T(cosk$X cosk^y ), (26) 

6) u=M(sin kfrX sin k^ y +cos k^ y cos k^X ), 
v=M(sink( ) X sink^y +cosk^y cosk^X ), 
7)u = S(e kl{x+y) ), v = -S(e kl(x+y) ), 

for three-dimensional displacements when there are no bulk forces 

8) M =D((/ 8;y -/ 8Z )/ 8X ),v=D((/ 8Z -/ 8X )/ 83, ), (27) 

w=D((e 8 *- e 8;y )e 8Z ), where coefficients k t = const, 

i = 1,2,3,4,5,6,7,8, are selected arbitrary from infinite interval 

-co<ki < +00. Standing here differentiable functions F,U,W,Q,T,M, 

S,D are also chosen arbitrary. Obviously, from the given list it is 
possible to form any new linear combinations of u=F+U, v=F+U 
etc. type. Fields 1-7 in (26) correspond to plane displacements, and 
are solutions of equations (25), since they transform divergence into 
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du dv 
zero 1 = 0, fields of displacement 8) in (27) turn three- 

dx dy 

,..,,. . du dv dw n 
dimensional divergence into zero 1 1 = U , symmetric 

dx dy dz 

tangential strains are equal to zero 

du- dUj 

7T n =ju( — L -\ ) = 0,i^ j,\/(i,j) , normal strains are nonzero 

dx, dx, 

J l 

7t u = 2jU£ u ^ 0, i = 1,2,3 . It is necessary to note that partial 

solutions (26), (27) of elastic equilibrium equations (24), are also 
partial solutions for equations with nonsymmetric strains tensor 

d . -, 
A divu + /uAu k =0,k =1,2,3 , but here, nonsymmetric 

8x k 

• , • dw. „ . . 

tangential strains are nonzero n = u — - ^U, n •• ^ n ,i ^ / , 

which is physical enough. 



§5. Equations of elasticity theory for nonsymmetric strain tensor 

Elasticity theory equations in the form of Navier-Cauchy-Lame for 
symmetric strain tensor are given in /1/-/3/: 

d 2 u 
p Q — t- = p Q F + (A, + u)graddivu + uAu , (1) 

dt 

u = ui +vj +wk - displacement vector, u x =u,u 2 =v,u 3 = w . 

Above it was justified, that if forces acting upon particles of 
deformed body move them to 

du=Kdx, 
then appropriate components of strain tensor must, in accordance 
with Hooke law, be proportional to components of nonsymmetric 
displacement tensor K. 

Hooke law for nonsymmetric strain tensor was given in §2 

syjj 

n jt = ASydivu + jus jt , s fi =— L ,i,j = 1,2,3 (2) 



dxj 
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Dynamics equations are result from substituting (2) into equations 

and for nonsymmetric strain tensor (2) it is derived 

p — - = p F + /Igraddivu + juAit (3) 

dt 

Projections (3) on coordinate axis form a system of three scalar 
hyperbolic -type equations 

d 2 u „ „ d ,du dv 8w. ,d 2 u d 2 u d 2 u. 

Pa — 7 = Po F x +A — ( + — + ) + M — 7 + — 7 + — t)> 

dt dx dx dy dz ' dx dy dz 

d 2 v „ . d ,du dv dw s ,d 2 v d 2 v 3V 
p ^r = p F+A — ( — + — + — )+ ju( — - + — - + — -), (4) 
y °dt 2 yo y dy dx dy dz dx 2 dy 2 dz 2 

d 2 w „ „ d ,du dv dw^ ,d 2 w d 2 w 3 2 h\ 

Po — t = Po F z +A — ( — + — + — ) + ju( — t + — 7 + — r) 

dt dz dx dy dz dx dy dz 

Further on given are numerical methods of solving problems of 
elasticity theory in displacements. 

Note. Consistency conditions of components of nonsymmetric 
displacement tensor are 

5% , d\ v ] d% | d 2 s JM 



dXjdx M dx^dxj dx t dx v dx v dx t 
d 2 £ Mi d 2 c iM d 2 e xj d 2 s jv =q 



dxjdx v dx v dxj dx t dx v dx v dx t 

§6. Explicit scheme 

Let's state the idea of building such a scheme in the example of 
two-dimensional initially-boundary problem 

d 2 u d ,du dv d 2 u d 2 u 

Po^y = P F x +U + M ) — (— + —) + p(— + —), 

dt ox ox dy dx dy 
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d 2 v 0n d .du dv. / d 2 v d 2 v 

i 3h i 

M Lr„ = d (*,)>), t-(=i =d(x,y),u \ s =p(x,y,t), 

10 ^ I 

Here, parameter // = // for symmetric strain tensor and // =0 for 
nonsymmetric strain tensor. The scheme has the following form 

M,.. -2m,.. +u ti „ 1 

Po- ^^=p F x+ U+ju°){—(u:-u;) + 

- n n n n , - 

1 , V i+lj+l ~ V i+lj-l V i-l;+l _V !-l;-lM , r /"„« „n\, (,.<> „»^ 



v "+ 1 _2v"+v"- 1 n 1 

a- 5 — 5 ! - ! -=Af,+aV){ r (v;-^)+ 

w 

+ ( — —)}+/A — (v x -v i )+ — (v -vj], 

i=i,...,N x -lj=l..,N y -i 

Initial conditions are specified 
Q h =^i = l,...,N x -l;y r j = l,...,N y -l: 

U =d n ,V =U'n-M- =dn +nl-,V- =a' n +T&'-, 
ij Oij ' ij Oij ' ij Oij y ' y Oij ij ' 

and boundary conditions of the first kind 

s h = i^o,i=N x ,j=W y ;j=o,j=N y ,i=W x ^i h =n h [js h , 

x = 0,< = <,< = $/, j = OAT 

x = a 1 ,u Nj =(p NxJ ,v NxJ =^ x .,j=0,N y , 
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errors of approximation are 

y% =0(T 2 ) + 0(h 2 xi ) + 0(h xM -hj + 0(h 2 yj ) + 0(h yj+l -h yj ), 

¥ n vi] =0(T 2 ) + 0(h 2 xi ) + 0(h xM -h xi ) + 0(h 2 yj ) + 0(h yj+l -h yj ), 

i = l,...,N x -l,j=l,...,N y -l; 

Condition of stability is the Courant criterion 

r 2 /4 ( + — ) + ( + — )] ^ 1, 

h X M+Ki Km K h yj+1 +h yj h yj+1 h yj 

i=l,...,N x -l,j=l,...,N y -l. 

Note. Combined coefficient /1 + jU is included with a special 

purpose of checking symmetry conditions of strain tensor 

du. du. 

= ,i, 7=1,2,3 for // = // in numerical computation of 

dx i dx: 

specific problems. 

§7. Solution of three-dimensional problems of elasticity theory on 

a uniform mesh 

For three-dimensional equations 

d 2 u „ „ dp / d 2 u d 2 u d 2 u s 

Pa — y = Po F x +A — + M — 2" + — 1 + — 2")' 
dt dx dx dy dz 

i2 a a2,. a 2.. a2. 



3 v „ .dp ,d v d v d v. 

— T = Po F v +A— + jU( + + 

dt 2 y dy dx 2 dy 2 dz 2 



d 2 w r, . dp / d 2 w d 2 w 3 2 w\ 
dt dz dx dy dz 
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du dv dw 

— + — H = p 

dx dy dz 



specified are meshes, in their nodes initial conditions 

Q h = ^i = l,...,N x -l;yj,j = l,...,N y -l;z m ,m = l,...,N z - 

U ijm ~ ^Oijm ' V ipn ~ ®Qijm ' W ijm ~ ®Qijm ' U ijm ~ ""Oijm ~"~ ^ ijm ' 

v 1 . = a / a + m / .. , w 1 . =au» +TC& , 

ijm ui/m *^ym ' ym Qijm * ym ' 

i = 1,^-lj = l,N y -l,m = IN Z -1 
and boundary conditions of the first kind 

S h ={i = 0,i = N x ,j=0J^„m = 0M~ z ;j = 0,j = N y , 

n h =n h {js h ; 

r\ n+1 /n n +l n+ l jG n+ ^ n+ ^ Jf? 1 ^ 

x — \j, u 0jm — (p Qjm , v 0jm — % ; . m , w 0jm — <p§ jm , 



T — n u n+1 - m n+1 v n+1 — /£" +1 w" +1 - /tt" +1 
A Wj , "^j m V^ym ' V NJm rNJm ' ^JVJm S S V e ;m ' 

~ — — — — — — — rv n+1 ,.«+l n+1 £n+l n+1 fcn+1 

]=Q,N y jn = 0, ^ , y = U, W, m = Worn ' V ;0m = #0m > W ,0m = $0m . 

_y — u 2 ,u iN — (f> m , v iN — <Pft„, w iN — <Pft„, 

y y y > > ? 

i = 0,N x ,m = 0,N z , 

On+l n+1 «+l ,C«+1 n+1 fera+1 

» M yO = <PijO » V yO = <%0 » W yO = %0 » 

~ ,-. ,."+1 ^."+1 , ,«+! >C" +1 -...,"+1 £"+1 

z=a 3 ,u ijNz =<p ijNz ,v ijNz =<Pi jNz ,w ijNz =ty Nz , 

i = oaT, j = oaT, 

global iterative process is built as follows: 

2 n+l,k n+\,k 

n+\,k r?n . t o ri+ijm "i-ljm 



um uwn ,-,7 

M> 2h 
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2 n+l,k n+\,k 

n+l,k r<n , T n rij+im t ij-lm 

V.. =P .. -\ /t 

ym vijm r*] 



Po 



y 



n+\,k 



n+Lk 



,,,n+l,fc r?n T n "ijm+l Pijm-1 
W,™ = r.. -\ A — 



ijm 



wijm 



Po 



2h. 



l<i<N x -l,l<j<N y -l,l<m<N z -l, 

n+l,i+l n+l,k n+l,k n+l,k n+l,k n+l,k 

a A /y. rijm rijm M im i-\ im ij+lm ij-lm 
+ AJS) — — = — — + — + 



Ijm 



e 



IK 



2h 



n+l,k n+\,k 

W —fit 

ijm+l ijm-l „+i £ 

H P 

r 11 



2h_ 



ijm 



1 < i < N x - 1,1 < j < N y - 1,1 < m < N z - 1 , 

n+l,k+l 

0+4, j r"-* ~ - — + 



rN x jm rN x jm 

1 V , \m U ' ^ 



h. 



d$ 






,n+l 



d$ 



\H+1 



n+l,k 



~dv w,7ffl ~&~ NJm m ' 



n+l,£+l n+l,k n+\,k sn n+l 

(1+4 #) oym ~ "o>" _ M i>" ~ tPojm 



6 



K 






+(^;+ec-Kr^7^-u<m<Ar z -i, 



a+w) 



n+l.i+l n+l,£ 

riN,,m riN.m 



e 



—(% 



2n+\ 



n+l,k 



n+l 



O < JV V 



+ ■ 



*2V.,m iN ,,-\m 



h„ 



+ 
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-d<fe nH 



./~Y~\n+\ __ n+\,k 
+ \~^~)iN„m PiN„m ' 



e ^dx"" m h 



n+l,k+l n+\,k pj «+l / C'" +1 '* : 
(-1 , A a\PiQm ~ PiOm /°^ \n+l , V ilm — ' ViOm 
U + Aom^ 7 = C— J.Om + \ + 



y 



oz 



n+l,k+l n+l,k _ fti 

PijN z -PijN z =( d<P t gfr 

V^ a 



+1 n+l,£ 

W,, 



n i A A\ * 'J"z * V" z /"f\n+L ''JN Z ijN z -l 

(i + V #) ^ = tV, + ; + 



ay ■ < 



n+l,k+l n+\,k _ n+1 fcn+l,£ 



P#o' -P#o' ,3#\„ 



+1 , W ij\ V^O 






H+l «+!,& 



+ (^);-PoM<^<^-i,i<;<iV v -i, 

oy 

where 

J p n Ih 2h 2 

ru X X 

1 l + sign(i-15), 1 r 1 l-sign(j-N +1.5) 

+ ^— ] + [ + 

2/i 2 2/?„ 2/r 2 

x y y 

1 l + «>n(/-1.5) n 1 r 1 l-sign{m-N +1.5) 

+ — - -] + [ ^^ ? - + 

2A y 2 2h z 2h z 2 

1 l + 5i#n(w-1.5),. 

+ — ]}, 

2h z 2 

i = l,...,N x -l,j = l,...,N y -l,m = l,...,N z -1; 
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1 jl l A 1 f 1 l 

j=\,...,N y -l,m = l,...,N z -l; 

2 1 2 i 

A*,- =i +— A—^,A iQm = i +— ;t-— , 

A> 2h 2 y A, 2h 2 y 

i = l,...,N x -l,m = l,...,N z -l; 

2 1 2 i 

T „ 1 „ T „ 1 



A =1 +— A— ^,A, =1 + — A- , 

Parameter > is selected from interval < 6 < 1/2. 

Process of calculations is stopped on such number of iterations k , 
with which the following inequality is true 

n+\,k"+\ n+\,k" 

maxl (l + A ijm 0)^ l£*!—te ef ff «0,ff*0 

0<i<N x ,0<i<N y ,0<m<N z 6 

While observing this criterion the last approximations 

n+\,k* n+l,k' n+\,k* «+U* 
Pirn ^iim >^iim ^ijm ' are taken as solutions 

n+l,k n+1 n+\,k n+1 n+l,k n+1 n+l,k n+1 

p « D ,U • ~U- ,V» «V- ,W- « W- , 

^ym *y'« ' ym ym ' ym ym ' ym ym ' 

/ = 0,...,^,7 = 0,...,^,m = 0,...,iV z . 
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